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ABSTRACT 


This  report  describes  In  detail  the  technical  progress  made  In 
studies  of  strong  Infrared  light  scattering  from  coherent  spin  wav 
during  the  second  six  months  of  a  three  year  Interdisciplinary  research 
program^ titled:  "Steerable  Volume  and  Surface  Spin  Waves  In  Ferrl- 
magnetlc  Films".  One  major  goal  of  the  program  Is  the  growth  of 
yttrlum-lron  garnet  (YIG),  single  crystals  of  high  quality  In  both  thin 
film  and  bulk  form.  Another  Is  the  development  of  novel  techniques  for 
controlling  and  studying  energy  propagation  characteristics  of  volume 
and  surface  spin  waves.  One  promising  technique  for  accomplishing  the 
latter  Involves  the  use  of  optical  techniques  to  probe  the  dynamics  of 
spin  wave  propagation.  Since  YIG  Is  opaque  to  visible  light  but  trans¬ 
parent  In  the  Infrared,  this  form  of  radiation  has  been  employed  In 
these  experiments. 

Coherent  light  scattering  from  coherent  microwave  spin  waves  In 
yttrium  Iron  garnet  (YIG)  Is  studied  theoretically  and  experimentally. 
The  method  Is  then  utilized,  In  conjunction  with  standard  microwave 
techniques,  to  probe  the  propagation  of  spin  waves  and  magnetoelastic 
waves.  In  spatlally-and/or-temporally-varylng  magnetic  fields. 

The  dispersion  of  magnetoelastic  waves  propagating  parallel  to 
saturation  magnetization  when  the  latter  Is  In  an  arbitrary  crystallo¬ 
graphic  direction,  first  presented  by  Morgenthaler ,  Is  rederlved.  For 
this  Important  case,  spin  wave  defocusslng  Is  reduced  and,  when  properly 
employed,  allows  variable  delay  of  magnetically  Injected  spln/longltud- 
Inal  elastic  waves,  and  conversion  from  longitudinal  elastic  waves  to 
spin  waves  In  spatially  varying  fields  to  be  studied.  Frequency  shifts, 
and  pulse  compression  of  longitudinal  elastlc/spln  waves  subjected  to 


pulsed  magnetic  fields  are  also  reported,  as  are  some  preliminary  exper¬ 
iments  on  the  spatial  conversion  of  shear/longltudlnal  elastic  waves  via 
spin  wave  coupling. 

■*v 

We  also  report  the  first  direct  observation  of  strong  Bragg  scattered 
Infrared  light  (1150  nm  wavelength)  from  coherent  microwave  spin  waves. 

The  latter  were  generated  In  YIG  through  efficient  spatial  conversion  of 
Injected  longitudinal  elastic  waves.  The  maximum  scattered  light  Inten¬ 
sity  Is  found  to  be  at  least  five  times  stronger  than  that  of  the  longi¬ 
tudinal  elastic  waves  of  comparable  power.  In  contrast  to  elastic  wave 
scattering,  spin  wave  light  scattering  Is  unique  In  that  the  diffracted 
light  Intensity  Is  observed  to  depend  on  the  Incident  light  polarization 
and  whether  the  diffracted  light  Is  upshlfted  or  downshifted  In  frequency. 
These  characteristics  are  In  complete  agreement  with  a  revised  theory 
that  Is  found  to  be  necessary  after  a  critical  review  of  previous  work. 

The  revision  Is  shown  to  be  capable  of  Interpreting  adequately  not  only 
our  experiments  but  also  observat-'ons  that  previous  authors  have  reported 
as  "anomalous." 


I  INTRODUCTION 


The  Microwave  and  Quantum  Magnetics  Group  and  the  Crystal  Physics 
Laboratory,  both  within  the  M. I.T.  Center  for  Materials  Science  and 
Engineering,  have  undertaken  a  three  year  Interdisciplinary  program  In 
the  area  of  "microwave  magneto-ultrasonics"  aimed  at  further  developing 
several  novel  concepts  which  may  lead  to  new  and/or  Improved  solid 
state  devices  employing  electromagnetlc/spln/elastlc  wave  coupling. 

Device  possibilities  Include  multl-tapped  delay  lines,  magnetoelastic 
beam  switches  and  pulse  compression  filters.  In  particular,  the  research 
program  will  concentrate  on  the  growth  and  exploitation  of  Improved 
single  crystal  yttrium  Iron  garnet  substrates  In  which  volume  and  surface 
spin  wave  propagating  at  microwave  frequencies  can  be  magnetically  steered 
and/or  otherwise  controlled. 

In  order  to  produce  crystals  of  yttrium  Iron  garnet  of  the  quality 
needed  for  the  research  envisioned  In  this  program,  two  major  constraints 
must  be  dealt  with.  First,  the  crystals  must  be  of  the  necessary  high 
quality  to  avoid  Introducing  extraneous  effects  due  to  grain  boundaries 
and  strain  Inhomogeneltles.  The  chemical  purity  is  of  major  Importance 
since  rare  earth  Impurities  drastically  increase  spin  wave  relaxation 
rates.  Second,  the  growth  technique  should  permit  ready  fabrication  of 
the  types  of  structures  desired,  such  as  multi-layer  configurations  and 
thin  slabs.  Consequently,  both  chemical  vapor  deposition  (CVD)  tech¬ 
niques  and  top  seeded  solution  (TSS)  growth  appear  attractive  and  both 
will  be  employed  throughout  the  entire  contract  period.  Thin  films  and 
bulk  crystal  substrates  will  be  evaluated  optically,  magnetically  and 
acoustically  for  purity,  homogeneity  and  low  losses.  Promising  specimens 
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will  be  used  for  magnetostatic  wave  and  exchange  dominated  spin  wave 
propagation  experiments  at  microwave  frequencies. 

A  surface  spin  wave  propagating  parallel  to  an  air-crystal  boundary 
of  a  thin  film  has  Its  energy  largely  confined  within  some  characteristic 
distance  that  under  certain  conditions  should  be  magnetically  controllable. 
In  one  Instance  causing  the  energy  to  be  closely  confined  to  the  surface 
and  In  another  allowing  It  to  spread  out  and  fill  the  entire  film.  Vol¬ 
ume  spin  waves  are  also  subject  to  magnetic  control.  For  example,  beam 
steering  Is  possible  due  to  dlpole-dlpole  Interactions  that  In  turn  can  be 
Influenced  through  directional  changes  In  the  magnetizing  field. 

A  high  power  (50  mw)  He-Ne  laser  operating  continuously  at  l.lSp 
has  been  obtained  for  use  In  producing  Bragg  scattering  from  elastic  and 
spin-elastic  microwave  packets  within  YIG  crystals.  This  high  resolu¬ 
tion  system  which  Is  now  operational  forms  the  basis  for  studies  of  spin- 
elastic  conversion  In  spatial  and/or  temporal  gradients  of  magnetic  field. 

The  doctoral  thesis  of  Benedict  Hu  has  been  completed  and  forms  the 
basis  of  this  report. 

Dr.  A.  Linz  has  been  assigned  responsibility  for  overseeing  the  top- 
seeded  solution  crystal  growth  and  chemical  vapor  deposition  program, 
Professor  D.J.  Epstein  for  crystal  evaluation  with  respect  to  magnetic 
properties  and  Professor  F.R.  Morgenthaler  for  microwave  spin  wave 
propagation  studies;  as  principal  Investigator  the  latter  has  overall 
responsibility  for  coordinating  the  various  phases  of  the  research. 


STUDIES  OF  MAGNETOSTATIC  WAVES  AND  MAGNETOELASTIC  WAVES 


IN  YIG  USING  OPTICAL  PROBING  AND  MICROWAVE  TECHNIQUES 


by 

Hung-Llang  Hu 


The  material  which  follows  is  Identical  in  content  to  a  thesis  submitted 
to  the  Department  of  Electrical  Engineering,  Massachusetts  Institute  of 
Technology,  on  January  29,  1971,  in  partial  fulfillment  of  the  require¬ 
ments  for  the  degree  of  Doctor  of  Philosophy  in  Electrical  Engineering. 
This  thesis  is  also  available  as  Microwave  and  Quantum  Magnetics  Group 
Technical  Report  No.  25,  February,  1971. 
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LIST  OF  SYMBOLS 

(other  ad  hoc  notation  is  explained  in  the  text  where  it  occurs) 


Symbol 

Definition 

Place  of  symbol  or 
first  occurrence 

a 

- 

half  length  of  major  or  minor  axis  of 

Eqs.  (2.16)  and 

the  elliptical  spin  wave  polarization 

(2.  17) 

a 

height  of  the  light  beam  cross  section 

Fig.  4.  2  and 

Eq.  (4.  14) 

k?  2 

A 

OJg  +  (Nii+^)  +  >^e*^ 

Eq.  (2.  lib) 

b 

— 

half  length  cf  major  or  minor  axis  of 

Eqs.  (2.16)  and 

the  elliptical  spin  wave  polarization 

(2.17) 

^ijk 

= 

first  order  (in  strain)  magnetoelastic 
constants 

Eq.  (A2.  1) 

b. 

ijk 

S 

small  signal  first  order  magneto¬ 

Eqs.  (1 .  3d)  and 

elastic  constants 

(A2.  8). 

h'ba 

= 

non-zero  first  order  magnetoelastic 
constants  of  cubic  crystals 

Eq.  (A2.6) 

z- directed  longitudinal  and  shear 

Eq.  (2.  37) 

Jt  o 

magnetoelastic  frequency  splittings, 
respectively 

b/i*  b 

V  s 

= 

z- directed  longitudinal  and  shear 
magnetoelastic  couplings,  respectively 

Eq.  (2.28) 

B.. 

ij 

= 

dielectric  impermeability  tensor 

Eq.  (4.  22) 

6B.. 

— 

small  change  of  dielectric  imperme¬ 

Eq.  (4.23) 

ability 

■  ^2  2 

B 

= 

^3  ^^22  ^  ^e*^ 

Eq.  (2.  11b) 

""ijkjC 

= 

second  order  elastic  constants 

Eq.  (1.4) 

^11'  ^12’ 

non-zero  second  order  elastic  con¬ 

Eq.  (1.4) 

= 

stants  of  cubic  crystals 

^44 

c 


velocity  of  light  in  free  space 


Eq.  (4.1) 


Symbol 


Definition 


Place  of  symbol  oi 
first  occurrence 


=  width  of  the  light  beam  cross  section 


=  a  proportionality  constant 


Fig.  4.  2  and 
Eq.  (4.14) 

Eq.  (4.  31b) 

Eq.  (2.  lid) 


=  mass  density  of  the  medium 


Eq.  (1.4) 


E  .  E 
o  o 


E(r') 

e‘v  Eli 


i-d  —d 

^H± 


/)»  S 


=  diameter  of  a  YIG  rod 

=  rf  electric  field  vector 

=  optical  electric  field  vector 

=  electric  field  of  the  incident  light 

=  electric  fields  of  the  incident  and 
diffracted  light,  respectively 

=  profile  of  the  light  beam 

=  vertically  polarized  (along  x-axis)  and 
horizontally  polarized  (along  z-axis) 
electric  fields  of  the  incident  light, 
respectively 

vertically  polarized  (along  x-axis)  and 
horizontally  polarized  (along  z-axis) 
electric  fields  of  the  diffracted  light, 
respectively;  the  subscript  (±)  denotes 
whether  the  diffracted  light  is  upshifted 
(+)  or  downshifted  (-)  in  frequency  with 
respect  to  the  incident  light 


magnetoelastic  force  density  (related  Eq.  (1.4) 


Eq.  (3.  la) 

Eq.  (1.1) 

Eq.  (4.1) 

Eq.  (4.  7a) 

Eq.  (4.  11 )  and 
Eq.  (4.30). 

Eq.  (4.11) 

Fig.  4.  2  and 
Eq.  (4.30) 


Fig.  4.  2  and 
Eq.  (4.  30). 


ijki’  ^mn 


to  h') 

=  small  signal  magnetoelastic  force 
density 

--  ^  frequency 

=  "linear"  magnetooptic  coefficients 


"quadratic"  magnetooptic  coefficients 


Eq.  (1.5) 


Eqs.  (4.40c) 
and  (A3.  9) 

Eqs.  (4.  40b) 
and  (A3. 13) 


Symbol 

Definition 

Place  of  symbol  or 
first  occurrence 

1?' 

small  signal  effective  magnetic  field 

Eq.  (1.2) 

ijme 

small  signal  effective  magnetic  aniso¬ 
tropy,  exchange  and  magnetoelastic 
fields,  respectively 

Eq.  (1.2) 

small  signal  Maxwellian  magnetic  fielc 

Eq.  (1.1a) 

-» 

H 

optical  magnetic  field  vector 

Eq.  (4.1) 

H.(r).  H.(z)  = 

internal  bias  magnetic  fields 

Eq.  (1.3a) 

H  ,  H 
o’  o 

external  bias  magnetic  fields 

Eq.  (1.  3a) 

„dem 

H  -  _ 

-V 

demagnetizing  field  along  the  static 
magnetization  direction 

Eq.  (1.  3a) 

„an  an 

effective  static  anisotropy  fields  along 
the  static  magnetization  direction 
(z-axis  or  x^-axis) 

Eqs.  (1 .  3a) 
and  (Al .  5c) 

(H®) 
c  c 

the  value  of  the  applied  field  for  which 
a  longitudinal  (shear)  magnetoelastic 
crossover  point  is  in  the  mid-point  of 
the  long  axis  of  an  axially  magnetized 
YIG  rod  or  bar 

Section  3.  3a 

H'(z) 

spatial  field  gradient  at  coordinate  z 

Eq.  (3.3) 

h'^  h'®  - 

crit’  crit  “ 

critical  spatial  field  gradients  for 
longitudinal  and  shear  magnetoelastic 
coupling,  respectively 

Eq.  (3.2) 

critical  temporal  field  gradient  for 
longitudinal  magnetoelastic  coupling 

Fig.  5.  13 

HT 

height  of  a  YIG  bar 

Eq.  (3.  lb) 

effective  current  density 

Eq.  (4.2b) 

t(k) 

wavevector  (wavenumber)  of  the 
elastic  waves,  spin  waves,  or  elastic/ 
spin  waves 

Eqs.  (2.7) 
and  (2 .  21 ) 

wavevector  (wavenumber)  of  the 
incident  light 

Eq.  (4.7) 

Symbol 


Definition 


Place  of  symbol  O] 
first  occurrence 


wavevector  (wavenumber)  of  the 
diffracted  light 

Eqs.  (4,12) 
and  (4.  13) 

II 

magnetic  anisotropy  constants  of  cubic 
crystals 

Eqs.  (Al.l) 
and  (2.15) 

L 

length  of  a  YIG  rod  or  bar 

Eq.  (3.1a) 

L  (S  ) 

cros  cros 

longitudinal  (shear)  magnetoelastic 
crossover  point 

Fig.  3.2 

m(r.  t) 

small  signal  magnetization  vector 
(in,  M  ^  0) 

Eq.  (1.1c) 

m  = 

small  signal  magnetization  amplitude 

Eq.  (4.41) 

M 

saturation  magnetization 

Eqs.  (1.2)  and 
(4.1) 

n  = 

index  of  refraction  in  an  isotropic 
medium 

Section  4. 1  and 

Eq.  (4.20) 

"i*  "d 

index  of  refraction  pertaining  to  the 
incident  and  diffracted  light 
polarization,  respectively 

Eqs.  (4.  6)  and 
(4.7) 

Nf.^ 

ij 

small  signal  anisotropy  tensor 

Eqs.  (1.3b)  and 
(A1.6) 

IJ 

small  signal  effective  anisotropy  tensor 

Eq.  (A1.9) 

Pa 

acoustic  power  density 

Eq.  (4.36) 

Pa 

power  of  a  circular  acoustic  beam 

Eq.  (4.37) 

-♦ 

P 

polarization 

Eq.  (4.1) 

-*eff 

effective  polarization 

Eq.  (4.3b) 

P  P  = 

ijkjf’  mn 

photoelastic  constants 

Eqs.  (4.23) 
and  (4.  24) 

Peff 

effective  photoelastic  constants 

Eq.  (4.36) 

r 


position  vector  from  an  origin  to  the 
observation  point 


Eq.  (4,  6)  and 
Fig.  4. 1 
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Symbol 


Definition 


Place  of  symbol  or| 
first  occurrence 


1  •  1  I. 

r  r ' 


il-»  -»i 

r-r 


position  vector  from  an  origin  to  the 
illuminated  volume 


unit  vectors 


Eq.  (4.  6)  and 
Fig.  4.  1 


Fig.  4.1 


total 

small  signal  total  Poynting  vector 

Eq.  (1.10) 

:*elas  ■^m 

s  » s  ,  _ 

-»me 

S 

small  signal  Poynting  vectors  of 
elastic,  magnetic  and  magnetoelastic 
waves,  respectively 

Eqs.  (1.11) 
and  (1 .  13) 

S,  S  = 

ki’  n 

symmetrical  strain  tensor 

Eqs.  (4.  23), 

(4.  24)  and  (A2.  1) 

S  (L  ) 

cros  cros 

shear  (longitudinal)  magnetoelastic 
crossover  point 

Fig.  3.2 

t 

time 

j 

Eqs.  (1.1)  and 
(4.  1) 

t' 

retarded  time 

Eq.  (4.6) 

transformation  matrix  and  its 
transpose 

Eqs.  (2. 4)  and 
(2.2) 

T  = 

Mjki 

’^EL'^  ^  "^ES^  ^  ’^SS 

Eq.  (3.4) 

T  T  = 

^EL'  ES 

time  taken  by  a  longitudinal  and  shear 
elastic  wave  packet  to  travel  the 
iength  of  a  YIG  rod 

Eq.  (3.4) 

T  'r  = 

^SL'  -^SS 

round  trip  time  of  flight  of  a  longi¬ 
tudinal  and  shear  elastic/spin  wave 
from  and  back  to  the  turning  point 

Eq.  (3.4) 

T 

turning  point 

Fig.  3.2 

V  = 

g 

group  velocity 

Eq.  (3.  5) 

group  velocity  of  spin  waves 

Fig.  2.  5 

V  = 

p 

group  velocity  of  (longitudinal  or  shear] 
elastic  waves 

Eq.  (3.  3) 

V/,,  V  = 

r  s 

group  velocities  of  longitudinal  and 
shear  elastic  waves,  respectively 

Eq.  (2.  26)  and 
Fig.  2.  5 

T 
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vv ' ,  w  = 

width  and  effective  width  of  the 
acoustic  or  spin  wave  in  the  plane  of 
incidence,  respectively 

Eqs.  (4.  14)  and 
(4.21) 

'WD 

width  of  a  YIG  bar 

Eq.  (3.  lb) 

w  = 

total 

small  signal  total  energy  density 

Eq.  (1.12) 

w  ,  w  , 
an  m 

w  ,  w  , 

me  elas 

small  signal  anisotropy,  magnetic, 
magnetoelastic  and  elastic  energy 
density,  respectively 

Eqs.  (1.3b), 

(1 . 6) ,  (1.  8)  and 

(1.7) 

W  ,  W 
an  me 

anisotropy  and  magnetoelastic  energy 
density,  respectively 

Eqs.  (A1 .  1 )  and 
(A2.1) 

X,  Xg), 

(x.  y,  z) 

a  coordinate  system  transformed  from 
the  "cubic-edge"  coordinate  system 
through  the  Euler  Angles  ((p,6.lp) 

Eq.  (2.1)  and 

Fig.  2.1 

z  .  = 

cr 

„s 

cr’  cr 

z*  = 

tp 

coordinate  of  either  longitudinal  or 
shear  magnetoelastic  crossover  point; 
coordinates  of  longitudinal  and  shear 
magneotelastic  crossover  points, 
respectively; 

coordinate  of  the  turning  point 

Eq.  (3.  6)  and 

Fig.  3.2 

Eq.  (3.3) 

m). 

M 

»y(=  tJ).  = 

M 

/  z  \ 

the  direction  ^osines  of  the  magnetiz¬ 
ation 

Eqs.  (4.40b), 

(4.  40c)  and  U.  41 ) 

a 

spin  wave  loss  in  db/^sec 

Eqs.  (3.7)  and 
(3.11) 

a 

the  angle  between  the  diffracted  light 
beam  and  the  wavefronts  of  elastic  or 
spin  waves 

Eq.  (4.  1  5)  and 

Fig.  4.2 

|3 

the  angle  between  the  diffracted  light 
beam  and  the  plane  of  incidence 
(y-z  plane) 

Eq.  (4.15) 

y 

gyromagnetic  ratio 

Eq.  (1.2) 
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-  eff 

6Cy 

= 

small  change  of  the  effective  dielectric 
tensor 

Eq.  (4.12b) 

-  eff 

6Xij 

lif 

= 

small  change  of  the  effective  dielectric 
susceptibility  tensor 

Eq.  (4.  8) 

6Ci^ 

= 

an  abbreviated  6c 
'it 

Eq.  (4.21) 

Ag 

= 

^11  ■  ^12  ■  ^^44 

Eq.  (4.46b) 

e 

= 

permittivity 

Eq.  (1.1) 

c 

= 

dielectric  constant 

Eq.  (4.1) 

Heff 

e 

= 

effective  dielectric  tensor 

Eq.  (4.  3c) 

c 

= 

an  abbreviated  c^^^ 

Eq.  (4.22b) 

V 

= 

spatial  or  temporal  magnon-phonon 
conversion  efficiency 

Eq.  (3.6) 

B. 

n 

= 

spatial  magnon-longitudinal  phonon 
conversion  efficiency 

Eq.  (3.7) 

— 

the  efficiency  of  conversion  from  shear 
elastic  waves  to  longitudinal  elastic 
waves  via  spin  wave  coupling 

Eq.  (3.  11) 

®B 

~ 

incidence  angle  or  Bragg  angle 

Eqs.  (4.  1  5), 

(4.  20)  and 

Fig.  4. 2 

= 

exchange  constant 

Eq.  (1.  3c) 

\ 

= 

wavelength  of  light  in  free  space 

Section  4. 1 

A 

= 

wavelength  of  spin  or  elastic  waves 

Section  4.  1  and 
Eq.  (4.19) 

= 

free  space  permeability 

Eq.  (1.1) 

M 


permeability  of  the  medium 


Section  4.  2 
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Symbol 
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Place  of  symbol  or 
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the  coordinate  system  in  which  the 
axes  are  along  the  cubic  edges  of  cubic 
crystals 

Eq.  (2.1)  and 

Fig.  2.1 

-4 

P 

small  signal  elastic  displacement 
vector 

Eq.  (1.4) 

((P.d.lp) 

the  Euler  Angles 

Eq.  (2.  2)  and 

Fig.  2.1 

«eff 

X 

effective  dielectric  susceptibility 
tensor 

Eq.  (4.  3a) 

CO  =  2irf  = 

microwave  angular  frequency 

Eq.  (2.21) 

‘"k 

spin  wave  angular  frequency 

Eq.  (2.7) 

Eq.  (1.2) 

i 

-yp  H. 

Eq.  (2.10a) 

2  2  2 

A  cos  i/j  +  B  sin  ip  -  C  sin  2^ 

Eq.  (2.18) 

2  2  2 

A  cos  -  B  sin  +  C  sin  2ip 

Eq.  (2.18) 

angular  frequencies  of  longitudinal  and 
shear  elastic  waves,  respectively. 

Eq.  (2.26) 

angular  frequencies  of  the  incident  and 
diffracted  light,  respectively 

Eqs.  (4.  7)  and 
(4.13) 

— 

-14- 


LIST  OF  FIGURES 


Fig.  No.  Page 

2.  1  The  Euler  Angles . 31 

2.2  Dispersion  diagram  for  spin  waves . 36 

2.  3  Longitudinal  and  shear  magnetoelastic  couplings  of 

z-directed  waves  in  a  {100}  piano . 43 

2.  4  Longitudinal  and  shear  magnetoelastic  couplings  of 

z-directed  waves  in  a  (1  lOj  plane . 45 


2.  5  Dispersion  diagram  for  z-directed  magnetoelastic 

waves;  the  spin/longitudinal  and  spin/shear  wove 
splittings  (bjf  and  bg)  are  shown  very  much 
exaggerated.  The  quantities  v^,  Vg  and  v^  are, 
respectively,  group  velocities  of  longitudinal  elastic, 
shear  elastic  and  spin  waves;  the  branch  numbers 
are  referred  to  in  the  text . 47 

3. 1  Block  diagram  of  L-band  microwave  circuitry.  50 

3.  2  Internal  field  profiles  of  a  YIG  rod  and  a  YIG  bar, 

showing  turning  points  (T)  and  crossover  points 

(Lcros  and  S^ros)*  distances  between  those 

points  are  very  much  exaggerated . 52 

3.  3  Measured  time  delay  of  the  first  longitudinal  (LOl) 
and  shear  (SOI)  echoes  versus  dc  magnetic  bias 
fields  at  1. 5  GHz . 55 

3.  4  Magnetoelastic  echoes  affected  by  a  pulsed  magnetic 
bias  field.  The  first  nulse  in  each  lower  trace  is 
due  to  leakage;  the  second  and  third  arc,  respec¬ 
tively,  the  longitudinal  (LOl)  and  shear  (SOI) 
magnetoelastic  echoes  (time  scale  1  ^scc/cm).  A 
full  explanation  of  cases  (a)  -  (c)  is  given  in  the  text.  57 

3.5  A  typical  magnetoelastic  echo  pattern  observed  with 
one-port  fine  wire  excitation.  The  first  pulse  in 
the  trace  is  due  to  leakage;  the  other  two  pulses 
are  due  to  the  longitudinal  and  shear  magnetoelastic 
coupling  (time  scale  1  ^scc/cm).  A  full  explanation 


is  given  in  the  text . 59 

3.  6  An  echo  pattern  of  Fig.  3.  5,  reconstructed  from 

theoretical  considerations . 61 


3.  7  Delay  of  elastically  injected  longitudinal  clastic/spin 
wave  packet  and  group  velocity  at  the  bar  center  as 
a  function  of  applied  magnetic  field  for  a  ”22.  5°”  YIG 
bar  which  measures  4  inm  x  4  mm  x  12  mm.  ...  64 


«.  No. 


I  f)- 


M  'l  lif  tlisl«)rt  Ion  of  a  Umfiituiiinal  olaHlic/spin  wavo 

pat  kt'l  wliifli  is  oxfitiMl  and  detci'lcd  by  a  ZnO 
transdiii'ci',  as  a  function  of  i?xtornal  applied 
inaj^iiflic  fii’ld  (time  scale  0.  I  psec/em).  ...  65 

;t. !)  Mafjnetoelastic  (ma^mon-phunon)  conver.sion 

efficiency  versus  inverse  spatial  magnetic  field 
gi-adient  ohtainetl  with  spin  wave  injection.  .  .  f58 

8.  10  Magnctoelastic  (niagnon>phonnn)  cor  'crsion 

efficiency  vei*sus  inverse  spatial  ma,.^  'tic  field 
gradient  obtained  with  acoustic  wave  injection.  71 

.1.  1 1  l*iiisc  compression  of  longitudinal  eln.stic/Rpin 

waves  subjected  to  pulsed  magnetic  fields.  The 

static  magnetic  field  here  is  high  enough  so  that 

tliei'c  ai'c  crossover  points  inside  the  YIG  bar 

(i.e.,  lit)  >  ^c)-  llorizontal  sweeps  arc  from 

left  to  riglit  witli  scale  0.  .‘ipsec/cm  (frequency 

1.5  011/.) .  74 

3.  12  l^ilsc  compression  of  longitudinal  elastic/spin 

waves  subjected  to  pulsed  ntagnetic  fields,  when 
llo  <  ile-  Horizontal  swetrps  are  from  left  to 
l  ight  with  scale  0.  5jisec/cm  (frequency  1. 5  GHz).  75 

3.  13  AiTangement  used  to  observe  spatial  conversion 

from  shear  clastic  waves  to  longitudinal  elastic 
waves  via  spin  wave  coupling . 7(1 

3. 14  'I'otal  delfiy  time  (tt)  and  spatial  conversion 

efficiency  (tjs+l,)  from  shear  elastic  waves  to 
longitudinal  elastic  waves  via  spin  wave  coupling 
as  a  function  of  field  gradient  at  the  longitudinal 
crossover  point . 7'J 

3.  1  5  Oscillograms  showing  the  spatial  conversion  from 

shear  elastic  waves  to  longitudinal  clastic  waves 

via  spin  wave  coupling.  The  horizontal  traces 

arc  from  left  to  right  (time  scale  1  ^sec/cm).  80 

4.  1  Radiation  from  an  illuminated  volume .  88 

4.2  Geometry  of  light  diffraction . 93 

4.3  Hragg  condition  for  + 1<. . 93 

4.4  Geometry  of  the  photoclastic  interaction.  .  .  .  100 

4.  5  Geometry  of  the  photoelastic  interaction  and 

orii'ntation  of  the  long  axis  (z-axis)  of  the  sample.  102 


-16- 


Fig.  No. 
4.6 


4.7 


Page 

The  geometry  of  light  diffraction  from  spin  waves; 

the  internal  field  profile  of  the  bar  is  also 

indicated .  107 

Light  diffraction  from  coherent  magnetostatic 

waves .  110 


5.  la  Block  diagram  of  the  optical  experimental  setup.  115 

5.1b  Optical  experimental  setup.  .  120 

5.  2  Diffracted  light  pulses  from  a  1 .  3  GHz  shear 

elastic  wave  pulse  propagating  along  a  (100> 
direction  and  polarized  along  a  (110)  direction  in 
YIG.  The  graph  was  recorded  with  the  signal 
from  the  PAR  boxcar  integrator.  The  insert  is 
the  scope  trace .  124 

5.  3  Diffracted  light  intensity  (in  arbitrary  units)  of 

shear  elastic/spin  waves  as  a  function  of 
external  magnetic  field;  the  light  beam  is  incident 
near  the  mid-point  of  the  long  axis.  The  insert 
is  the  scope  trace  swept  from  right  to  left  when 
Ho  -  H§,  indicating  small  group  velocity 
dispersion.  A  full  explanation  is  given  in  the 
text .  126 

5.  4  Diffracted  light  intensity  (in  arbitrary  units)  of 

shear  elastic/spin  waves  as  a  function  of  e.xternal 
magnetic  field;  the  light  beam  is  incident  near 
the  mid-point  of  the  long  axis.  The  insert  is  the 
scope  trace  swept  from  right  to  left  when  Hq  - 
indicating  the  peak  diff.'acted  light  intensity 
(marked  "spin"  on  the  curve)  and  a  slight  decrease 


of  the  group  velocity.  .  128 

5.  5  Diffracted  light  intensity  of  shear  elastic/spin 

waves  as  a  function  of  external  magnetic  fields, 
illustrating  the  effect  of  homogeneity  of  the 
applied  field .  130 

5.  6  Diffracted  light  intensity  of  shear  elastic /spin 

waves  as  a  function  of  external  magnetic  fields, 
illustrating  the  effect  of  homogeneity  of  the 
applied  field .  131 


5.  7  Diffracted  light  pulses  from  a  1 . 5  GHz  longi¬ 

tudinal  elastic  wave  pulse  propagating  along  a 
"22.  5°"  YIG  bar.  The  graph  was  recorded 
with  the  signals  from  the  PAR  boxcar 
integrator.  The  insert  is  the  scope  trace. 


135 


Fig.  No. 


-17- 


Page 


5.  8 


5.  9 


5.  10 

5.  11 


5.  12 


5.  13 


Scope  traces  indicated  with  signals  from  the  micro- 
wave  circuitry  and  the  photomultiplier  as  a  function 
of  external  magnetic  field.  The  upper  and  lower 
traces  are  synchronized;  the  upper  traces  are 
from  the  microwave  receiver  and  the  lower  traces 
are  from  the  photomultiplier.  The  diffracted  lig'.t 
signal  is  polarized  atff/2  with  respect  to  the 
incident  light  (E^,  E^_).  In  (a),  the  upper  trace 
indicates  the  first  leakage  pulse,  and  two  longi¬ 
tudinal  elastic  echoes  whose  amplitudes  saturate 
the  receiver.  In  (b),  the  longitudinal  echoes  are 
strongly  attenuated  due  to  the  magnetoelastic 
coupling.  In  (c),  the  longitudinal  echoes  are  still 
attenuated;  the  lower  trace  now  indicates  a  strong 
diffracted  light  signal  from  coherent  spin  waves 
(compare  the  diffracted  light  pulse  with  that  in 
Fig.  5.7) .  138 

Scattered  light  intensity  (in  arbitrary  units)  of 
spin  v.'aves  and  longitudinal  elastic  waves  as  a 
function  of  external  magnetic  field.  The  insert  is 
the  oscilloscope  trace  showing  the  peak  intensity 
of  the  diffracted  light  pulse  (compare  with  the 
diffracted  light  pulse  of  the  scope  trace  in  Fig.  5.  7).  140 

Diffracted  light  intensity  as  a  function  of  input 
microwave  power.  .  141 

Diffracted  light  intensity  (E^^,  E^+)  from  coherent 
spin  waves  as  a  function  of  the  timing  of  the 
pulsed  bias  field;  the  spin  waves  are  generated 
through  temporal  conversion  from  longitudinal 
elastic  waves  (time  scale  0.  5(^sec/cm).  A  full 
explanation  is  given  in  the  text.  .  144 

(a)  diffracted  light  pulses  from  longitudinal 
elastic  waves  (E^,  E^)  (time  scale  0.  5  pisec/cm). 

(b) .  to  (f)  -  an  expanded  version  of  Fig.  5.  11 

(E,\,E^,)  (time  scale  0.  2  tpec/cm).  ....  146 

V  li-\-  ^ 

Magnetoelastic  (phonon -♦  magnon)  conversion 

efficiency  versus  inverse  temporal  field  gradient. 

H-f  .  . . ■  .  147 

crit 


-18- 

INTRODUCTION 

In  the  recent  past,  the  propagation  of  magnetoelastic  waves  in  low 

loss,  saturated  ferrimagnetic  crystals  has  been  a  subject  of  great 

interest  because  of  their  potential  us>es  in  compact  microwave  signal 

processing  devices.  For  example,  frequency  filtering,  variable  delay, 

gating,  pulse  compression  and  expansion  are  all  possible  operations. 

However,  their  practical  development  requires  additional  information 

about  the  spin  wave  channel  of  power  flow,  together  with  its  coupling  to 

83 

the  well  understood  elastic  wave  channel.  In  1961,  Schlbmann 

predicted  that  spin  waves  could  be  excited  in  non-uniform  magnetic  bias 

84 

fields.  Shortly  afterwards,  Eshbach  reported,  on  experiments  with 

a  normally  magnetized  yttrium  iron  garnet  (YIG)  disk  in  which  the 

magnetic  field  is  non-uniform,  the  c;xcitation  of  spin  waves  by  rf  magnetic 

fields.  He  further  indicated  that  the  excited  spin  waves  travel  toward 

lower  magnetic  fields  with  increasing  wavenumbers  and  eventually 

convert  into  shear  elastic  waves.  Since  then,  extensive  studies  in  YIG 

have  been  carried  out  on  the  shear  elastic/spin  waves  propagating  in 

both  the  spatially  varying  magnetic  fields^’  and  temporally  varying 

8  18  20 

magnetic  fields.  '  ’  It  was  later  found  that  when  the  applied  field  is 

not  along  the  propagation  direction  of  the  shear  elastic/spin  waves,  the 
losses  increase  tremendously.  This  led  to  work  on  the  stability  of  spin 
wave  trajectories.  These  studies  were,  however,  handicapped  by  the 
lack  of  a  powerful  experimental  tool.  Light  probing  appears  to  be  a 
convenient  and  potentially  powerful  method,  since  the  technique  has 
already  been  widely  developed  to  in/estigate  elastic  waves.  However, 
in  the  case  of  YIG,  which  has  a  transparent  window  between  1 .  1  and  6 
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niicrons,  the  technique  is  of  much  reduced  power  due  to  the  compar¬ 
atively  weak  light  sources  and  low  sensitivity  detectors  available  in  this 

wavelength  range.  Besides,  in  contrast  to  the  exhaustive  work  already 

56  74-79 

done  on  elastic  waves,  very  little  '  has  been  done  with  optical 

spin  wave  interactions.  Accordingly,  one  of  the  prime  goals  of  this 

thesis  is  to  provide  improved  understanding  of  such  interactions. 

A  study  of  the  coupling  of  the  longitudinal  elastic  waves  and  spin 

waves  is  also  a  major  topic,  since  in  practice  it  is  easier  to  fabricate 

an  efficient  longitudinal  piezoelectric  transducer  than  a  shear  one.  The 

lack  of  prior  experimental  investigations  in  this  case  is  attributable  to 

the  fact  that  in  a  cubic  crystal  the  coupling  vanishes  identically  whenever 

the  wave  propagation  is  simultaneously  parallel  to  the  magnetization 

vector  and  any  one  of  the  three  commonly  employed  principal 

<100,  110,  111)  crystallographic  directions,  along  which  the  elastic  waves 

propagate  as  pure  modes.  Transverse^*^  or  obliquely  magnetized^ ^ 

YIG  rods  of  appropriate  orientations  do  provide  the  longitudinal  coupling 

1 3 

but  the  defocussing  effects  limit  this  approach. 

1  5 

Morgenthaler,  in  view  of  the  fact  that  YIG  is  nearly  elastically 
isotropic,  first  proposed  that  the  longitudinal  magnetoelastic  coupling 
for  waves  propagating  along  the  dc  magnetization  vector  may  be  obtained 
if  th  propagation  direction  is  not  chosen  along  any  one  of  the  principal 
directions  (100,  110,  111)  but  rather  along  a  direction  chosen  to  optimize 
the  longitudinal  magnetoelastic  coupling.  He  found  that  the  direction 

is  in  a  {110}  plane  at  an  angle  of  25.  52°  from  a  (100)  axis.  Subsequently 

12  o 

in  our  experiments,  using  a  rod  oriented  in  a  {100}  plane  at  (22.  5°) 

from  the  (100)  axis  which  provides  about  90%  of  the  maximum  frequency 

splitting,  we  have  indeed  observed  the  spin/longitudinal  elastic  wave 
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propagation.  Similar  observations  were  also  reported  by  Lewis,  who 
used  a  YIG  rod  oriented  in  a  {100}  plane  at  20°  from  the  <100)  axis.  This 
opened  up  the  possibility  of  studying  not  only  the  various  phenomena  of 
longitudinal  elastic/spin  waves  in  spatially  and/or  temporally  varying 
fields,  similar  to  the  shear  wave  case,  but  also  the  conversion  of  longi¬ 
tudinal  elastic  waves  to  shear  waves,  and  vice  versa,  via  spin  wave 
coupling.  This  is  the  second  topic  to  which  this  thesis  is  devoted. 

In  Chapter  I,  the  equations  of  motion  for  the  coupled  magnetic  and 
elastic  system  are  reviewed,  in  order  to  provide  all  the  necessary 
equations  needed  in  the  following  chapter.  In  Chapter  II,  Morgenthaler' s 
results  are  rederived  and  summarized.  In  Chapter  III,  extensive 
experimental  studies  of  longitudinal  and  shear  elastic /spin  waves,  in 
spatially  and/or  tempor  ally  varying  magnetic  fields,  are  described. 
Topics  include  variable  delay  of  magnetically  injected  spin/longitudinal 
elastic  waves,  and  conversion  from  the  longitudinal  elastic  to  spin  waves 
in  spatially  varying  fields.  Frequency  shifts,  and  pulse  compression  of 
the  longitudinal  elastic /spin  waves  subjected  to  pulsed  magnetic  fields 
have  also  been  obtained.  Some  preliminary  results  on  shear /longitudinal 
elastic  wave  spatial  conversion  via  spin  wave  coupling  are  presented 
and  discussed. 

In  Chapter  IV,  a  review  of  light  diffraction  is  given.  It  is  then 

7  8 

followed  by  a  new  theory  describing  the  phenomena  of  coherent  light 

diffraction  from  coherent  spin  waves.  The  theory,  in  contrast  to  prior 

56 

work  of  Auld  and  Wilson,  has  appropriately  taken  into  account  all  the 
terms  linear  in  the  small  signal  spin  amplitudes  for  the  perturbed 
dielectric  tensor;  the  latter  is  the  one  that  gives  rise  to  the  light 
scattering  from  spin  or  elastic  waves.  As  a  result,  the  theory  predicts 
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that  the  intensity  of  the  light  diffracted  from  spin  waves,  nnlike  that  from 
elastic  waves,  depends  on  the  incident  light  polarization  and  whether  the 
diffracted  light  is  upshifted  or  downshifted  in  frequency.  This  appears 
to  explain  satisfactorily  not  only  our  experiments  described  in  Chapter  V 
but  also  observations  that  previous  authors^^'  reported  as  "anomalous. 

In  Chapter  V,  experimental  work  on  light  diffraction  from  elastic 
and/or  spin  waves  is  reported.  Smith's  work^^'  on  light  diffracted 
from  shear  elastic /spin  waves  was  repeated  here,  but,  unlike  his  work 
in  which  he  used  large  microwave  power  input  and  therefore  non-linear 
interaction  probably  occurred,  the  microwave  power  we  used  was  much 
reduced  and  within  the  linear  interaction.  From  the  experiments,  we 
showed  that  the  anomaly  he  reported  indeed  comes  from  the  interaction 
of  light  and  spin  waves. 

In  section  5.  3,  we  report  the  first  direct  observation  of  strong 

Bragg  scattered  infrared  light  (1150  nm)  from  microwave  coherent  spin 
7  8 

waves.  The  experiments  utilized  coherent  spin  waves  generated 
through  efficient  space-gradient  conversion  from  longitudinal  elastic 
waves;  the  spin  waves  have  negligible  change  of  wavenumber,  are 
converted  more  rapidly  (therefore  with  less  loss)  and  the  scattered  light 
can  easily  be  distinguished  from  that  due  to  longitudinal  elastic  waves. 

The  scattered  light  intensity  is  found  to  be  dependent  on  the  incident 
light  polarization  and  the  sign  of  frequency  shift,  again  in  accord  with 
our  theory. 

Finally,  light  scattering  from  spin  waves  has  been  used  to  measure 
the  temporal  conversion  efficiency  of  longitudinal  elastic  waves  to  spin 


waves. 


CHAPTER  I 

SMALL  SIGNAL  EQUATIONS  OF  MOTION  OF  A 
MAGNETOELASTIC  FERROMAGNET 

In  this  chapter,  small  signal  equations  of  motion  for  the  electro¬ 
magnetic,  spin,  and  elastic  system  in  a  lossless  magnetoelastic 
ferromagnet  are  reviewed  briefly  so  as  to  provide  all  the  necessary 
equations  used  in  the  later  chapters.  Conservation  theorems  for  small 

signal  power  flow  are  also  discussed.  For  more  detailed,  in-depth 

1  _  6 

discussion  one  may  refer  to  many  excellent  articles, 

4-6 

1 . 1  Small  Signal  Equations  of  Motion 

In  a  lossless  saturated  ferromagnet  insulator,  small  signal  equations 
of  motion  governing  the  magnetization  vector  and  elastic  displacement  of 
the  medium  are  described  in  terms  of  the  interaction  among  electro- 


magnetic,  spin,  and  elastic  systems. 

For  the  electromagnetic  system  Maxwell's  equations 

are: 

7  X  li  -  c  ^  =  0 

(1.1a) 

7xe+U5T-(h+m)  =  0 

”0  0 1 

(1.1b) 

7  •  (li  +  in)  =  0 

(1.1c) 

7  •  (ce)  =  0 

(1.  Id) 

where  the  conductivity  has  been  assumed  to  be  zero,  the  permittivity 

C  is  assumed  to  be  a  strain  independent  scalar,  and  the 

-♦  5  6 

dilational  variation  of  the  small  signal  magnetization  m  is  neglected.  ’ 
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'rho  vectors  h,  e  are  the  smal)  signal  Maxwellian  magnetic  and  electric 


field  vectors  and  is  the  free  spac(?  permeability. 

The  linearizc'd  torque  equation  for  the  small  signal  magnetization  m 


is  given  by; 


?t  "  \  M  m  +  h  +  hj  (1.2) 


where  M  =  is  the  static  magnetization  vector,  m(r,  t)  the  small 

deviation  from  equilibrium  (and  m*  M  =0),  w^^  =  and  y 

gyromagnetic  ratio  (negative  for  electron),  H^(r)  the  total  static  magnetic 
field  (applied  H^,  anisotropy  ,  plus  demagnetizing  -N^^M)  inside  the 
medium  and  ii  Maxwellian  rf  field;  the  vector  ii'  is  an  effective  field 
which  gives  rise  to  torque  density  from  magnetic  anisotropy  n  , 

^0X  ^1X10 

exchange  h  ,  magnetostriction  h  '  and  anything  else.  The  effective 
magnetic  field  components  have  been  discussed  in  great  detail  by 


Morgenthaler. 


In  the  simplest  case,  where  isotropic  exchange 


interaction  is  assumed  and  high  order  interactions  neglected,  they  can 
be  written  as  (with  the  static  magnetization  vector  in  the  +z  direction); 


II.  = 
1 


H  -  N,.M  +  H':. 
o  M  M 


(1.3a) 


Mo  »«^i 
1] 


I  (-N^>.)  (i  =  l,2) 
j  =  l  ^  ^ 


r"  (i  j) 


(1. 3b) 


, ex  .  _2 

h.  =  \  V  m. 
1  e  1 


(i  =  1,2) 


(1.  3c) 


,  aw  1  3 

±  _JEe  =  .  _L  j  b. 


u  am. 
'^o  1 


o  ,  k=l 


ijk  ax. 


^jk  "  ^kj 


(1.3d) 


The  quantity  Nr.  is  the  small  signal  anisotropy  tensor,  \  the  isotropic 


/ 
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exchange  constant,  the  small  signal  magnetoelastic  energy  density, 

Wan  the  small  signal  anisotropy  enei'gy  density,  p  the  displacement 
vector,  and  the  first  order  (in  strain)  magnetoelastic  tensor.  The 
small  signal  anisotropy  tensor  and  the  first  order  magnetoelastic  tensor 
are  discussed  in  detail,  respectively,  in  Appendix  I  and  Appendix  II. 

For  the  elastic  subsystem,  use  of  the  linearized  Newton  force  law 
gives: 


4 


3 

£ 

j,  k,i=l 


h  /  ^PkV  -t 
ax.  \  ijki  ax^/  1 


(i  =  l,2,  3)  (1.4) 


where  is  the  mass  density  of  the  medium,  second  order  elastic 

constant  tensor,  and?'  the  magnetoelastic  force  density  (related  to  ii') 
that  includes  first  order  electromagnetic  and  interaction  forces.  In  the 
cubic  ferrimagnetic  YIG,  terms  such  as  the  magnetic  body  force  are 
negligible;  therefore  f  is  dominated  by  the  linear  magnetoelastic 
coupling  and  can  be  well  approximated  as: 

3  ^  ,  aw _ V  3  am. 


I  I  j.k.l 


(1.5) 


Equations  (1.1),  (1.  2),  and  (1.4)  are  a  set  of  time-dependent 

coupled  differential  equations,  and  in  general  are  very  complicated  to 

solve.  Since  in  this  thesis  we  are  primarily  concerned  with  the 

magnetoelastic  region  of  the  spectrum,  values  of  the  wavenumber  k  are 
4  -1 

on  the  order  of  10  cm  .  In  this  range  electromagnetic  effects  are 
small  and  the  electric  field  e  in  Eq.  (1.1a)  can  be  neglected  to  the  first 
approximation.  This  results  in  the  so-called  "magnetostatic 

7 

approximation",  which  we  will  use  throughout  the  thesis. 


5  6  8 

1 . 2  Small  Signal  Power  Energy  Conservation  Theorem^*  ' 

The  contributions  to  the  small  signal  energy  djnsity  of  the  system 
described  by  Eq.  (1.2)  arise  from  magnetic,  elastic,  and  magnetoelastic 
terms.  The  magnetic  energy  results  from  Zeeman,  anisotropy,  exchange, 
and  dipolar  contributions. 


w 


m 


^1^0 


2 

I 

i.  j=l 


Tvran 

N..  m.m. 
iJ  1  1 


(1.6) 


The  elastic  term  contains  the  kinetic  and  potential  energies 


''^elas  ~  2  ^o 


at 


(1.7) 


The  magnetoelastic  component  is 


w 


me 


2  3  ap. 

=  L  Z  b. .  m.  ^ 

i=lj,k=l 


(1.8) 


Notice  that  the  energy  associated  with  the  small  signal  electric  field 
e  is  of  higher  order  and  is  not  included  here.  In  the  Poynting  vector, 
however,  the  electromagnetic  term  is  of  first  order  and  must  be  taken 
into  account.  Therefore  Maxwell's  equations  (1.1)  give  rise  to  the 
small  signal  electromagnetic  Poynting  theorem 


and  (1. 4)  to  the  small  signal  elastic  Poynting  theorem 


where 
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(1.11) 


^elas 

i 


r 

ijWSx^M 


In  addition  to  the  above  two  Poynting  theorems  governing  the  electro¬ 
magnetic  and  elastic  energy  subsystems,  one  may  also  introduce  other 

12  5 

Poynting  theorems  *  *  for  a  material  and  a  magnetoelastic  sub.system. 
These  conservation  theorems  will  not  be  discussed  here  because  no  use 
of  them  will  be  made  later. 

Finally  the  Poynting  theorem  for  the  total  system  is  found  by 

summing  the  contributions  of  all  the  subsystems.  In  a  time- varying 

5  6 

but  spatially  uniform  magnetic  bias  field,  the  result  is;  ’ 


-♦total  aw  , 

,.s  .  „,.2, 


where 


•♦total  ^  ^  -♦elas  -*m  •♦me 

S  =(exh)  +  S  +S  +S 


(1.13) 


W,  ,  ,  =  W  +  w  ,  +  w 

totSi.  an  elas  me 


(1.14) 


The  other  small  signal  conservation  law,  namely  the  Stress- 

5  6 

Momentum  Conservation,  has  been  discussed  by  Morgenthaler^'  °  and 
will  not  be  repeated  here.  It  is  sufficient  to  summarize  the  conclusions 
that  one  may  draw  from  these  conservation  laws  .  In  a  time-varying 
but  spatially  uniform  magnetic  bias  field  the  propagation  of  either 
magnetostatic  or  magnetoelastic  waves  occurs  at  constant  wavenumber 
and  constant  average  small  signal  momentum  but  with  variable  frequency, 
average  power  and  average  energy.  On  the  other  hand,  in  a  time 
invariant  but  spatially  varying  magnetic  bias  field  the  propagation  of  the 
waves  occurs  at  constant  frequency,  average  power  and  average  energy 


7 


but  witli  variable  waveiiuniboi'. 
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and  c.verage  small  signal  momentum. 


6.  8 
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CHAPTER  II 

PROPAGATION  OF  Z- DIRECTED  MAGNETOELASTIC  WAVES 
IN  A  GENERAL  CRYSTALLOGRAPHIC  DIREC'ITON 

In  the  past,  numerous  studies  of  shear  elastic/spin  wave  interactions 

in  single  crystal  YIG  have  been  reported;^'  comparatively  few 

seem  to  have  been  reported  of  the  longitudinal  elastic/spin  wave  coupling, 

1  2 

and  none  until  the  recent  work  of  Morgenthalcr  et  al.  for  the  important 

case  of  wave  propagation  parallel  to  the  magnetization  vector  (z-directed 

magnetoelastic  waves).  This  dearth  of  experiments  is  attributable  to 

the  fact  that  whereas  in  a  cubic  crystal  the  shear  magnetoelastic  coupling 

is  present  regardless  of  the  crystal  orientation,  this  is  not  true  for  the 

longitudinal  magnetoelastic  coupling.  In  fact,  the  latter  vanishes 

identically  whenever  the  wave  propagation  is  simulatneously  parallel  to 

the  magnetization  vector  and  any  one  of  the  three  commonly  employed 

principal  crystallographic  directions  <100,  110,  111).  Transversely 

magnetized  rods  of  appropriate  orientations  can  provide  longitudinal 

magnetoelastic  coupling  but  do  not  appear  to  be  very  useful  because  of 

1  3 

strong  defocusing  effects.  Although  slight  tilting  of  a  (100)  rod  with 

respect  to  the  magnetic  field  (magnetization  vector)  does  allow  the 

observation  of  a  longitudinal  elastic/spin  wave  interaction,  defocusing 

1 3 

fields  limit  this  approach;  in  addition,  such  experiments  are  not 
easily  analyzed. 

1  5 

Morgenthaler  first  pointed  out  that  these  problems  may  be  over¬ 
come  if  the  waves  are  allowed  to  propagate  along  the  magnetization 
vector  (z-directed)  when  the  latter  is  parallel  to  some  non-principal 
crystallographic  direction  chosen  so  as  to  optimize  the  longitudinal 
magnetoelastic  coupling.  Although  in  this  case  the  acoustic  branch 
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('lastic  waves  are  in  general  neither  purely  longitudinal  nor  transverse 

1  6 

(shear),  tlu'  I'lastie  anisotropy  of  YIG  given  by 

2  c^, 

1  -  - =  0.053 

^11-^12 

is  small  enough  so  that  they  may  be  approximated  as  pure  modes. 

1 2 

Subsequently  we  reported  the  first  experimental  observations  of 
longitudinal  elastic/spin  wave  interaction.  These  detailed  (Experimental 
results  will  be  given  in  Chapter  III. 

In  this  chapter  we  will  rederive  the  dispersion  relations  of  magneto¬ 
elastic  waves  propagating  in  an  arbitrary  crystallographic  direction 
assuming  elastic  isotropy  but  magnetoelastic  anisotropy.  This  greatly 
simplifies  the  mathematical  calculation  and  also  provides  a  very  good 
approximation  insofar  as  longitudinal  elastic /spin  wave  and/or 
degenerate  shear  elastic/spin  wave  interactions  are  concerned. 

However,  in  the  case  where  two  non-degenerate  shear  elastic  modes  are 
involved,  the  calculation  of  shear  elastic/spin  wave  coupling  neglecting 
elastic  anisotropy  is  quite  poor  owing  to  the  comparable  magnitude  of 
elastic  and  magnetoelastic  frequency  splittings. 

We  begin  by  studying  the  frequency  spectrum  of  the  uncoupled  spin 
waves  propagating  in  an  arbitrary  crystallographic  direction.  After 
introducing  magnetoelastic  coupling  and  obtaining  the  dispersion 
relations  of  z-directed  magnetoelastic  waves,  we  look  for  the  particular 
crystallographic  direction  which  will  give  the  maximum  value  of  the 
square  of  the  longitudinal  magnetoelastic  splitting  of  frequencies 
(henceforth  referred  to  as  longitudinal  magnetoelastic  coupling).  This 
direction  turns  out  to  lie  in  a  {110}  plane,  at  an  angle  of  (25.  52®)  from 
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a  (100)  axis.  However,  the  experiments  to  be  described  in  the  next 
chapter  were  done  using  a  different  crystallographic  direction,  i.  e. ,  in 
a  {100}  plane  at  an  angle  of  (22.  5°)  from  a  (100)  axis.  In  this  direction 
the  longitudinal  magnetoelastic  coupling  is  approximately  80%  of  the 
maximum  value.  The  latter  direction  was  chosen  because,  for 

-4 

simplicity,  the  derivation  was  first  carried  out  for  the  magnetization  M 

constrained  to  lie  in  a  {100}  plane;  in  this  case  (22.  5°)  is  the  optimum 

angle.  Finally  it  is  shown  that  the  longitudinal  magnetoelastic  coupling 

2 

is  proportional  to  (b^  -  b2)  (b^  and  b2  being  the  ordinary  magnetoelastic 

constants),  and  therefore  vanishes  identically  in  the  case  of  magneto¬ 
elastic  isotropy,  i.  e. ,  =  b2. 

2 .  1  The  Euler  Angles^ 

Throughout  the  entire  thesis,  we  always  use  a  coordinate  system 

•4 

((x,  y,  z )  or  (x^ ,  X2,  x^))  in  which  the  magnetization  vector  M  is  along  the 
positive  z-axis  (or  x.^-axis  as  the  case  may  be).  This  coordinate 
system  is  then  related  to  the  usual  coordinate  system  (^,77,^)  iti  which  the 
cubic  edges  are  chosen  as  the  coordinate  axes  by  the  Euler  Angles 
((p,0,i/;).  The  latter  are  defined  in  Fig.  2.1  as: 

X  =  R^  ,(^/;)  R^,,(e)  R^(<p)  4 

where  x  =  (x,  y,  z)  or  (x^ ,  X2,  x^^) 

i  =  ii.v.O 


The  operator  R^((p)  denotes  a  rotation  through (p  about  the  j^-axis,  and 
thus  gives  rise  to  a  coordinate  system  (x",y",z");  the  operator  RyM(6) 
denotes  a  rotation  through  0  about  the  y"-axis  and  thereby  transforms  to 
a  coordinate  system  (x',y',z');  finally,  the  operator  R^,(i/;),  a  rotation 
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through  ip  about  the  z'-axis  takes  ont  to  the  coordinate  system  (x,  y,  z). 
Corresponding  to  these  three  operators,  one  defines  a  transformation 
matrix  as; 

X  =  (2.1) 

where 

f^<p,e,0)  =  R^,(ii))  Ry,(e)  R^((p) 


COS0  sin0  0 

COS0  0  -sin0 

coS(p  sincp  L 

-sin0  cosij)  0 

0  1  0 

-sin(p  c0S(p  0 

_  0  0  1  _ 

i_  sin0  0  COS0 

_  0  0  1 

(2.2) 

cos0cos0cos^- sin0sin(p  cos^!jcos6sin^+ sin0cos(p  -cosj/jsine 
-cosj/j  sin^- sini/;cos  6cos<p  +cosj()coS(p- sinj/jcosQsincp  sin^sinQ 
sinftcos^  sinesin(p  cos0 


Since  T  is  orthonormal,  the  inverse  of  T  is  equivalent  to  the  transpose 
of  i.  e.,  T. 

4  =  f(<p.  a.^)x  (2.3) 


where  for  ^  =  0; 


t(<p,e.o)  = 


COS0  COS^P 

cos  0  sin<p 
-sin0 


-  sin<p  sin0  cos(p 
costp  sin0  sinip 
0  COS0 


(2.4) 


and  in  general  for  ^  0 ; 


T(<p,0,j/))  = 

cos</«os0cos<p-  sinj/jsin<p 
cosj/)cos0sin(p+  sin</;cos<p 
-cosj/)sin0 


-sin0cos0cos<p-  cosi/)sin<p 
-sin^cos0sin(p+  cos0cos<p 
sinj/)sin0 


sin0cos<p 

sin0sin<p 

COS0 


(2.5) 
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2.2  Eigonvaiues  and  Eigenraodes  of  Spin  Waves  Propagating 
along  an  Arbitrary  Crystallographic  Direction 

In  this  section,  the  niagnetoelastic  interaction  is  neglected  so  that 
one  may  concentrate  on  the  spin  waves.  From  Eq.  (1.2),  one  may  write 
(neglecting 


ftm^ 

at 

/“i 

^  V  M  ^^2 

am„ 

/  “i 

at 

,  ,  an  .  ex> 

“  ^2  ■  ^2  ■  ^2  > 


(2.6a) 


(11.  6b) 


The  magnetization  M  is  along  the  positive  x^-axis  (or  z-axis)  and  the 
coordinate  system  (x^,  x^)  used  here  is  related  to  the  "cubic-edge" 

coordinate  system  ( ^, tj, (^ )  by  the  Euler  Angle  ((p,6,0)  (i.  e. ,  with^  =  0). 
The  quantities  h?*^  and  hf^  are  given  before  in  Eq.  (1. 3) 

(with  -♦  Ng  and  Maxwellian  field  can  be  found 

through  Eqs.  (1.1).  Since  we  are  concerned  here  only  v/ith  the  eigen¬ 
value  problem,  we  assume; 


m  (r,  t) 

ii(r,  t) 
e(r, t) 


Re  m  e 


j(w.t-k- 
Rel  he  ^ 

r- 

ReLee 


(2.7a) 

(2.7b) 

(2.7c) 


where  is  the  angular  frequency  and  k  =  (kj^,  k^,  k^)  the  wavevector. 

Substituting  Eqs.  (2.  7)  into  Eqs.  (1, 1)  and  (2.  6),  while  using  the 

"magnetostatic  approximation"  mentioned  earlier,  one  obtains; 

(m.k, +m„k„) 


(2.  8) 


-:u- 


} 


^  =  -  "m<N21  +  ^>'"1  -[“3  ^  <^22  +  ’^''"“m] 


(2.9a) 


[w3  +  (Nii +-^)C0M  +  Xk 


(2.  9b) 


where 


“3  =  -Wo«i  =  -Wo(Ho-N3M+Hf) 


(2. 10a) 


.an 


All 

In  Appendix  I, both  Ht"  and  N^.  are  derived  in  terms  of  the  Euler  Angles 
((p,6j  0)  and  the  anisotropy  constant  K^.  Furthermore,  it  is  shown  there 
that  one  may  write 


“3  '  '  'Wo  <«o  ■  '^3“> 


(2.10b) 


,an 


All  *  ctn 

provided  that  N?j“  is  replaced  by  defined  in  Appendix  T.  We  rewrite 
Eq.  (2.10)  into  a  matrix  form: 


d 

dt 


«  ^ 

mi 

nio 

2_ 

-C 

A 


-B 

C 


m. 


m. 


(2.11a) 


where 


A  =  "3*<Nii";i>‘>’M^^‘''“M 

k 


(2.  11b) 


B  =  Wg  +  (N22  +  ^2^  ‘*^M  '*'^6*^  ‘*^M 


(2. 11c) 


k.  k_ 


(2.  lid) 


In  order  to  find  the  major  and  minor  axes  of  the  elliptical  spin 
polarization,  it  is  helpful  to  introduce  the  third  Euler  Angle  tjj,  i.  e. ,  a 
rotation  through  ^  about  the  z-axis.  We  thus  write: 


5- 


cos^ 

-sin^^ 

m„ 

L 

sin^^ 

cos^ 

(2.12) 


where  and  are  along  the  major  and  minor  axes  of  the  ellipse. 
Substituting  Eq.  (2.12)  into  (2.  11a),  one  obtains: 


- 

— 

“ 

r 

“  — 

COSi/; 

+sini^ 

-C 

-B 

COSfp 

-sinj/j 

“2 

-  sin0 

cosilt 

A 

C 

sin^[) 

COS^jl) 

“2 

- 

— 

— 

.. 

L  J 

-  C  cos  2}li  -  '2  sin  2  0  (A  -  B) 
2  2 

A  sin  Bcos  +  C  sin  2^ 


2  2 

-(A cos  Bsin  C  sin2|/) 
C  cos  2^  +  1  sin2^(A  -  B) 


(2.13) 


We  now  demand  that  the  diagonal  elements  of  the  above  matrix  be  equal 
to  zero.  Thus: 


2P 

tan2(^  =  = 


^^^12+  ) 


l(r 


,2  ,2 

k„  k. 


(^22’Nll)'^(rl'4) 


k“  k‘ 


(2.14) 


where,  from  Appendix  I, 

2K 


Njj  =  - \  ^-|sin^2e^l  -  ;|sin^2(p^J 


u  M 

^o 


(2.  15a) 


^^22  =  — -4  [-fsin^esin^Zip] 


u  M 

^o 

2K 


1  r  3 

Nj2  “  - 2  |_"g  sinesin2  0sin4(p 


u  M 


(2. 15b) 


(2. 15c) 


From  Eqs.  (2. 14)  and  (2. 15),  one  realizes  that  the  angle  in  general  is  a 
complicated  function  of  Q.tp,  and  the  direction  of  propagation  It/  |lt| . 


Solving  Eq.  (2.  13)  subject  to  the  condition  of  Eq.  (2,  14),  one  obtains 


where 


and 


‘^1 

A  2 

=  A  cos  0 

+  I>  sin^0  -  C 

sin  2  0 

(2.  18a) 

‘^2 

=  A  sin  0 

2 

+  Bcos  0  +  C 

sin  2  0 

(2.18b) 

k 

II 

CM 

3 

3 

11 

AB  - 

(2.18c) 

A  plot  of  vs  k  (with  and  as  parameters)  is  shown  in 
Fig.  2.  2.  It  can  be  shown  that  the  lower  bound  of  the  spin  wave  manifold 
is  the  z-directed  waves  (i.  e. ,  kj  =  kg  =  0).  For  YIG,  ^22 

and  Nj2^]y[  usually  small  compared  to  at  microwave  frequencies. 

Therefore  one  may  deduce  from  Eq.  (2.  18c)  that  the  upper  bound  of  the 
spin  wave  manifold  is  approximately  the  waves  with  wavevector 
perpendicular  to  the  dc  magnetization  M. 


Fig.  2.  2  Dispersion  diagram  for  spin  waves. 
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2.  3  Dispersion  of  z- Directed  Magnetoelastic  Waves  in  an  Elastically 
Isotropic  but  Ala^netoelastic  Anisotropic  "Cubic"  Medium 

W('  assume  here'  that  the  coordinate  system  (xp  X2,  x^)  used  here  has 
been  iransformed  through  the  Euler  Angle  (<p,  from  the  "cubic-edge" 
coordinate  system  (4,77,  magnetization  is  along  the 

z-axis  and  the  major  and  minor  axes  of  the  elliptical  spin  wave  polariza¬ 
tion  are  along  the  x^  and  X2  axes.  Then  from  Eqs,  (1. 2)  and  (1.  3),  the 
coupled  spin  wave  equations  of  motion  are  written  as; 


at 


-60,  m„+  I  - b„.,  -r-^ 

j,  k=l  ^  k 


(2.19a) 


am2 

at 


0^2  m 


--Mb, 


1  .  ,  ,  u  Ijk  ax, 

J,  k=l  ^o  k 


(2.19b) 


where  oOj^  and  6O2  ^fe  given  by  Eq.  (2. 18)  and  b^^j^  are  given  in  Appendix  11. 
For  the  isotropic  elastic  case,  we  can  write  from  Eq.  (1.4) 


=  ^44^  ^i  <^12  ^44)^Va3^  ^ 


(2.20) 


+  E  b,  ..  r— ^ 
j.k=l 


(i  =  1,  2,  3) 


where  c^2  are  the  elastic  constants  and  c^^^^  2:  c^2  2c^^for  YIG 

As  in  Eq.  (2,  7)  we  assume 


m(r,  t)  =  Rej(m  e^^^  ^  (2.21a) 

p(r,  t)  =  Re|(£  ~  (2.21b) 


where  k  =  (k^,  k2,  k^). 
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Substituting  Eqs.  (2.21)  into  Eqs.  (2.19),  we  obtain; 


S  j  b, ,  .k. 
j=l  J  J 


I  ^  .  I  ^  . 

“JW  ^  la  ^llk^ki^  lu  ^12k^k 
k=l  Mo  f*o  ^ 


?,i^21k^j 

J=1  •’ 


FJWj 


/  J'’22k''kl  S  jV“‘>23k''k  El 


2  2  ' 

-D^w  H-c^^k  !  (Cj2+C44)k^k2 


+(Ci2-+C44)ki 


|+(Ci2+C44)k2 


^^12'^^44)*'i*^3|  ^^12'^^44^*^2'^3 


3  OOjy, 

f  ,1  JI“'^13k'^k  —2 
k=l  ^o 


^^12''’^44^*^1^2  fi-1 


*‘^12'^‘^44*''2'‘3  ^2 


■°o“^+‘'44‘‘^  e-3 

+(Cj2+C44)k^ 


(2.22) 


The  eigenvalues  can  be  found  by  equating  the  determinant  of  the  square 
matrix  in  Eq.  (2.  22)  to  zero.  Here  we  are  interested  only  in  z -directed 
waves;  therefore  k^  =  k2  =  0,  but  k^  =  k,  and 

.^M,  .^M  .**^M 

^2  1^^213*^  ijr"^223*^  1J“^233*^ 

^o  ^o  ^o 

.^M  , 

^1  ’1^  1^^113^  ijr"^123*^  l|r~^133^ 

O  "o  _  n 


jw 

^^2 

AVli  1 

1^  ^213^ 

i-^b  k 
1^^  ^223*^ 

-jco 

^O 

lu  ^123^ 

'^o 

1^13*^ 

1^213‘^ 

0 

1^123^ 

1  ^223^ 

0 

-D  ui^+c.  .k^ 
0  44 

1^133^ 

1^233^ 

0 

0 

-D  w^+c,  k^ 
o  11 


(2.23) 


\ 
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Hcrall  an  idinitity  ol  clolcM'nunants 


Cll.'t 


W  X 

Y  Z 


dot  (W  -  XZ'^Y)  det  Z  (2.24) 


whot'o  W  and  Z  in  the  diagonal  must  be  square  matrices.  Using  (2.24) 
one  obtains: 


.2  .  -  2  2.  ,  2  2.2 

(u)  -00^0:2)  (u)  -CCg)  (ui  -COg) 

.2 

,2  2,;,  2  2w  ,2  ,  ^2  ,  ,  ,  2  2, 

=  fy-— (w  -tUg)|(tu  -a)s)(Wibj33+uJ2b233)  +  (ui  -ui^)* 
O'f) 


/W,,k 


2  2 


L‘^l(bii3'*  bi23) +^2^b2i3  +  b^23)J  }  ("d^)  ‘  (2.25) 

{(tu*  -  Wg)[(bj33b2^3  -  b^  ^^b233)  +  <b233bj  33  -  b233b^  23^  ] 

+  (uJ  (b^23^^21  3  ■  ^113^223^  } 


where 


00. 


V  k 
s 


/  D_ 


/Wn;ik 


[5 

I  D, 


2  2 


44 


(2.26a) 

(2.26b) 


/  M  \ 

Note  that  the  term  containing  (-= - j  in  Eq.  (2.  25)  is  of  fourth  power  of 

0^0^ 

the  magnetoelastic  constants  b...  and  usually  negligible.  For  YIG, 
even  near  the  region  where  either  w  —  (longitudinal  crossover)  or 
u)  —  oOg  (shear  crossover)  the  value  of  this  term  is  at  most  5%  of  that  of 
other  terms  at  microwave  frequencies.  From  Eq.  (2.  25),  one  therefore 


has: 
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2  2.,  2  2  2v  /  2  2.  ,  2  .  2  2v  ,2  0^,^ 

W  -uij^)(ui  -C0jg)(w  “Wg)  =  (w  "^s  ^jP  ^  ^  '^jP  ‘‘^s^k^s  (2.27) 


where 


2 


Dm  ^133  ‘*^2'^233^ 

0^0  jP  '"k 


(2. 28a) 


1  ^ 

^2  “^M  1  r  /.  2  ,  ,  2  ,  .  „  2  .  ^2 

^s  ■  D  M  ^  W.  ^“^1^^113  ^123^  ‘^2^‘^213  ^223^^ 

0^0  s  k 


(2.  28b) 


Notice  in  Eq.  (2.27)  that  when  m  ^  uij^,  i.  e. ,  near  the  longitudinal  cross- 

2 

over,  the  term  containing  is  dominant  on  the  right-hand  side,  while 

2 

when  uJ  —  00  ,  i.e.,  near  the  shear  crossover,  the  term  containing  b  is. 

s  s 

2  2 

The  quantities  b.  and  b  are,  therefore,  defined  as  the  longitudinal  and 
Jl  s 

shear  magnetoelastic  couplings,  respectively. 

Now  let  us  proceed  to  maximize  the  longitudinal  magnetoelastic 
2 

coupling  b^  by  choosing  an  appropriate  propagation  direction  with  respect 

to  the  crystallographic  axis  (i.e.,  a  particular  set  of  the  Euler  Angles 

((p, 6, 0 ) where  0  is  a  function  of  (<p, e )  given  in  Eq.  (2.]  4)).  In  general, 

this  problem  is  quite  involved.  HoAvever,  the  simple  approximation 
'•^l 

—  2:  1  may  be  made  for  z-directed  v'aves  in  YIG  at  microwave 

2  00^ 

frequencies.  The  ratio  is  the  square  of  the  ellipticity  of  the  spin 
waves  (Eq.(  2.  17))  which,  for  the  z-directed  waves,  is  due  to  the  aniso¬ 
tropy  field  alone  and  therefore  very  close  to  1.  Using  this  approxi¬ 
mation  one  obtains; 

oOj^  —  Wi  —  ujg  (2.  29) 


2  2 

The  quantities  b^  and  bg  can  be  rewidtten  as; 


±  2  2  . 
Wp  '"'13.3  °233^ 


(2.  30a) 


-  u- 


.2  ^  1  r,,2  ^  ,  2  .  ,.2  ^  .2 

ciJj,  ^^^113  ^123^  ^^213  ^223^^ 


(2.  30b) 


2  1  ,^2  ,^2 


Thus  we  are  left  with  the  problem  of  maximizing  b«  =  rr - —  (b,  oo+b.,„„) 

a  Lj^  1  J  *3  ^  o  o 

which,  according  to  Appendix  II,  (A2.14),  is  independent  of  i.  e.  ,  of 
rotations  around  the  dc  magnetization  (z-axis).  From  (A2.14)  in 

Appendix  II  we  have:  2 

,2  _  /  1  N  .  ^ 

■  \Du  ujJ^ 


2 

^eff 


(2.  31a) 


where 

O 

'-af 


r  2  2  2 

^sin‘'4e+3in  2^(^  sin  6sin2esin46 

2  4  2  6  2  ”1 

+  |sin  2(psin  0sin  20)  +  |- sin  fisin  4(p 


(2.31b) 


A  numerical  analysis  to  maximize  b^^^  in  Eq.  (2.  31b)  by  choosing  angles 
(<p,B)  gives: 


when 


^^eff^max  ^ 

(b„  -  b,)^ 

1 . 26  X  ^ 

4M 

(2.  32a) 

45° 

(2.  32b) 

^®^max 

25.52° 

(2.  32c) 

F’rom  Eqs.  (2.32b)  and  (2.32c),  one  realizes  that  the  propagation 

direction  that  optimizes  the  longitudinal  magnetoelastic  coupling  for 

z  -  directed  waves  is  in  a  {1 1 0]  plane  ((«p)j^^jj  45°)  at  25.  52°  from  a 

<100)  axis;  the  coupling  is  independent  of  the  magnetic  anisotropy 

(within  the  approximation  of  small  spin  wave  ellipticity)  and  directly 

2 

proportional  to  the  (bj  -  bg)  ,  which  is  a  measure  of  magnetoelastic 
anisotropy.  When  b^  -  b2,  corresponding  to  magnetoelastic  isotropy, 
then  the  coupling  vanishes  identically. 
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Let  us  conclude  this  chapter  by  working  out  explicitly  both  the 
longitudinal  and  shear  magnetoelastic  coupling  in  the  {100}  and  {110} 
planes. 

First  of  all,  we  rewrite  the  angle  ^  in  Eq.  (2.  14)  in  terms  of  ((p,9) 
explicitly  from  Eq.  (2.  15)  and  under  the  condition  of  z-directed  waves 
(1<  =  igk)  as; 

tan  2i/,  =  - _  cose  sin  4(p - -  (2.  33) 

4  cos^e  -  sin  2(p(l  +  cos  0) 

Then  we  have  for 


(1)  {100]  plane  i(p  =  0) 


From  Eq.  (2.  33),  =  0,  letting  w  or  to  in  the  appropriate 

places,  and  obtaining  b... 's  from  Appendix  II,  Table  2A.  2;  the  result  is; 

IJK 


2  2 

(w-t0|^)(00-Wjg)(w-t0g)  =  -H(w-t0g)b^  +  (to-to^)bg] 


(2.  34) 


where 


bA<p  =  0,9)  =  ■ 


(2. 35a) 


M 


^2,  _  {b2  +  [^2  ^  ^^1  "  ^2^  /n  ncrux 

b  ((p  —  0,0)  —  n  (2 ,  3 51d) 

<=44  "0“ 


and  use  has  been  made  of; 


and 


L  to  n 
of 


D  to^ 
o  s 


(2.  36a) 


(2.  36b) 


2  2 

A  plot  of  b^(  0,6)  and  bg(0, 0)  is  shown  in  Fig.  2.  3  where  values  of  b^  and 
bg  for  YIG  at  room  temperature  are  used  for  Fig.  2.  3b.  Notice  the 
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2  2  TT  o 

maximum  in  the  {100]  plane  occurs  at  sin‘^46  =  1  or  0  =  ^  ((22.  5°)  from 

2 

a  (100)  axis)  and  the  value  of  in  this  case  is  80%  of  the  absolute 
maximum  (b?)  in  Eq.  (2.  32a).  Because  of  availability,  this  "22.  5°" 
direction  is  chosen  in  the  experiments  described  in  Chapter  III.  For  shear 
magnetoeiastic  coupling,  the  minimum  value  in  the  {100}  plane  occurs  at 

2  ff 

sin  20  =  1  or©  =  ^((110)  direction).  The  numerical  values  of 
(longitudinal  magnetoeiastic  frequency  splittings)  and  b,  (shear  magneto- 


elastic  frequency  splittings)  for  YIG  at  room 

temperature  in  the  "22. 

direction  are 

given  by: 

277b^(0,-^) 

=  5.  53/f 

(MHz) 

(2.  37a) 

2rbJ0l) 

=  45.  5v/f 

(MHz) 

(2.  37b) 

where  f  is  the  frequency  (in  GHz). 


(2)  {110]  plane  (<p  =  J) 


From  Eq.  (2.  33),  ih  =  0,  letting  o)  Wj  or  a)  ^  w,  in  the  appropriate 

Jl  is 

places  and  obtaining  b... 's  from  Appendix  II,  Table  2A.  2;  the  result  is: 


2  2 

(w-u:,,)(oo-ajjg)(u;-oOg)  =  4[(w-c0g)b^  +  (w-w^)bg] 


(2.  38) 


where 


b^((P=f.e>  = 


M 


—  (sin  4  0  +  sin^0sin2  6)^ 


(2.  39a) 


bs(<p=|.e)  = 


'44Mo 


M 


T  ^^^2  f^^l  ’  ^2^  sin^2  0)^ 

+  +  (bj  -  bg)  sin^0)^] 


(2.  39b) 


2  ff  2  ff 

A  plot  of  b^  (416)  and  b^  ^,0)  is  shown  in  Fig.  2.  4,  where  values  of  b^ 


and  bg  for  YIG  at  room  temperature  are  used  for  FTg.  2. 4b.  The 


(TiO)  plane  |CllO]  (TiO)  plane 
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Fig.  2.4  Longitudinal  and  shear  magnetoelastic  couplings  of 
z- directed  waves  in  a  {110}  plane. 
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numerical  values  of  b„  and  b  for  YIG  at  room  temperature  at  0  =  25.  52° 

X. 

are  given  by: 

2TTb^(|.  25.  52°)  =  6.20v/f  (MHz)  (2.40a) 

2nbg(|,  25.52°)  =  36.1v  f  (MHz)  (2.40b) 

Finally,  the  disoersion  of  z-directed  waves  (Eqs.  (2.  34)  or  (2.  37))  is 
shown  in  Fig.  (2.  5),  where  both  longitudinal  and  shear  magnetoeli^stic 
splittings  are  indicated.  Note  from  p'ig.  (2.  4)  or  Eq.  (2.  31b)  that  the 
longitudinal  magnetoelastic  coupling  (or  frequency  splittings)  vanishes 
identically  when  the  propagation  direction  of  z-directed  waves  is  along 
any  of  the  three  principal  crystallographic  directions  (100,  110,  111). 


3 


Fig.  2,  5  Dispersion  diagram  for  z  directed  magnetoelastic 

waves;  the  spin/longitudinal  and  spin/shear  wave 

splittings  (b,  and  b  )  are  shown  very  much 
i.  s 

exaggerated.  The  quantities  v.,  v  and  v,  are, 

il  S  K 

respectively,  group  velocities  of  longitudinal 
elastic,  shear  elastic  and  spin  waves;  the  branch 
numbers  are  referred  to  in  the  text. 
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CHAPTER  III 

EXPERIMENTAL  RESULTS  OF  Z-DIRECTED  MAGNETOELASTIC 
WAVES  USING  MICROWAVE  TECHNIQUES 

In  this  chapter,  microwave  experimental  results  are  reported  for 
both  longitudinal  and  shear  elastic/spin  waves  propagating  in  spatially 
and/or  temporally  varying  magnetic  fields.  The  calculations  made  in 
the  last  chapter  are  used  to  interpret  the  results. 

In  the  past,  a  number  of  authors  have  studied  shear  elastic/spin 

waves  in  YIG  using  either  spatially^  or  temporally  varying  magnetic 
8  18 

fields.  *  Here  similar  studies  on  the  previously-unexplored  longi¬ 
tudinal  elastic/spin  waves  were  carried  out.  (Some  results  have  been 
reported.  '  )  Both  magnetic  injection  from  "fine  wire"  and 

acoustic  injection  from  piezoelectric  transducers  were  employed. 
Variable  delay  and  conversion  of  longitudinal  elastic/spin  waves  in 
spatially  varying  fields  are  reported.  Frequency  shifts  and  pulse 
compression  of  the  longitudinal  elastic/spin  waves  subjected  to  pulsed 
magnetic  fields  have  also  been  obtained.  Finally,  some  preliminary 
experiments  on  shear /longitudinal  elastic  wave  spatial  conversion  via 
spin  waves  are  presented  and  discussed.  Results  of  laser  probing  of 
magnetoelastic  waves  through  infra-red  Bragg  scattering  are  reported 
laster  in  Chaptei-  V. 

In  all  the  experiments  described  here,  either  a  YIG  rod  or  a 
rectangular  bar,  with  the  symmetry  (or  long)  axis  oriented  in  a  {100} 
plane  at  an  angle  of  (22.  5°)  from  a  (100)  direction,  hereon  referred  to 
as  the  "22.  5°"  direction,  was  used.  From  the  results  of  the  last 
chapter  it  is  shown  that  this  crystallographic  direction  provides  nearly 
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80  %  of  tlie  niaxinuini  z-dirccted  longitudinal  magnetoelastic  coupling 
(nearly  90 To  of  the  maxinuim  longitudinal  magnetoelastic  splittings  in 
frequem:y). 

8.  1  Experimental  Setup 

Static  bias  magnetic  fields  were  provided  by  an  electromagnet  and 
the  pulsed  magnetic  fields  were  produced  by  means  of  a  coil  wound 
around  the  sample  and  driven  by  a  current  pulser.  All  experiments 
were  done  at  room  temperature. 

The  L-band  microwave  apparatus  used  is  identical  to  that  of 
20 

RoTende  and  is  shown  in  Fig.  3.1.  It  consists  basically  of  a  signal 
generator  (Alfred  650  main  frame  and  651 A  rf  head,  1-2  GHz)  followed 
by  two  attenuators  (HP  354A  and  Alfred  El 01)  in  tandem,  and  a  PIN 
modulator  (HP  8731 B).  The  PIN  modulator  is  capable  of  producing  rf 
pulses  as  short  as  100  nsec.  A  circulator  (Sperry  D52L21)  was  used 
when  reflection  measurements  were  made.  A  stub  tuner  (Weinschel 
DS109)  which  follows  the  circulator  for  matching  load  impedances 
was  taken  out  of  the  circuitry  when  short  rf  pulses  were  employed.  The 
receiver  system  consists  of  either  a  narrow-band  superheterodyne  or  a 
wide-band  system.  Th(!  former  is  composed  of  a  local  oscillator 
(GR  1021 -P2),  a  balanced  mixer /preamplifier  (LEL  LAC-3-60-12-50, 
bandwidth  12  MHz)  and  a  60  MHz  IF  amplifier  (LEL  IMM-2-7480,  band¬ 
width  10  MHz),  The  Wj^de-band  system  is  formed  by  a  transistor 
amplifier  (Avantck  AM-1000-0,  G  =  25  db,  NF=6  db),  followed  by  a  TWT 
amplifier  (HP  489A,  G  =  35  db,  NF  =  30  db),  and  finally  a  tunnel  diode 
detector  (Aertech  D104B).  The  output  of  the  tunnel  diode  is  amplified 
by  40  db  before  being  fed  into  a  Tektronix  585  scope. 
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Fig.  3  1  Block  diagram  of  L-band  microwave  circuitry. 
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3.  2  Experiments  with  Spin  Wave  In  ections 

The  propagation  of  shear  elastic/spin  waves  in  an  axially  magnetized 
YIG  rod  has  been  studied  theoretically  by  SchlBmann  and  Joseph, 

24  9  10 

Auld  et  al. ,  and  expei  imentally  by  Eshbach  and  Strauss.  The 
internal  field  H^(z)  in  a  YIG  rod  or  a  rectangular  bar  is  nonuniform  and 
given  before  by  Eq.  (1.  3a): 

H.(z)  =  +  Hg”  +  (1.3a) 


&11 

where  H  is  the  externally  applied  field,  H.,  is  the  effective  anisotropy 

field  and  is  the  demagnetizing  field.  The  variation  of  the  latter 

along  the  axis  of  the  rod  or  the  long  axis  of  the  rectangular  bar  is  given, 

2  2 

respectively,  to  first  order  in  M/H^  by  the  Sommerfeld  formula  for 
the  rod; 


H 


dem 


L-z 


[{L-z)2  + 


(3.  la) 


where  z  =  0  is  taken  at  one  of  the  end  faces,  L  is  the  length  of  the  rod, 

23 

and  D  the  diameter,  and  for  the  bar': 


with 


jjdem  ^  _  M[cot"^f(0,  0,  z)  +  cot’^f(0,  0,  -z)] 


f(0,  0,z) 


[(WD)^  +  (HT)^  +  (L-2z)^l(L-2z) 
2(WD)  X  (HT) 


(3.  lb) 


where  z  -  0  is  taken  at  the  center  of  the  long  axis,  L  is  the  length  of  the 
long  axis  of  the  bar  and  HT,  WD  the  transverse  dimensions.  The  axial 
internal  field  profiles  in  each  case  are  shown  in  Fig.  3.2. 

1 2 

3.  2a  Spin/Longitudinal  Elastic  Wave  Conversion 

The  room  temperature  experiments  were  carried  out  at  1. 5  GHz 


-52- 


Fine  wire  _  ^0 


ZnO  Ho 

transducer 


(b)  Axial  internal  field  profile  of  o  YIG  bar 

Fig.  3.  2  Internal  field  profiles  of  a  YIG  rod  and  a  YIG  bar, 

showing  turning  points  (T)  and  crossover  points 

(L.,  and  S  The  distances  between  these 

cros  cros 

points  are  very  much  exaggerated. 
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using  a  single  crystal  YIG  rod  (L  =  0.  99  cm,  0  =  0.3  cm)  with  its 

symmetry  axis  oriented  along  the  "22.  5°"  direction  and  both  end  surfaces 

optically  polished;  fine  wire  excitation  was  used.  The  results  are  best 

explained  with  the  use  of  the  dispersion  curve  shown  in  TTg.  2.  5  and  the 

internal  field  profile  shown  in  Fig.  3.  2a.  Within  a  certain  range  of 

values  of  the  applied  field,  there  exists  a  point  on  the  axis,  the  internal 

field  of  which  corresponds  to  =  0  in  the  spin  wave  dispersion  curve. 

This  is  the  so-called  "turning  point"  marked  T  in  Fig.  3.  2a,  and  is  given 

by  H.(z^p)  =  w/  ;  the  quantity  co  is  the  frequency  of  the  input  micro- 

wave  power.  The  pulsed  rf  magnetic  field  generated  by  the  input  fine 

wire  couples  electromagnetic  energy  to  spin  wave  energy  at  the  turning 

24 

point  via  magnetostatic  waves,  giving  rise  to  a  spin  wave  packet 
(branch  1  with  k  -» 0  in  Fig.  2.  5).  This  wave  packet  then  travels  from 
the  turning  point  toward  the  end  of  the  rod  with  constant  frequency  and 
increasing  wave  number,  accelerating  as  it  encounters  lower  bias  field 
values.  As  the  wave  packet  reaches  the  region  of  coupling  with  longi¬ 
tudinal  elastic  waves  (longitudinal  crossover  point  marked  in 

Fig.  3.2a),  a  fraction  of  the  energy  is  converted  into  phonon  form 
(i.  e.,  stays  on  branch  1  in  Fig.  2.  5)  and  the  rest  stays  in  the  magnon 
state  (branch  2).  According  to  SchlBmann  and  Joseph,  the  fraction  of 
energy  which  stays  in  the  magnon  state  decreases  when  the  ratio  between 

the  field  gradient  at  the  crossover  region  H'(z^^^g)  and  the  critical 

'i  12 

gradient  increases.  The  latter  is  given  by: 


(3.2a) 


Similarly  fcr  the  shear  waves. 


the  critical  gradient  H 


's 

crit 


is  given  by: 
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H 


'  s 

crit 


(3.  2b' 


For  values  of  the  applied  field  such  that  the  longitudinal  crossover 
region  is  near  the  center  of  the  rod  (small  gradient),  most  of  the  magnon 
energy  excited  at  the  turning  point  is  converted  into  longitudinal  phonon 
energy,  which  is  reflected  at  the  end  cf  the  rod  and  therefoi'e  returns  to 
the  turning  point.  The  energy  which  I'emains  in  the  magnon  state  is 
later  converted  almost  completely  (because  of  the  large  shear  critical 
gradient)  into  shear  elastic  form  (at  the  shear  crossover  region  marked 
S  in  Fig.  3.  2a)  to  give  rise  to  another  reflected  echo.  If  the 
applied  field  is  increased  so  that  the  crossover  regions  move  toward  the 
rod  end,  the  spatial  field  gradient  at  the  longitudinal  crossover  region 
increases  and  the  fraction  of  the  energy  which  is  scattered  to  branch  2 
increases,  resulting  in  a  gradual  trade-off  between  the  amplitudes  of  the 
longitudinal  and  shear  wave  echoes.  Over  the  range  of  echo  time  delay 
(0.  8  to  6.  5  sec)  we  have  been  able  to  employ,  the  spatial  gradient  of  the 
field  at  the  longitudinal  crossover  region  varies  from  about  the  value  of 
the  critical  gradient  to  ten  times  larger.  This  allows  the  measurement 
of  the  rnagnon-phonon  conversion  efficiency  ve  -sus  field  gradient  in 
ranges  not  attainable  before.  The  results  are  reported  in  section  3.  3b. 
Typical  data  of  delay  timer  versus  dc  magnetic  field  of  the  first  longi¬ 
tudinal  and  shear  echoes  obtained  at  1.  5  GHz  are  shown  in  Fig.  3.  3 
together  with  the  theoretical  curves  calculated  with  the  assumption  of 
zero  frequency  splittings  at  the  crossover  points.  The  latter  is  given 
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Fig.  3.  3  Measured  time  delay  of  the  first  longitudinal  (LOl) 
and  shear  (SOI)  echoes  versus  dc  magnetic  bias 
fields  at  1 .  5  GHz. 
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where  v  (phonon  velocity)  is  either  v  .  (longitudinal  phonon  velocity)  or 
P  ^ 

Vg  (shear  phonon  velocity)  as  appropriate,  is  the  distance  measured 
from  the  end  surface  to  the  turning  point,  f(=  ^)  is  the  frequency  of  the 
input  microwave  power,  and  H'(z  )  is  the  field  gradient  at  the  turning 

tp 

point  T.  As  can  be  seen,  the  agreement  is  very  good  for  the  longi¬ 
tudinal  elastic  waves. 

Further  endorsement  of  the  interpretation  given  to  the  observed 
echoes  is  provided  by  experiments  with  pulsed  bias  magnetic  fields. 

The  oscillograms  of  Fig.  3.4  illustrate  some  typical  results.  The  lower 
traces  represent  the  video -detected  microwave  pulses  from  the  super¬ 
heterodyne  receiver  and  the  upper  traces  show  the  current  steps  applied 
to  the  coil  wound  around  the  sample,  and  both  traces  are  synchronous. 

In  Fig.  3.4(a)  the  first  longitudinal  and  shear  echoes  can  be  seen  with 
delays  of  2.  0  and  3.  6  sec,  respectively,  corresponding  to  a  static 
applied  field  of  735  Oe.  (Notice  for  this  applied  field  strength,  at  the 
unsaturated  rod  end,  the  magnetoelastic  wave  is  dominated  by  elastic 
wave  characteristics  and  hence  hardly  affected  by  the  unsaturated 
magnetic  medium.  However,  in  the  neighborhood  of  turning  point  and 
crossover  points  where  the  medium  is  essentially  saturated,  the 
magnetoelastic  wave  is  dominated  by  spin-wave  characteristics  and 
hence  the  theory  is  still  applicable.  )  The  insertion  loss  of  the  longi¬ 
tudinal  echo  is  50  db.  Higner-order  echoes  from  longitudinal  and  shear 


branches  which  are  discussed  in  the  next  section  can  also  be  seen  at 


small  delays.  In  Fig.  3.  4(b)  a  10  Oe  positive  field  step  is  applied  when 
the  branch  1  wave  is  on  a  return  trip  as  a  magnon  (i.  e. ,  toward  the 
point  T  in  FTg.  3.  2a)  and  the  branch  2  wave  is  still  in  the  shear  clastic 
state  (i,.  e.,  near  the  end  surface  in  F'ig.  3.2a;  notice  also  that  the 
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Fig.  3. 


4  Magnetoelastic  echoes  affected  by  a  pulsed  magnetic 
bias  field.  The  first  pulse  in  each  lower  trace  is 
due  to  leakage;  the  second  and  third  are,  respectively, 
the  longitudinal  (LOl)  and  shear  (SOI)  magnetoelastic 
echoes  (time  scale  1  ^sec/cm).  A  full  explanation  of 
cases  (a)  -  (e)  is  given  in  the  text. 
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receiver  delays  the  echoes  by  0.  3  ^sec,  so  that  the  scope  traces  show  the 
field  steps  advanced  with  respect  to  the  echo  time  scale).  The  result  is 
that  the  field  step  produces  a  frequency- shifted  echo,  because  the  spin 
waves  subjected  to  a  pulse  field  undergo  a  frequency  shift  with  constant 

g 

wavenumber,  while  the  phonon-like  branch  2  wave  is  virtually  unaffected. 
In  Fig.  3.  4(c)  the  local  oscillator  is  retuned  to  receive  only  the  frequency- 
shifted  echo.  In  Fig.  3.4(d)  the  field  step  is  applied  when  the  branch  1 
wave  is  in  the  longitudinal  elastic  state  (away  from  the  point  in 

Fig.  3.2(a))  and  the  branch  2  wave  is  in  the  magnon  from  (at  the  position 
between  L  and  S  )  resulting  in  a  frequency  shift  of  the  latter.  In 
Fig.  3.4(e),  the  recovery  of  the  frequency- shifted  pulse  is  shown. 


3.2b  Identification  of  Echoes 

In  addition  to  the  two  echoes  due  to  longitudinal  and  shear  elastic/ 

spin  waves  shown  in  Fig.  3.4,  there  are  a  number  of  echoes  whose 

spacings  change  in  an  orderly  manner  as  a  function  of  the  strength  of  the 

2  fi 

applied  field  H^.  This  is  quite  analogous  to  what  Strauss  and  Eggers, 

27 

using  an  experimental  arrangement  similar  to  that  of  Spencer  et  al. , 
have  observed;  they  reported,  in  a  (111)  YIG  rod,  pulse  sequences 
associated  with  both  shear  elastic  and  shear  elastic/spin  wave  propa¬ 


gation.  In  our  case,  however,  there  are  many  more  pulse  sequences 

because  of  the  presence  of  both  longitudinal  and  shear  magnetoelastic 

coupling.  A  typical  echo  pattern  from  reflection  measurements,  using 

the  wide-band  detection  system,  is  shown  in  Fig.  3.  5.  Following 
2  6 

Strauss,  let  T...  -  denote  the  arrival  time  of  reflected  echoes  for  an 
'  ijki 

input  pulse  at  t  =  0.  The  observed  pulses  then  arrive  at  the  times 


'rEL^jTgj^^2kTgs  + 


i  T 


ss 


(3.4) 
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REFLECTED  PULSES 


Fig.  3.5  A  typical  magnetoelastic  echo  pattern  observed  with 
one-port  fine  wire  excitation.  The  first  pulse  in  the 
trace  is  due  to  leakage;  the  other  two  pulses  are  due 
to  the  longitudinal  and  shear  magnetoelastic  coupling 
(time  scale  1  ^sec/cm).  A  full  explanation  is  given 
in  the  text. 
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where 

i  =  0,1.2.  3. . 

j  =  1,2.3, . 

k  =  0,  1,2,  3. . 

-  1,  2,  3, . 


The  quantity  T^g  is  the  time  taken  by  a  shear  elastic  wave  to  travel  the 
length  of  the  rod,  Tgg  the  round  trip  time  of  flight  of  a  shear  elastic/spin 
wave  from  and  back  to  the  turning  point,  and  Tg^  the  corresponding 

times  for  the  longitudinal  waves.  There  are  two  sequences  of  echoes, 
one  originating  from  the  longitudinal  crossover  point,  the  other  the  shear 
crossover  point,  as  indicated  in  Fig.  3.  5  and  Fig.  3.  6.  Each  sequence 
consists  of  a  primary  pulse  (j  =  1  or  ^  =  1)  followed  by  satellite  pulses  of 
smaller  amplitudes  (j,  jf  =  2,  3,  •  •  • ).  As  the  external  magnetic  field  is 
varied  at  a  fixed  frequency,  the  spacing  between  the  primary  pulses 
rerriains  constant  (i.  e. ,  2Tgg  and  2Tg^),  but  the  spacings  between  the 
satellite  pulses  and  the  primary  pulses  in  different  sequences  change. 
Indeed,  the  behavior  of  the  pulse  spacings  as  a  function  of  permits  the 
identification  of  the  various  echoes.  In  particular,  in  the  sequence 
originating  from  the  shear  crossover  point,  tJie  echo  (0110)  (the  dashed 
line  in  Fig.  3.  6)  starts  as  a  shear  elastic/spin  wave  from  the  shear 


crossover  point  S  ,  converts  to  a  longitudinal  elastic  wave  at  the 
longitudinal  crossover  point  ori  its  return  trip  from  the  end  surface, 

and  finally,  after  traversing  the  length  of  the  rod,  returns  back  to  the 
turning  point  (Fig.  3.2(a)).  Since  the  internal  field  profile  is  concave 
downward,  the  corresponding  echo  (1001)  does  not  exist. 


The  echo  sequence  (OOkjf)  can  be  interpreted  as  due  to  incomplete 
conversion  between  branch  1  and  branch  2  in  Fig.  2.  5  at  the  longitudinal 


PULSE  AMPLITUDE 
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REFLECTEO  PULSES 
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KIr.  3. 6  An  echo  pattern  of  Fig.  3,  5,  reconstructed  from 
theoretical  conaiderations. 
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crossover  point.  The  echo  sequence  (ijOO)  corresponds  to  what  Strauss 
26 

and  Eggers  have  reported.  While  they  suggested  that  the  origin  of 
the  sequence  they  observed  is  due  to  the  polarization  reversal  of  the 
clliptically  polarized  shear  elastic  waves  at  the  end  surface  of  the  (111) 
YIG  rod,  we  do  not  understand  very  well  about  the  origin  of  the  sequence 
(ijOO)  we  observed  in  the  "22.  5°"  YIG  rod. 


3.  3  Experiments  with  Acoustically  Injected  Longitudinal 
Elastic /Spin  Waves 

3.  3a  Propagation  of  Longitudinal  Elastic /Spin  Waves 
in  Spatially  Varying  Magnetic  Fields 

The  room  temperature  experiments  were  performed  using  a  YIG 


bar  which  measured  4  mm  x  4  mm  x  12  mm  and  had  its  long  axis 
oriented  along  the  "22.  5°"  direction  mentioned  earlier.  A  longitudinal 


elastic  wave  packet  at  1. 5  GHz,  100  nsec  duration,  was  launched  using  a 

c 

sputtered  ZnO  thin  film  piezoelectric  transducer.  For 
where  is  the  applied  field  value  for  which  the  longitudinal  crossover 

point  L  in  Fig.  3.  2(b)  is  at  the  middle  of  the  long  axis  (z  =  0  in  this 
case),  the  wave  character  becomes  more  spin-like  as  the  wave 


propagates  toward  the  center  of  the  bar.  The  delay  difference  “^SL 
between  the  longitudinal  elastic/spin  waves  and  the  longitudinal  elastic 
wave  can  be  calculated  by: 


where  the  group  velocity  of  the  longitudinal  elastic/spin  wave  v  (of 

o 

branch  1  in  Fig.  2.  5)  can  be  derived  from  the  dispersion  relation  in 
Eq.  (2.  34).  Numerical  integration  of  (3.  5)  for  the  YIG  bar  yields  the 
curve  for  the  delay  difference  and  group  velocity  versus  applied  magnetic 
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fiolds  shown  in  Fig.  3.  7.  Notice  the  extremely  small  range  of 
i-oquircd  for  the  cliange  of  wave  character  from  longitudinal  elastic  to 
spin-like,  compared  to  that  of  shear  waves.  This  is  expected  because 
the  magnetoelastic  coupling  of  the  longitudinal  waves  is  approximately 
two  orders  of  magnitude  smaller  than  that  of  the  shear  waves.  Experi¬ 
mentally  it  was  also  observed  that  the  longitudinal  elastic/spin  waves, 
unlike  the  shear  ones,  did  not  experience  appreciable  changes  of  the 
delay  time  because  the  change  of  group  velocity  is  too  abrupt  to  be  easily 
observed  (Fig.  3.  7).  In  fact,  the  wave  packets  only  showed  some  wave¬ 
form  distortion  because  of  dispersion,  as  shown  in  Fig,  3.  8.  In 
Fig.  3.  8(g)  and  the  lower  traces  of  Fig,  3.  8(h),  (i)  and  (j),  the  input 

power  has  been  increased  by  13  db  to  show  the  distorted  wave-form. 

H 

The  quantity  pertaining  to  Fig.  3.  8  is  estimated  to  be  716  Oe  and  is 

D 

within  6%  of  the  theoretical  calculation.  For  2  H^,  i.  e. ,  when  the 

longitudinal  crossover  point  is  inside  the  bar,  strong  Bragg  infrared 

light  scattering  from  spin  waves  has  been  observed  and  will  be  described 

0 

in  Chapter  V.  For  >  H^,  branch  hopping  between  the  upper  branch 
(branch  1  in  Fig.  2.  5)  and  lower  branch  (branch  2  in  Fig.  2.  5)  takes 
place  (Fig.  3.  8,  (h)  to  (p).  We  have  taken  advantage  of  this  fact  to 
measure  the  spatial  spin/elastic  wave  conversion  to  be  described  in  the 
next  section. 

3.  3b  Spatial  Longitudinal  Elastic/Spin  Wave 
Conversion  Efficiencies^^ 

In  this  section  we  present  the  results  of  direct  measurements  of  the 
longitudinal  elastic/spin  wave  conversion  efficiencies  at  room  tempera¬ 
ture  in  a  single  crystal,  axially  magnetized  YIG  rod;  the  rod  used  here 
is  identical  to  the  one  described  in  section  3.  2. 
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Fig.  3.  7  Delay  of  elastically  injected  longitudinal  elastic/ 
spin  wave  packet  and  group  velocity  at  the  bar 
center  as  a  function  of  applied  magnetic  field  for 
a  "22.  5°"  YIG  bar  which  measures  4  mm  x  4  mm 


X  12  mm. 


Vo/V, 


66  7  Oe 
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transducer,  as  a  function  of  external  applied 
magnetic  field  (time  scale  0. 1  /^ec/cm). 


-66- 

Although  shear  elastic/spin  wave  conversion  has  been  studied 

11  21 

theoretically  by  previous  authors,  '  '  experimental  investigation  has 

10  25  28 

been  limited  to  the  case  of  nearly  complete  conversion.  '  '  This 

I  g 

is  due  to  the  fact  that  critical  field  gradients  of  shear  elastic/ spin 

waves  (Eq.  (3.  2b))  in  YIG  are  much  larger  than  those  obtainable  in 
practice,  while  that,  of  the  longitudinal  magnetoelastic  splitting 

(Eq.  (3.2a))  are  much  smaller  and  well  within  the  reach  of  our  experi¬ 
ments. 

We  have  used  two  different  ways  to  measure  the  conversion 
efficiencies,  one  with  magnetic  injections  in  which  only  the  weak  coupling 
region  can  be  investigated,  and  another  with  acoustical  injections,  in 
which  case  nearly  the  entire  range  of  interactions  was  covered;  the 
results  in  both  cases  are  in  reasonable  agreement  with  the  previous 
theory. 


I.  Experiments  with  Magnetic  Injections 

The  experimental  arrangement  is  identical  to  that  described  in 
section  3.  2a.  The  mechanism  of  the  wave  propagation  has  also  been 
described  in  detail  there.  It  is  sufficient  to  say  here  that  when  the  spin 
wave  originating  at  a  turning  point  first  reaches  the  longitudinal  cross¬ 
over  point  (marked  L  in  Fig.  3.2a),  a  fraction  of  the  energy  is 
converted  into  phonon  form  and  the  rest  remains  in  the  magnon  state. 

The  conversion  efficiency  was  first  studied  by  Schltimann  and  Joseph^  ^ 

21 

and  later  also  by  Kirchner  et  al.  Their  results  indicate  that  the 
conversion  efficiency  tj  is  related  to  the  field  gradient  H'(z^^)  at  a 
crossover  point  in  the  following  manner; 


»'crlt 


ri 


1  -  e 


(3.6) 
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whcr('  is  the  critical  field  gradient  for  either  longitudinal  and  shear 

waves  (Eqs.  (3.  2))  and  is  the  field  gradient  at  either  longitudinal 

ft  S 

(z  )  or  shear  crossover  point  (z  ).  In  our  YIG  rod,  at  1. 5  GHz  and 
cr  cr 

room  temperature,  the  critical  gradients  for  longitudinal  and  shear 
waves  are,  respectively,  225  and  28,  500  Oe/cm.  The  trade-off  between 
the  amplitudes  of  the  longitudinal  and  shear  wave  echoes  described 
earlier  in  section  3.  2a  can  best  be  described  in  terms  of  longitudinal 
elastic/spin  wave  conversion  efficiency  r/  in  the  following  equation: 

"  p7  *=^Pt2(Tgg  -  Tg^^)  3^3  ]  (3.7) 

TJ  L 


where  P3/P]_^  is  the  power  ratio  between  shear  and  longitudinal  elastic/ 

spin  wave  echoes,  2Tgg  and  2Tg^  are,  respectively,  the  total  time 

spent  in  the  magnon  state  for  shear  and  longitudinal  elastic/spin  waves, 

and  a  is  the  magnon  loss  (in  db)  per  unit  time.  We  have  neglected  the 

phonon  loss  which  is  very  small.  Figure  3.  9  shows  the  measured 

results,  using  a  =  12  db//isec  (full  spin-wave  linewidth  0. 16  Oe)  to  fit 

21 

the  theoretical  curve.  Here  the  field  gradients  have  been  determined 
using  the  well-known  Sommerfeld  formula  in  Eq.  (3.  la).  Although  the 
field  strengths  applied  in  the  experiment  were  fairly  low  (from  650  to 
1200  Oe)  the  assumption  that  the  magnetization  of  the  YIG  rod  is 
saturated  by  the  applied  fields  is  in  reasonable  agreement  with  the 
argument  given  previously  in  section  3.  2a.  Note  that  the  experimental 
points  fluctuate  about  the  theoretical  value.  The  fluctuation  of  the 
amplitudes  of  the  magnetoelastic  echoes,  both  longitudinal  and  shear, 
has  been  observed  experimentally,  but  its  origin  is  not  yet  understood. 


The  trade-off  between  the  amplitudes  of  the  longitudinal  and  shear 
elastic/spin  wave  is  vividly  demonstrated  in  Fig.  3.9.  .  However,  the 
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magnetic  field 
injection. 
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oxcessive  magnon  loss  prevents  us  from  going  to  the  strong  coupling 
region  in  this  type  of  experiment,  because  a  small  spatial  field  gradient 
(therefore  a  strong  coupling)  causes  excessive  attenuation  due  to  long 
magnon  delays.  Nevertheless,  we  have  managed  to  measure  the 
conversion  efficiency  in  a  strong  coupling  region  by  using  acoustically 
injected  longitudinal  elastic/spin  waves. 

II.  Experiments  with  Acoustically  Injected  Longitudinal 
Elastic/Spin  Waves 

In  order  to  obtain  the  conversion  efficiency  throughout  the  entire 
range  of  the  interaction,  experiments  using  acoustical  injections  were 
carried  out.  The  room  temperature  experiments  were  also  performed 
at  1.  5  GHz  using  a  ZnO  transducer  bonded  to  one  of  the  two  polished  rod 
end  surfaces.  The  longitudinal  acoustic  echoes  were  observed  as  a 
function  of  the  strength  of  applied  field  H^.  In  both  the  low  and  high 
magnetic  field  values,  where  the  magnetoelastic  interaction  is  small, 
the  amplitude  of  the  echo  is  marked.  As  the  magnetic  field  increases 
from  zero  to  a  point  where  the  longitudinal  crossover  region  is  placed  in 
the  middle  of  the  rod,  i.  e. ,  =  H^,  the  elastic  echo  disappears 

(or  rather  is  attenuated  by  more  than  30  db)  as  a  result  of  (almost) 
complete  conversion  into  spin  waves.  A  further  increase  of  the 
magnetic  field  moves  the  crossover  region  toward  the  end  of  the  rod, 
thus  also  increasing  the  spatial  field  gradients.  The  longitudinal 
elastic /spin  wave  conversion  efficiency  starts  to  decrease,  as  observed 
by  the  increase  of  the  amplitude  of  the  longitudinal  elastic  wave  echoes 
(refer  to  traces  (h)  to  (p)  in  Fig.  3.  8).  In  addition,  very  weak  spin/ 
longitudinal  elastic  echoes  were  also  observed.  Therefore,  by 
measuring  the  attenuation  of  the  first  longitudinal  elastic  echo  as  a 
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function  of  the  magnetic  field  strength,  one  obtains  the  conversion 
efficiencies  in  nearly  the  entire  range  of  the  interaction.  Since  there 
are  actually  two  longitudinal  crossover  points  in  an  axially-magnetized 
rod,  a  longitudinal  elastic  disturbance  makes  conversion  to  spin  waves 
at  four  occasions  in  making  a  round  trip  throughout  the  rod.  The 
experimental  results  are  plotted  in  Fig.  3. 10,  where  use  has  been  made 
of 

=  P.  (1  -  (3.  8) 

Again  the  field  gradients  have  been  determined  using  the  Sommerfeld 
formula  in  Eq.  (3.1a).  The  quantity  P^  is  the  power  of  the  longitudinal 
elastic  waves  in  the  absence  of  the  magnetoelastic  interaction  and  P^  is 
the  one  with  interaction.  The  result  indicated  in  Fig.  3.10  is  in  fair 
agreement  with  the  theory. 

3.  3c  Propagation  of  Longitudinal  Elastic /Spin  Waves 
in  Temporally  Varying  Magnetic  Fields 

Frequency  shifts  and  pulse  compression  of  the  wave  packets  subject 
to  pulsed  magnetic  fields  are  discussed  below. 

Measurements  of  the  frequency  shifts  of  longitudinal  elastic/spin 
waves  subjected  to  both  positive  and  negative  pulsed  fields  were  made 
using  the  axially  magnetized  YIG  bar  described  earlier  in  section  3.  3a. 
A  polycrystal  YIG  bar  with  identical  saturated  magnetization  was 
attached  to  one  end  of  the  single  crystal  YIG  bar  to  reduce  the  field 
gradients  inside  the  sample,  while  a  thin  film  ZnO  transducer  was 
sputtered  on  the  other  polished  end  surfaces.  Positive  frequency 
shifts  of  up  to  3  MHz,  and  negative  frequency  shifts  of  down  to  1  MHz 
have  been  observed  with  a  95  Oe/Amp  coil  wound  around  the  composite 
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efficiency  versus  inverse  spatial  magnetic  field 
gradient  obtained  with  acoustic  wave  injection. 
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structurc.  The  frequency  shifts  were  observed  with  a  Tektronix  IL  20 

spectrum  analyser  used  to  study  the  first  longitudinal  elastic  echo  which 

was  gated.  The  small  longitudinal  magnetoelastic  coupling  makes  the 

frequency  modulation  effect  of  the  spatially  non-uniform  internal  fields 

(with  this  composite  structure)  much  stronger  than  in  the  case  of  shear 

waves.  Consequently  no  very  meaningful  quantitative  data  have  been 
2Q 

obtained. 

O 

Pulse  expansion  and  pulse  group  velocity  modulation  with  the 
longitudinal  clastic/spin  again  was  not  successful.  The  failure  is 
believed  to  be  due  to  the  same  causes  that  complicate  the  frequency  shift 
measurements,  as  stated  above,  i.  e. ,  non-uniform  internal  fields  and 
the  comparative  smallness  of  the  longitudinal  magnetoelastic  coupling. 
Pulse  compression  of  up  to  10  db  (considerably  larger  than  with 

O 

shear  waves  )  has  been  achieved  when  the  strength  of  the  applied  field 
is  either  H  <  H  or  H  >  H  .  The  mechanism  of  pulse  compression 

w  C  w  C 

with  acoustic  wave  injections  can  be  described  in  the  following  manner. 

A  longitudinal  elastic  wave  packet  of  spatial  length  (and  duration 
AT.  =  ALj/Vjg)  is  generated  by  the  ZnO  transducer  at  the  rod  end.  As 
it  propagates  toward  the  center  of  the  rod  of  increasing  magnetic  fields 
at  fixed  frequency  but  variable  wavenumbers,  its  character  becomes 
more  spin-like,  its  spatial  length  shortens  because  of  its  decreasing 
group  velocity,  but  its  duration  ATj  remams  unchanged.  When  a 
negative  pulsed  field  that  opposes  the  dc  magnetic  field  is  applied,  then 
the  duration  of  the  wave  packet  suddenly  decreases  (hence  pulse 
compression)  as  its  group  velocity  increases  abruptly  to  that  of 
longitudinal  elastic  waves  v^j.  The  duration  AT^  of  the  compressed 
pulse  can  be  found  according  to; 
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(3.9) 

wluM'i*  is  tiu*  fii  oup  velocity  of  the  elaKtic/spin  wave  packet  at  the  time 
thi‘  pulsiril  field  is  applied.  Kquation  (3.9)  can  also  be  derived  from 
eonsidei  ation  o'  the  Fourier  spectra  of  the  wave  packet  in  wavenumbers 
and  frequencies.  Notiee  that  Eq.  (3.9)  indicates  that  the  pulse 
compression  ratio  (ATj/^Tj.)  is  always  larger  for  the  case  of  longi¬ 
tudinal  waves  than  for  the  shear  wave  case. 

A  typical  case  in  which  >  11^  is  shown  in  Fig.  3. 1 1 .  The  lower 
trace  of  Fip.  3. 1  l(n)  indicates  a  longitudinal  elastic  echo  which  has  gone 
through  both  crossover  points  inside  the  sample.  The  upper  trace  of 
Fig.  3. 1 1(a)  indicates  a  negative  field  pulse  of  80  Oe  that  opposes  the 
dc  bias  field.  In  Fig.  3. 1 1(b),  the  field  pulse  is  applied  after  the  wave 
train  has  been  reflected  from  the  free  end  of  the  bar  and  returned  to  the 
crossover  point  between  the  center  and  the  free  end  of  the  bar.  In  this 
case  only  the  trailing  edge  is  compressed  since  the  trailing  edge  is  still 
spin- like  (on  branch  1  in  Fig.  2.5),  while  the  leading  edge  has  hopped 
(to  branch  2  in  Fig.  2.  5)  and  is  therefore  elastic-like.  Similarly  in 
P'ig.  3. 1 1  (c),  the  field  pulse  is  applied  when  the  waves  have  left  the 
transducer  for  the  first  time  and  reached  the  crossover  point  between 
the  transducer  end  and  the  center  of  the  bar.  Again  the  trailing  edge 
(on  branch  1  in  P^ig.  2.5)  is  more  spin- like  than  the  leading  one  (on 
branch  2  in  Fig.  2.5)  and  therefore  is  compressed.  This  kind  of 
compression  again  illustrates  the  small  longitudinal  magnetoelastic 
coupling  and  the  comparatively  larger  effect  of  non-uniform  internal 
fields  on  the  longitudinal  waves  as  compared  to  the  shear  waves. 

When  a  typical  case  is  shown  in  Fig.  3.12.  In 
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(C) 

Fig.  3.11  Pulse  compression  of  longitudinal  elastic/spin 

waves  subjected  to  pulsed  magnetic  fields.  The 
static  magnetic  field  here  is  high  enough  so  that 
there  are  crossover  points  inside  the  YIG  bar 
(i.  e.,  >  H^).  Horizontal  sweeps  are  from 

left  to  right  with  scale  0.  5^sec/cm  (frequency 
1.  5  GHz). 
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(a) 
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(b)  Ho  «  713  Oe 
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(c) 


Fig.  3. 12  Pulse  compression  of  longitudinal  elastic/spin 
waves  subjected  to  pulsed  magnetic  fields,  when 
Hq  <  H^.  Horizontal  sweeps  are  from  left  to 
right  with  scale  0.  5^sec/cm  (frequency  1. 5  GHz). 
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Fig.  3.12(a),  the  upper  trace  indicates  the  first  longitudinal  elastic  echo 
(marked  with  an  arrow)  which  emerged  after  passing  through  the  length 
of  the  bar  in  the  absence  of  longitudinal  magnetoelastic  coupling,  and  the 
lower  trace  shows  the  negative  field  pulse  of  10  Oe  that  opposed  the  dc 
bias  field.  In  Fig.  3. 12(b),  is  increased  so  that  the  longitudinal 
elastic  echo  (marked  with  an  arrow)  has  been  strongly  attenuated.  In 
Fig.  3. 12(c),  the  negative  field  pulse  is  applied  just  as  the  wave  packet 
reaches  the  middle  of  the  rod  for  the  first  time;  therefore  the  longi¬ 
tudinal  echo  is  compressed  as  a  whole  for  a  compression  ratio  of  10. 

3.  4  Preliminary  Results  on  Longitudinal /Shear  Elastic  Wave 
Conversion  via  Spin  Waves  in  Spatially  Varying  Fields 

30 

Comstock  and  Wigen  first  proposed  the  spatial  conversion  from 

shear  to  longitudinal  elastic  waves,  and  vice  versa,  via  spin  waves  in 

31 

transversely  magnetized  rods.  Morgenthaler  also  suggested  longi¬ 
tudinal/shear  elastic  conversion  via  spin  waves  by  means  of  time-varying 
magnetic  fields.  Experimental  studies  of  these  two  cases  have  not  yet 
appeared  in  the  literature.  Here  we  report  some  preliminary  results  on 
spatial  longitudinal/shear  elastic  conversion  via  spin  waves  in  axially 
magnetized  single  crystal  YIG  rod  oriented  along  the  "22.  5°"  direction. 
The  room  temperature  experimental  arrangement  is  shown  in  Fig.  3.  13, 
where  a  polycrystalline  nickel  ferrite  rod  with  a  saturation  magnetization 
of  3000  Gauss  is  placed  at  one  end  of  the  YIG  rod  to  obtain  a  mono- 
tonically  increasing  internal  field  profile.  A  decreasing  field  profile, 
however,  exists  at  the  end  of  the  YIG  rod  near  the  air  gap  between  the 
rods.  A  CdS  thin  film  on  each  end  of  the  rod  served  as  an  input  shear 
wave  and  output  longitudinal  elastic  wave  transducer;  the  input  trans¬ 
ducer  unfortunately  does  not  have  good  discrimination  against  the 
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excitation  of  longitudinal  elastic  waves  and  therefore  generates  both 
shear  and  longitudinal  waves  at  the  same  time.  A  shear  elastic  wave 
packet  (and  also  a  longitudinal  one),  launched  at  the  input  transducer, 
travels  toward  the  center  of  the  rod.  When  it  reaches  the  shear  cross¬ 
over  point  (marked  in  Fig.  3.13  and  corresponding  to  branch  2  in 

Fig.  2.  5),  it  is  nearly  completely  converted  into  spin  waves;  the  latter 
then  proceed  to  the  longitudinal  crossover  (marked  in  Fig.  3. 13), 

whereupon  part  of  the  energy  is  converted  into  longitudinal  elastic  waves 
(corresponding  to  staying  in  branch  2  in  Fig.  2.  5),  which  then  move 
forward  and  reach  the  other  transducer.  The  delay  timer^  can  be 
deduced  from  Eq.  (3.  3)  as: 
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where  Tj,  is  the  time  duration  spent  in  the  elastic  wave  state,  while  Tg  is 
the  time  duration  spent  in  the  spin  wave  state.  The  conversion 
efficiency  can  be  written  as; 
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where  the  phonon  loss  is  neglected  and  a  is  the  magnon  loss  (in  db/nsec)t 
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transmitting  receiving 


Fig.  3. 13  Arrangement  used  to  observe  spatial  conversion 
from  shear  elastic  waves  to  longitudinal  elastic 
waves  via  spin  wave  coupling. 


'  ’  s 

H  the  field  gradient  at  the  longitudinal  crossover  point,  H  the  corres¬ 
ponding  shear  one.  Both  the  delay  times  and  the  conversion  efficiencies 

'St 

as  a  function  of  H  are  plotted  in  Fig.  3. 14,  using 


a  =  5  db/^sec 


(3.  13) 


and  the  values  of  the  critical  field  gradient  are  given  in  Eq.  (3.2). 

Notice  that  as  the  external  field  (and  hence  the  field  gradient)  increases, 
the  delay  time  approaches  the  time  required  for  the  elastic  waves  to 
traverse  the  length  of  the  rod.  Also  there  is  a  peak  in  the  conversion 
efficiencies  corresponding  to  an  optimum  H'  (or  H^).  The  experi¬ 

mental  observations  are  shown  in  Fig.  3. 15.  Figure  3,  15(a)  shows  the 
reflection  echoes  from  the  input  transducer,  where  echoes  marked  L 
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as  a  function  of  field  gradient  at  the  longitudinal 
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(c)  554  Oe 
L  S^l  L 


(e)  627  Oe 


(g)  756  Oe 
U  2  *  L 


(h)  873  Oe 
L  S  L  L 


( i)  1125  Oe 


Fig.  3.15  Oscillograms  showing  the  spatial  conversion  from 
shear  elastic  waves  to  longitudinal  elastic  waves 
via  spin  wave  coupling.  The  horizontal  traces 
are  from  left  to  right  (time  scale  1  ^sec/cm). 
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are  longitudinal  elastic  ones  while  the  echo  marked  S  is  the  shear  one. 

In  Fig.  3. 15,  from  (b)  through  (i),  are  the  echo  patterns  picked  up  at  the 
receiving  transducer  as  a  function  of  the  external  applied  field  H^.  All 
the  transmitted  pulses  marked  either  L  or  S  are  the  purely  elastic  longi¬ 
tudinal  or  shear  echoes.  In  Fig.  3.15(c),  the  crossover  points  come  in 
the  center  of  the  rod  and  therefore  almost  all  echoes  except  that  corres¬ 
ponding  to  the  shear  one  disappear  as  a  result  of  conversion  to  spin 
waves.  In  Fig.  15,  (d)  through  (e),  there  are  two  echoes  marked  by 
arrows;  the  first  one  (with  1  on  the  arrow)  is  believed  to  come  from 
shear/longitudinal  conversion  and  does  show  a  maximum  amplitude  in  (f) 
as  predicted  in  Fig.  3. 14;  the  second  echo  (with  2  on  the  arrow),  which 
is  always  spaced  at  exactly  one  round  trip  time  required  for  longitudinal 
elastic  waves,  is  thought  to  come  from  the  reflection  of  the  first  echo 
at  the  receiving  transducer.  In  addition,  notice  from  (h)  and  (i)  that  the 
first  echo  (with  arrow  marked  1 )  eventually  approaches  the  transmitted 
longitudinal  echo  and  these  two  echoes  disappear  at  the  same  time  when 
the  fields  are  increased  further,  as  in  (i).  Notice  that  the  decreasing 
field  profile  at  the  receiving  transducer  end  merely  inflicts  some 
additional  loss.  The  loss  is  caused  by  the  conversion  from  the  trans¬ 
mitted  longitudinal  waves  to  spin  waves  at  the  longitudinal  crosspoint 
near  the  receiving  end. 

We  are  unable  to  make  any  quantitative  comparison  with  the 
theoretical  curve  in  Fig.  3. 14  because  of  the  unknown  internal  field 
profile  for  this  composite  structure  and  the  poor  qualities  of  the 
transducers.  However,  we  do  not  believe  that  these  two  echoes  could 
come  from  the  input  longitudinal  elastic  waves  for  the  following  reason. 
In  order  to  have  variable  delay  times  as  a  function  of  the  external  fields. 
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the  longitudinal  wave  would  have  to  make  a  round  trip  to  and  from  the 
turning  point  plus  a  one-way  trip  as  a  longitudinal  elastic  wave  from  the 
transmitting  transducer  to  the  receiving  one.  In  this  case  the  waves 
would  spend  more  time  as  spin  waves,  go  through  three  times  the  longi¬ 
tudinal  conversion  point  and  therefore  suffer  more  loss  than  the  shear/ 
longitudinal  elastic  waves  mentioned  earlier.  Furthermore,  we  do  not 
expect  to  see  a  maximum  amplitude  of  the  echo,  such  as  in  Fig.  3.  15(f), 
for  this  case.  For  all  these  reasons,  we  believe  the  first  echo,  whose 
delay  time  varies  with  the  external  applied  field,  does  come  from  shear/ 
longitudinal  elastic  wave  conversion  via  spin  waves. 
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CHAPTER  IV 

THEORY  OF  LIGHT  DIFFRACTION  FROM  COHERENT 
ELASTIC  AND  SPIN  WAVES 


4. 1  Introduction 

Scattering  of  light  by  means  of  elastic  waves  in  liquids  and  solids 

has  been  extensively  studied  since  the  effect  was  first  predicted  by 
32  33 

Brillouin  in  1922  and  observed  by  Gross  in  1930  in  liquids  where  only 

34 

longitudinal  elastic  vibrations  are  allowed.  In  1932,  Debye  and  Sears 

35 

and,  independently,  Lucas  and  Biquard,  performed  the  first  succer-isful 

diffraction  of  light  from  coherent  ultrasonic  waves.  This  phenomenon, 

light  inelastically  scattered  from  an  acoustic  elastic  wave,  has  come  to 

be  known  as  Brillouin  scattering,  while  light  inelastically  scattered  from 

36 

optical  phonons,  discovered  in  1928  by  C.  V.  Raman,  is  known  as 

37 

Raman  scattering.  Early  work  on  Brillouin  scattering  was  used  to 

determine  elastic  constants  and  to  study  thermal  phonon  distributions  in 

various  materials.  The  advent  of  the  laser  in  the  sixties  has  greatly 

revived  the  intei  est  in  light  scattering  and  made  the  technique  useful  to 

probe  not  only  elastic  waves  but  many  other  elementary  excitations  such 

as  plasma  waves  and  spin  waves. 

For  Brillouin  scattering,  it  is  necessary  to  distinguish  between 

diffraction  at  the  Bragg  angle  (Bragg  scattering)  and  Raman- Nath  type 

3  8 

scattering,  after  Raman  and  Nath  who  explained  the  spatial  form  of  the 
scattered  light  in  the  experiments  of  Debye  and  Sears  and  Lucas  and 
Biquard.  In  Raman-Nath  type  scattering  the  light  enters  the  medium 
parallel  to  the  wavefronts  of  the  elastic  waves  and  the  diffracted  light 
appears  on  both  sides  of  the  primary  beam  in  the  form  of  equally  spaced 
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liiit?s.  \\  ilhiril  ''  dis('iiss(><l  the  distinction  between  these  two  cases  and 

showi'd  that  siailliM’itig  I’esults  whenever  the  acoustic  column  is 

suffii;ienlly  u  ide  such  that  the  product  of  the  optical  wavelength  in  the 

iiiediiim  (\/n)  and  the  acoustic  beamwidth  (w)  is  greater  than  the  square 

of  the  acoustic  wavelength  (A)  (i.  c. ,  ^w  > /i).  Indeed,  Rytov^^  has 

shown  that  the  diffracted  light  will  appear  in  a  single  spot  and  will  be  a 

maximum  when  the  light  is  incident  at  the  Bragg  angle.  This  is  also 

indicated  in  Brillouin's  original  work,  and  it  is  the  only  case  with  which 

this  thesis  will  be  concerned.  There  have  been  many  excellent  review 

41-48 

articles  on  Brillouin  scattering,  on  its  theories,  its  application  as 

49-5') 

a  research  tool,  '  ‘  and  its  practical  applications  such  as  in  optical 

51-55 

signal  processing. 

In  this  chapter,  the  main  objective  is  to  review  the  theory  of  light 

scattering  and  generalize  it  to  include  the  theory  of  coherent  light 

scattering  from  coherent  spin  waves.  The  generalized  theory  presented 

56 

here  has  essentially  revised  Auld  and  Wilson's  (A-W)  calculations  by 
using  the  perturbed  dielectric  tensor,  which  is  linear  in  small  signal 
spin  ave  amplitudes  rather  than  in  total  crystal  magnetization.  It 
appears  to  explain  adequately  not  only  our  experimental  results 

57  58 

doscriV'd  in  Chapter  V,  but  also  observations  that  previous  authors  ' 
reported  as  "anomalous,  "  unexplained  by  A-W  theory. 

Our  now  theory  indicates  that  the  characteristics  of  light  scattering 
from  coherent  spin  waves  are  quite  different  from  scattering  from 
elastic  waves.  Unlike  the  case  of  coherent  elastic  waves,  the  scattered 
light  intensity  from  coherent  spin  waves  is  found  to  bo  dependent  on  both 
the  incident  light  polarization  and  the  sign  of  the  frequency  shift  (Stokes 
or  anti-Stokes  lines).  Also  the  intensity  is  independent  of  the  direction 
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of  the  applied  magnetic  field. 

The  chapter  begins  with  a  rather  extensive  review  of  the  macroscopic 
theory  of  light  scattering,  followed  by  a  description  of  the  characteristics 
of  light  scattering  from  coherent  elastic  waves.  Finally,  the  theory  of 
coherent  spin  wave  light  scattering  is  presented. 

Because  of  convention  and  convenience  in  comparing  with  the  optics 
literature,  the  c.g,  s.  unit  system  is  used  throughout  this  and  the  next 
chapter. 

4.2  The  Macroscopic  Theory  of  Light  Diffraction 

Before  starting  our  discussion,  a  few  words  regarding  magnetooptics 
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in  general  are  appropriate  here.  Following  Landau  and  Lifshitz, 
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Pershan  discussed  in  great  detail  that  fact  that,  especially  at  optical 
frequencies,  the  propagation  of  an  electromagnetic  wave  in  any  material 
can  always  be  described  in  terms  of  an  effective  dielectric  permittivity 
C  and^  =  1.  This  c  •  which  incorporates  whatever  magnetic 
characteristics  a  particular  material  has,  can  be  defined  in  the  following 
manner.  Maxwell's  equations  in  c.g.  s.  unit  system  for  a  nonconducting 
and  current-free  medium  are- 


n  V 

^  _  1  ai!  47r 

(■  .la) 

V  A 

c  at  c  at 

7  X 

jt  _  1  a^  4ff  a^ 

c  at  c  at 

(4.1b) 

7  • 

(2+  4ff?)  =  7.  (c2)  =  0 

(4. Ic) 

7  • 

(il+4ffi^)  =  0 

(4.  Id) 

where  i^,  2  are  the  magnetic  and  electric  field  vectors.  The  vectors  ^ 
and  ?  are  the  magnetization  and  polarization  of  the  medium  and  c  is  the 
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velocity  of  light  in  free  space;  the  quantity  c  is  the  dielectric  constant  of 
the  medium. 

From  Eqs.  (4.1a)  and  (4.1b),  one  obtains 


1 

2 

c 


^  -  (7x7x5) 


c2  at 


(4.2a) 


where 


eff 


c  7  X  M 


(4.2b) 


Equations  (-1.  2)  indicate  that  the  effect  of  the  magnetization  can  be 

-♦  eff 

incorporated  with  an  effective  current  density  J  .  Therefore  one  may 

♦♦eff 

define  an  effective  dielectric  susceptibility  x  which  takes  account  of 
the  effect  of  magnetization  M  as 


^  (jeff  .  g) 


(4.  3a) 


An  effective  polarization  then  is  defined  as: 


^eff  ^  jeff  .  g 


(4.3b) 


♦♦eff 

Similarly,  an  effective  dielectric  tensor  f  is  then  defined  as; 


♦♦eff  .  ♦♦eff 

e  =  1  +  4'irx 


(4.  3c) 


Substituting  Eq.  (4.3a)  into  Eq.  (4.2a),  one  obtains 


(7^ --^-^)  5(?.t)  =  +  (4.4) 

c  3t  c  at 


where  use  has  been  made  of  7  •  5  =  -5*  7lne.  Notice  that  Eq.  (4.4) 
can  also  be  obtained  from  Eq.  (4. 1)  directly  by  setting  M  =  0  (or  ^  =  1) 
and  P  =  P  .  This  means  that  whatever  effect  the  magnetization  has  is 
included  in  the  effective  dielectric  susceptibility  x  or  tensor  c 
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A  more  detailed  account  has  been  given  by  Pershan. 

If  we  consider  2(r,  t)  in  Eq.  (4.4)  as  a  superposition  of  pk  ne  waves 
of  wavelength  \/n,  and  the  elastic  or  spin  waves  are  of  wavelength  A, 

the  term  7  (S  •  In  r)  is  of  the  order  6(r  times  the  term  7^  Under  usual 

A 

experimental  conditions,  6c  (the  change  of  ordinary  dielectric  constant) 
is  much  smaller  than  unity,  and  ^  at  most  of  the  order  of  unity.  This 
term  can  usually  be  safely  neglected.  Then  we  have 


(4.5) 


c  at 


It  is  important  to  recognize  that  the  medium  contains  as  many  as 

9 

10  atoms  in  a  region  as  small  as  a  cube  of  the  light  wavelength  and 
therefore  can  be  regarded  essentially  as  a  continuum.  A  light  wave 
passing  through  such  a  medium  produces  an  effective  polarization 
at  the  point  r'.  This  polarization,  in  turn,  radiates  or  scatters 
electromagnetic  energy  in  all  directions.  Suppose  we  are  interested 
only  in  the  radiation  in  the  far  field  (Fraunhofer  region).  From  Eq.  (4.  5), 
the  electric  field  scattered  (or  diffracted)  at  point  r  by  the  effective 
polarization  within  a  volume  (d^r')  «  (X/n)^  is:^^’ 

.Ad,-.,  1  .3  ,1 

dE  (r.t)  '-jL-' - ' - ' - - »  --'J  k-r'|n, 

c  Ir-  r'l  t'=t - - - 

(4.6) 

where  r,  r'  and  the  unit  vectors  it-*  -♦,|  and  i  are  shown  in  Fig.  4.  1. 

The  quantity  t'  (=  t - - - )  is  the  retarded  time,  c  is  the  light  velocity 

in  free  space  given  before  and  n^  is  the  index  of  refraction  pertaining  to 
the  diffracted  electric  field  in  the  medium.  For  simplicity  we  have  taktm 
the  field  point  r  to  be  within  the  medium.  If  it  is  outside  the  medium, 
one  may  find  the  field  there  by  using  the  laws  of  refraction  at  the 
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illuminated  volume 
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boundaries.  However,  one  must  be  aware  of  the  fact  that  sometimes 

there  may  be  additional  contributions  to  the  scattered  electric  field  dE*^ 

62 

from  the  multiple  reflection  of  the  undiffracted  light  at  the  boundaries. 


The  electric  field  of  a  monochromatic  incident  light  of  frequency 
and  wavevector  within  the  medium  can  be  written  as; 


L 

L 

t) 

(4.7a) 

where 

2iTn.  n-co. 

'‘l 

_  i  _  i  L* 

X  c 

(4.7b) 

and  n.  is  the  index  of  refraction  pertaining  to  the  incident  electric  field 

inside  the  medium.  In  general,  when  the  polarization  of  the  incident 

light  is  different  from  that  of  the  diffracted  light,  n.  ^  n^.  The  quantity 

in  Eq.  (4.  6)  is  to  be  found  by  substituting  the  total  (incident  and 

diffracted)  electric  field  E  into  Eq.  (4.  3b).  However,  since  x  is 

modulated  slightly  by  coherent  elastic  waves  or  spin  waves  and  the 

interaction  between  light  and  coherent  elastic  waves  and  spin  waves  is 

very  weak,  one  may  approximate  S  in  Eq.  (4.  3b)  with  the  incident  electric 

field  in  Eq.  (4.  7a)  only  (the  "Born  approximation,  "  well  known  in 
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scattering  theory  ).  From  Eqs.  (4.  3b)  and  (4.  7a),  one  obtains  the 
polarization  at  point  r'  as: 


(4.8) 


where  )  i®  the  time  average  part,  and  6x  the  modulated  part  of 
the  effective  susceptibility.  Since  6x  (r’,t')  is  usually  modulated  by  a 
sinusoidal  elastic  (or  spin)  wave,  it  can  be  decomposed  into  the  Fourier 
spectra  of  the  elastic  (or  spin)  wave: 
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(4.9) 


where  |k|  =  ^  and  uj(k)  are,  respectively,  the  wavelength  and  frequencies 
of  the  coherent  elastic  or  spin  modes;  the  quantity  is  the  trans¬ 

verse  profile  of  the  mode  1<;  the  ±  sign  is  placed  in  front  of  a)(k)  to 
account  for  the  degeneracy  in  the  dispersion  relation  for  positive  and 
negative  angular  frequencies.  Equation  (4.  8)  and  Eq.  (4.9)  are  to  be 
substituted  in  Eq.  (4.  6)  to  evaluate  the  diffracted  electric  field  A 

furth'r  simplification  may  be  made  when  one  realizes  that  the  frequencies 

of  microwave  elastic  waves  (or  spin  waves)  are  very  small  (<  10^ ^Hz) 

1 4 

compared  to  the  light  frequency  ( ~  10  Hz).  Therefore  one  may 
2'*eff  2 

approximate  a  P  /at  in  Eq.  (4.  5)  as; 


a2Aeff  2 

LP_  -  -a)^?(^,t') 


(4.10) 


On  substituting  Eqs.  (4.  8),  (4.9)  and  (4.10)  into  Eq.  (4.6),  and 
integrating  over  the  illuminated  volume  V,  one  obtains: 

2  j[w,  ±co(k)t-lcl  - r] 

e‘*(?,  t)  =  -  4  V  {  [v  ^ - ) ' 


^  Xj^(r')  E(r')  e  d  r'j- 


(4.11) 


where  the  term  arising  from  <x)  is  dropped  and  =  E^E(r');  E(r') 
accounts  for  the  profile  of  the  light  beam.  We  have  used  the  fact  that 
if  lr|  »  Ir*!,  then 
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c 


(co^^  ±  co(k)) 


d  L  I-*  -♦. 


r-r'  ^ 


0,0.'.  .  .  . 

d  Li  ■*  /■*  ■*  I  \ 

- =-  I  •  (r -r') 

c  r 


r- r'l  s:  |r|  in  the  denominator  of  Eq,  (4,  6) 


^0ff  1  ^  cff 

Using  the  fact  that  6c  "  Trr  from  Eq.  (4.3c),  one  rewrites 


Eq,  (4. 1 1 )  as; 


2‘*(r.  t)  =  V  Pr 


where 


3  = 


°  =  S 


o 

M  -4.  -♦ 


-j(k^-k^+k).r'  ^  ^ 


E(r')Xj(r')  d'^r' 


(4.  12b) 


(4,  12c) 


The  vector  product  (i^xi^x)  merely  ensures  that  the  diffracted  light  is 

polarized  perpendicular  to  the  vector  the  vector  determines  the 

scattering  light  amplitude;  the  diffraction  integral,  D,  the  diffraction 
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pattern.  The  diffraction  pattern  has  been  discussed  by  McMahon, , 
Gordon”^  and  others.  In  our  case,  assuming  both  E{r')  and  X]^  (i  ' )  arc 
constant  throughout  a  volume  much  larger  than  optical  wavelength,  we 
immediately  obtain  the  phase  matching  (momentum  conservation) 
conditions; 

±  ^  (wj^=co^±w)  (4,13) 


However,  assuming  a  small  illuminated  rectangular  volume  a  x  w'  x  c 
(the  interaction  geometry  is  shown  in  Fig.  4.  2),  and  that  the  incident 
light  is  almost  parallel  to  the  wavefronts  of  the  elastic  or  spin  waves 
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(9g  —  0  in  Fig.  4.  2),  the  diffraction  pattern  5  at  the  observation  point  r 
can  be  written  as; 


D  =  a  w'c  Dj  DgDg 


(4.14) 


where  the  observation  point  r  is  at  a  small  angle  a  from  the  y-axis  and 
]8  from  y-z  plane,  and  use  has  been  made  of; 


L 


1! 


V  ‘'l  \  ''l  Sb 

1  k'  sin/3  +  t  k'  cosi3c6sa+^  k'  cosjjsina 

X  Lj  Lj  Z  Lj 


(4. 15a) 

(4. 15b) 

(4. 15c) 


The  quantities  D^,  and  Dg  are  given  by; 

sin  dk^asin^) 

^1  ”  ^  k^  a  sin)3 

sin[|(kl  cos/3cosa-k^  cos9„)w'] 

T"'  -  Lj  Lj 

^2  "  t(k^ cosjScosa  -k^ cos9g)w' 

sin[Hkl  cos)3sina  +  k,  sin9„-k)c] 

3  "  ^(k^  cos/3 sino  +  k^sin9g-  k)c 

As  expected,  the  Fraunhofer  diffracted  field  pattern  in  Eq.'  (4. 14) 
resembles  the  Fourier  transform  of  the  illuminated  rectangular  volume. 
In  Eqs.  (4.  16),  the  familiar  term  ^  is  maximum  when  E=  0.  One 
then  obtains: 

k'^Of  sin/3  =  0 

k|xos/8cosa  -  k^  cos9g  =  0  (4.17) 

k'  cos/Ssino  +  k.  sin 


(4. 16a) 

(4. 16b) 

(4.  16c) 


k 
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f»-w'  -H 


Fig.  4.  2  Geometry  of  light  diffraction. 


Fig.  4. 3 


Bragg  condition  for  k  =  k.  +  k. 
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or 


^  =  0 


cosa  - 
sina  +  k. 


cosOg 

sinSg 


0 

k 


(4.18) 


Equations  (4. 18)  can  be  represented  geometrically  in  Fig.  4.  3  and  can 
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easily  be  solved  as; 

-"72  <"r"d>] 

1  \ 

\  r,  /  2  2,1 

2ffn.  ^  ^^^d  2ff 

where  k  =  — r — ,  k'  =  — r —  and  k  =  have  been  used.  Notice  that 
L  ^  L  \  A 

we  have  neglected  the  small  difference  in  length  between  k^  and  k^.  In 
the  case  when  the  usual  Bragg  condition,  Eq.  (4, 19)  reduces 

to: 

sin  a  =  sin  0g  =  (4.  20) 


(4.19a) 

(4.19b) 


A-W  were  able  to  find  the  ratio  of  the  intensity  of  the  diffracted  light 
and  the  incident  light  by  integrating  Eq.  (4. 12a)  and  (4. 16)  approxi¬ 
mately.  Under  the  conditions  of: 

(1)  n.  =  n,  =  n 

1  d 

(2)  small  Bragg  angle  =  e„  s:  0 

(3)  a  rectangular  illuminated  volume  a  x  w'  x  c, 


they  obtained: 


^d  _  /ffw'  -  effV 

r  =  \^\  ^^ij  ; 

o 


(4.21) 
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where  is  the  appropriate  component  of  6?^^^  determined  by  the 

^  J 

vector  ^  in  Eq.  (4. 12b)  and  w'  is  tlie  illuminated  width  of  elastic  (or  spin) 
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wave  in  the  y-z  plane  in  Fig.  4.  2.  As  shown  by  McMahon,  one  may 
always  replace  w'  by  an  effective  w,  so  that  the  effect  of  arbitrary 
cross  sections  of  light  and/or  elastic  waves  (or  spin  waves)  can  be 
accounted  for. 

In  the  following  sections,  we  will  concentrate  on  the  quantities  6c- 

and  ^  in  Eq.  (4. 12b)  to  study  the  characteristics  of  light  scattering  from 

elastic  and  spin  waves.  The  intensity  calculation  can  then  be  calculated 

from  Eq.  (4.  21).  Since  it  is  understood  that  the  effective  dielectric 
•*eff 

tensor  c  is  the  one  we  are  concerned  with,  the  superscript  "eff"  is 
dropped. 


4.  3  Characteristics  of  Light  Diffraction  from  Coherent  Elastic  Waves 

The  wave  propagation,  at  optical  frequencies,  inside  a  crystal  is 

best  described  by  the  well-known  indicatrix,  which  is  an  ellipsoid  whose 

coefficients  are  the  components  of  the  relative  dielectric  impermeability 
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tensor  B..,  namely 
^  J 


where 


3 

I  B..  x.x. 

i.j=l  ^  J 


B..  =  B.. 

ij 


1 


(r^). 


ij 


(4.22a) 


(4.22b) 


andc  is  the  effective  dielectric  tensor.  Thus,  in  general,  a  small 

change  of  refraction  index  or  dielectric  impermeability  6B..  produced  by 

strains  in  the  elastic  waves  (photoelastic  interaction)  can  be 

described  in  terms  of  a  change  in  the  shape,  size  and  orientation  cf  the 

indicatrix.  The  change  is  conveniently  written  with  photoelastic 

constant  P..,  „  as: 
ijki 


-96- 


where  6B..  =  6B..,  and  the  quantities  6B..,  P..,  S,  -  are  all  dimension- 
ij  Ji  ^  ij  i]ki'  ki 

less.  Since  both  6B..  and  S,  .  are  symmetrical  tensors,  the  Eq.  (4.23) 

IJ  ki! 

can  be  contracted  into  the  matrix  form; 

6 


6B  =  I  P  S 
m  ,  mn  n 

n=l 


(m  =  1.2.  ...6)  (4.24) 


where  the  indices  are  contracted  according  to  the  Voigt  notation; 

11  -*  1  23  -♦  4 

22  -»  2  13  -♦  5 

33  -»  3  12  6 


P  =  P 
mn  ijki 


(m,  n  =  1, 2,  3  •  •  •  6)  (4.25a) 

(i.j,k,je  =  1,2,  3. ..) 


and 


6Bij 

‘^12 

6Bj 

'^6 

6B5 

‘^12 

'^22 

‘^23 

-* 

6B5 

6B2 

6B. 

4 

.‘^13 

*^32 

‘B33. 

'^5 

6B4 

6B3^ 

■s,. 

S2I 

S3l' 

'  Si 

iSe 

Si  2 

®22 

^23 

-♦ 

S2 

ic 

2^4 

_Sl3 

®23 

®33 

■ 

iSg 

ic 

2^4 

S3 

(4.25b) 


(4.25c) 


For  a  cubic  crystal  (class  m3m),  like  YIG,  the  photoelastic  matrix  P 


mn 


has  the  following  form  v/hen  referred  to  the  axes  along  the  cubic  edges 
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(principal  axes); 


6 


O 


u 


Pll 

Pi  2 

Pi  2 

0 

0 

0 

Pi  2 

Pll 

Pi  2 

0 

0 

0 

Pi  2 

Pi  2 

Pll 

0 

0 

0 

0 

0 

0 

P44 

0 

0 

0 

0 

0 

0 

P44 

0 

o 

0 

0 

0 

0 

P44 

(4.26) 


Since  we  are  interested  in  the  small  change  of  dielectric  tensor, 
6c.j,  we  want  to  express  6cjj  in  terms  of  From  Eq.  (4,  22b),  6?.^ 

can  be  written  as: 


6B. 


(i.j  =  1.2,  3) 


(4. 27) 


For  a  cubic  crystal,  the  unperturbed  dielectric  tensor  is  a  scalar  or 


where  6..  is  the  kronucker  delta  function, 
ij 

Combining  Eqs.  (4.24),  (4.26),  (4.  27)  and  (4.  28),  the  photoelastic 
interaction,  when  referred  to  the  principal  axes  of  a  cubic  crystal 
(class  m3m),  can  be  described  as; 
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(Sj)] 

i 

(Sj) : 

1 

(SJ  ! 

!  (4.29) 

(S4) ; 

(Sj) 

(Sg) 


In  terms  of  Eqs.  (4. 12a),  (4, 12b)  and  (4.  29),  assuming  the  Bragg 
angle  0^  e:  0  (Fig.  4.  2),  the  interaction  between  an  elastic  wave 
propagation  along  a  principal  axis  and  an  incident  light  polarized  also 
along  a  principal  axis  can  be  easily  determined  in  the  following  manner 
(Fig.  4.2). 

Since  0^  =  o  s:  0,  and  if  we  assume  that  the  y-z  plane  of  incidence  is 

horizontal,  the  incident  and  diffracted  ,  electric  fields  can  be 

o  o± 

written  as; 

s  1  E^  +  1  eJj  (4.30a) 

o  X  V  z  H  ' 


1  E^  +  1  E^ 

X  v±  z  H± 


(4.  30b) 


where  the  (±)  sign  indicates  whether  the  frequency  of  the  diffracted  light 
is  upshifted  (+)  or  downshifted  (-)  from  the  frequency  of  the  incident 
light  as  determined  by  the  phase  matching  conditions  (4. 13).  From 
Eqs.  (4.12a),  (4. 1 2b)  and  (4.  30),  we  obtain 


(4.31a) 
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E^ 

V 

v± 

=  c 

^*xz 

^H± 

^*zx 

L 

^^zz 

1 

i- 

(4.  31b) 


where  c  is  to  be  determined  by  Eq.  (4.  12a).  For  our  purpose,  the 
explicit  form  of  c  is  not  important  because  we  want  to  find  out  6f..  in 
Eq.  (4.  21)  only.  Note  that  Eq.  (4.  31b)  is  valid  even  in  a  coordinate 
system  which  is  oriented  in  an  arbitrary  crystallographic  direction,  as 
long  as  =  a  s  0.  Here,  since  x,  y,  z  axes  are  assumed  to  be  the 
principal  axes,  Eq.  (4.  31)  and  Eq.  (4.  29)  reveal  that; 

(a)  when  longitudinal  elastic  waves  propagate  along  a  principal 
axis  (thereby  producing  in  Fig.  4.  2),  the  incident  light 
and  diffracted  light  are  polarized  along  the  same  principal 
axes; 

(b)  when  shear  elastic  waves  propagate  along  a  principal  axis 
(thereby  producing  either  or  Sgg  in  Fig.  4,2),  the 
incident  and  diffracted  lights  are  polarized  along  two 
perpendicular  principal  axes.  Besides,  the  shear  strain 
$22  does  not  scatter  any  light  at  all. 

Notice  that  as  6c  ••  is  real  for  the  photoelastic  interaction,  the  diffracted 
^3 

light  intensity  is  independent  of  the  sign  of  the  frequency  shift.  As  we 
shall  see,  this  is  not  the  case  for  the  coherent  spin  wave  scattering. 

In  case  elastic  waves  do  not  propagate  along  a  principal  axis  or  the 
light  is  not  polarized  along  a  principal  axis,  one  merely  uses  the  method 
of  superposition  to  find  the  nature  of  the  interaction.  Two  examples 
are  used  to  illustrate  this. 
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Examplc  1 

The  interaction  geometry  of  the  photoelastic  interaction  is  shown  in 
Fig.  4.  4.  Longitudinal  (with  the  strain  S  )  or  shear  elastic  waves 


zz 

(with  the  strains  and  excited  by  a  transducer.  Again  notice 

from  Eq.  (4.11)  or  Eq.  (4.  31b)  that  the  strain  does  not  scatter  any 
light. 


'B 


[hO] 


Fig.  4.  4  Geometry  of  the  photoelastic  interaction. 


Since  the  plane  of  incidence  is  not  a  (100)  plane  (therefore  the  light  is,  not, 
in  general,  polarized  along  a  principal  axis),  one  may  either  work  out 
the  photoelastic  interaction  in  the  principal  axis  coordinate  system  by 
transforming  all  the  quantities  expressed  in  the  present  (x,  y,  z) 
coordinate  into  the  principal  axis  coordinate  or  transform  the  photo¬ 
elastic  constants  into  the  present  coordinate  system  and  use 

(4.  31b).  In  the  present  case  we  do  the  latter.  The  transformation 
from  the  principal  coordinate  system  ("cubic-edge"  coordinate  system) 
is  discussed  in  Chapter  II  in  terms  of  the  Euler  Angles  (<p,$,f(j).  From 
Kq.  (4.  29)  and'Eq.  (2.  4)  with  0  =  0,  =  0,  butp  ^  0,  we  obtain; 
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independent  of<p 

=  0 


(4.  32) 


We  therefore  conclude  that  for  a  small  Bragg  angle,  the  photoelastic 
interaction  here  is  independent  of  the  rotation  around  the  {10'^]  axis  in 
YIG.  In  this  case,  the  relation  between  the  polarizations  of  the  incident 
and  the  diffracted  light  can  also  be  readily  obtained  from  Eqs.  (4.  32)  and 
(4.  31b) 

Assuming  that  the  polarization  of  the  incident  light  is  at  an  angle  Q 
from  the  x-axis  (Fig,  4.4),  then,- 

(a)  for  shear  elastic  waves  (the  strain  S  ),  the  polarization  of 

xz 

the  diffracted  light  is  at  an  angle  (^  -  8)  from  the  x-axis; 

(b)  for  longitudinal  elastic  waves  (the  strain  S  ),  the  polarization 

z  z 

of  the  diffracted  light  is  at  an  angle  tan"^((Pj ^ /P^ 2)  tan  Q)  from 
the  x-axis. 

Example  2 

The  geometry  is  shown  in  Fig.  4.  5.  Again  when  we  use  Eq.  (2.2  ) 
withp  =0,  0,  but  Q  0,  we  obtain  for  the  strain 


p 

xxzz 

=  |sin^2ePjj  +  (1  -  isin^2e)Pj2  "  sin^26P^^ 

p 

zzzz 

=  (1  -  i  sin^2e)  Pj  1  +  i  sin^20  Pj  2  sin^2e  P^^ 

(4.33) 

=  p 

xzzz 

=  sin  40  (-  Pj  j  +  Pj  2  2P^^) 

and  for  the  shear  strain 
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Fig.  4.  5  Geometry  of  the  photoelastic  interaction  and  orientation 
of  the  long  axis  (z-axis)  of  the  sample. 
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=  isin4e(-Pj,  + 

Pxxiiz  =  j=i“^"’ll-Pl2-2P44> 

2  2  2 

P  =P  =  4sin  20P, ,  -  T  sin  20  P,  „  +  cos  20  P . . 
zxxz  xzxz  2  ''11^  ®12  ^  AA 


m 

In  particular,  for  a  longitudinal  strain  and  0  ='g: 


“xx 

P11  +  3P12-2P44 

■^11  ■^^12  "^^^44 

W 

1 

II 

(4.  35) 

^*zx 

;Pi1  +  Pi2''2P44 

^^11  ■^^12  "^^^44 

Equation  (4.  35)  is  used  in  Chapter  V'. 

From  these  two  examples,  one  realizes  that  the  intensity  of  the 
diffracted  light  can  be  calculated  by  using  Eq.  (4.21); 

6«2 


^  (^Op, 


(4.  36) 


v'^D 
P  o 


3  2 

where  P^  =  S  ;  the  quantity  P^^^  is  the  appropriate  photoelastic 

constant,  Vp  the  velocity  of  elastic  wave  (longitudinal  or  shear), 

P^  acoustic  power  density,  the  density  of  the  medium  and  S  is  the 
appropriate  strain. 

For  a  circular  acoustic  beam  of  an  effective  diameter  w,  Eq.  (4.  36) 
becomes; 


(4.37) 


where  P^  =  ^irw' 


P^  and  is  the  total  acoustic  power. 


For  YIg'^^’ 


at  \  =  1150  nm; 
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yi  ^  1. 67  X  10‘®  P^(mW)  (4.  38) 

o 

for  the  longitudinal  elastic  waves  propagating  along  a  {100}  direction  and 
off  *12- 

I .  f. 

~  ^  4.0x10  P^(mW)  (4.  39) 


for  the  shear  elastic  waves  propagating  along  a  {100}  direction  and 

P  1’ 

'  oil  '■  44- 
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4.  4  Characteristics  of  Light  Diffrac.ion  from 
Coherent  Spin  Waves 

5  6 

In  previous  work,  Auld  and  Wilson^*^  (A-W),  following  the  derivatior 

59 

of  Landau  and  Lifshitz,  calculated  the  intensity  of  light  diffracted  from 
coherent  spin  waves,  using  a  pertur  oed  dielectric  tensor  restricted  to 
terms  linear  in  the  total  crystal  mangetization.  For  many  magneto¬ 
optic  experiments,  this  is  sufficient  to  interpret  the  observations. 
However,  when  the  dielectric  tensor  is  modulated  by  a  coherent  spin 
wave,  one  must  look  for  all  terms  that  are  linear  in  the  small  signal  spin 
wave  amplitude,  not  just  terms  linear  in  the  total  magnetization  M.  For 

example,  when  the  magnetization  is  along  the  plus  z-axis,  terms  like 
2 

m  M  and  m  M  ,  although  of  higher  order,  are  linear  in  the  small 
X  z  y  ^ 

signal  spin  wave  amplitudes  m^  and  m^.  Our  derivation  of  the  perturbed 

dielectric  tensor,  detailed  in  Appendix  III,  is  straightforward  and 

essentially  follows  that  of  Landau  and  Lifshitz,  with  the  exception  that  ^ 

linear  terms  in  m  and  m  are  retained.  The  result  for  a  cubic  crystal 
X  y 

(class  m3m),  like  YIG,  is  of  the  form,  when  expressed  in  the  "cubic 
edge"  coordinate  system; 


6<..iX)  =  6fI.(X)  +  j6c'.'(X) 
ij  ij  ■’  ij 


(i,  j  =  1,2,  3) 


(4.  40a) 


XX 

01  iW 

0 

0 

0 

1 - 

CM  X 

6c' 

yy 

gl2(^) 

gll(X) 

gj  2^^  ^ 

0 

0 

0 

2 

Ot 

y 

<z 

gii(x) 

gll(X) 

0 

0 

0 

2 

^‘z 

6c' 

yz 

0 

0 

0 

844 
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0 

“y^'z 

'4x 

c 

0 

C 

0 

g44(M  0 

a  a 

z  x 

^^xy 

0 

0 

0 
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0 

844^^) 
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X  y 

(4.  40b) 
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where  6c'.  -  6c'-. 

ij  Ji 


1 1 

XX 

0 

0 

0 

1 1 

6c 

yy 

0 

0 

0 

1 1 

6c 

zz 

0 

0 

0 

1 1 

6c 

yz 

f(x) 

0 

0 

1 1 

6c 

zx 

0 

f(x) 

0 

1 1 

6c 

xy 

0 

0 

f(x) 

(4.  40c) 


where  6c-  •  =  -6c  -- 
ij  Ji 

For  simplicity  we  have  assumed  that  the  medium  is  lossless  and 

therefore  6c  (X)  is  Hermitian.  The  quantity  f(>)  is  the  linear  coefficient, 

g  (X)  the  quadratic  coefficients,  a  (=M  /M),  a  (=M  /M),  and 
mn  XX  y  y 

a^(-  M^/M)  the  direction  cosines  of  the  magnetization;  both  f(X)  and 

gmn^^)  ^  function  of  optical  frequencies  or  wavelength  X.  When 

g  -  0,  the  tensor  6c  reduces  to  that  of  A-W.  Note  that  g  is 

^mn  '  ^mn 

identical  in  form  to  in  Eq.  (4.29). 

The  calculation  of  the  intensity  of  light  diffracted  from  spin  waves 

can  be  made  using  Eqs.  (4.21),  (4.  31b)  and  (4.  40).  Notice  that  in 

Eq.  (4.  40)  the  perturbed  dielectric  tensor  is  now  modulated  by  z- directed 

spin  waves  (wavevector  ^  parallel  to  the  dc  magnetization  vector,  which 

is  along  the  z-axis),  with  frequency  u).  Furthermore,  for  simplicity, 

let  us  assume  that  the  z-directed  spin  waves  are  circularly  polarized  so 

that  a  anda  can  be  written  from  Eq.  (2. 1 6)  with  (uJ.  =co)  as; 

X  y  K 

a. 


M 


X  m  /  ^  1  V 

-X  '  IT  = 


(4.  41a) 
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M„  ^ 

a  =  ^  =  -^sinCuft-kz)  (4.41b) 

Here  the  x,  y,  z  axes  do  not  necessarily  coincide  with  the  principal  axes. 

We  are  now  ready  to  do  the  calculation.  Just  as  in  the  last  section, 
we  use  two  examples  to  illustrate  this. 

Example  1 

The  geometry  of  light  diffraction,  identical  to  Fig.  4.4,  is  shown  in 
Fig.  4.6;  the  spin  wave  amplitudes  are  given  in  Eq.  (4.41). 
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Fig.  4.  6  The  geometry  of  light  diffraction  from  spin  waves; 
the  internal  field  profile  Hj  of  the  bar  is  also 
indicated. 


Since  the  plane  of  incidence  is  not  a  {100}  plane,  we  obtain  from 
Eq.  (4.  40)  and  Eq.  (2.2),  with  iji  =  6  =  0,  but  (p  0: 


-loa- 


independent  of( 


(4. 42)' 


6f  =  6c  =  fi .  >  a  a  +  i  f  a 

*xz  *zx  '’44  2.  X  •’  ^ 


wiu'iH?  only  tonus  linear  in  a  ,  a  are  kept,  and  a  —  1.  Substituting 

X  y  z 

Eq.  (4.  41 )  into  Eq.  (4.42),  we  obtain; 

^^zx  "  [®44“z  M  J  ^  M  ■  *^2)] 


(4.  43a) 


JILf/p  ft  ^  .  np-j(u)t-kz)] 

2M  L'®44“z^‘^  ®  +tg44“z  ^  J 


6c  =  6c  ' 
xz  zx 


(4.  43b) 


i[<i44“z  -  f)  -  <«44“z  +  f)  e- 


Substituting  Eqs.  (4.  43)  into  Eq.  (4.  31b),  and  noting  the  phase  matching 
condition  in  Eq.  (4.  13),  we  obtain; 


§44^2  +‘ 


-  P-JH 

^  2M 


(4.  44) 


g44“z±f 


J  L 


where  a  subscript  v  or  H  denotes  whether  the  polarization  of  the  electric 
field  vectors  is  perpendicular  (x-axis)  or  parallel  to  the  horizontal  (y-z) 
plane  of  incidence.  The  plus  or  minus  sign  indicates  whether  the 
frequency  shift  is  upward  (anti-Stokes  line  w^+u))  or  downward  (Stokes 


The  asterisk  on  6c  '  denotes  complex  conjugation. 

X  z 
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line  -  OJ). 

From  Eq.  (4.  44)  one  can  deduce  that  the  maximum  diffracted  light 
intensity  depends  on  the  incident  light  polarization  and  the  sign  of  the 
frequency  shift.  If  f  and  g^^  have  the  same  sign,  will  give  the 
maximum  diffracted  light  intensity  for  the  upshifted  case,  whereas  E^ 
will  for  the  case  of  downshifted  diffracted  light.  Furthermore, 
reversing  the  direction  of  the  magnetic  field  changes  the  sign  of  and 
the  sense  of  the  spin  wave  circular  polarization.  Consequently,  the 
above  situation  remains  unchanged,  because  the  ratio  (g^.a  0.  ./fCK„)  .">till 
has  the  same  sign. 

The  "anomalous"  experimental  observations  that  Collins  and 
57 

Wilson  have  reported  fit  precisely  this  description,  provided  one 
remembers  that  low  to  medium  wavenumber  magnetostatic  waves  are 
backward  waves  with  the  wavevector  direction  opposing  the  direction  of 
the  power  flow.  They  found  that,  with  a  laser  beam  polarized 
perpendicular  to  the  plane  of  incidence,  the  scattered  beam  was  always 
strongest  on  the  left  side  of  the  beam  when  using  the  left-hand  antenna, 
but  was  always  found  on  the  right  side  of  the  main  beam  when  using  the 
right-hand  antenna.  These  two  cases  correspond  to  frequency  upshifted 
diffracted  light,  as  indicated  in  Fig.  4.  7.  The  converse  effects, 
i.  e. ,  frequency  downshifted  diffracted  light  when  the  input  laser  polar¬ 
ization  was  parallel  to  the  plane  of  incidence,  are  also  similarly 
explained. 

Other  experimental  evidence  is  furnished  by  Smith's  work^^  in 
which  he  found  similar  situations  for  light  scattering  from  shear 
magnetoelastic  waves  with  a  large  mixture  of  spin  waves.  Again  he 
reported  it  as  "anomalous.  "  Section  5.2  gives  a  more  detailed  account. 
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Fig.  4.  7  Light  diffraction  from  coherent  magnetostatic  waves. 
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Example  2 

The  geometry  of  light  diffraction  is  identical  to  Fig.  4.  5  where  the 
applied  magnetic  field  is  along  the  z-direction.  Again  when  we  use 
Eq.  (2.  2)  with  (p  =  0  =  0,  but  0^0,  we  obtain; 


-gi2-2«44>“z“x 

=  isin4e(-gjj  +g,2  +2g^^)a^a^  (4.45) 

''c  2  2  2 

2egjj-isin  2egj2  +  cos  2Qg^^]a^a^  + 


tr 

In  particular,  for  6  ='0»  we  obtain; 


‘‘xx 

^^xz 

_  1 

- 

(gll-gl2+2g44>“z“x-^j''‘“y 

^^zx 

^^zz 

“  4 

(gll-gl242«44>“x“z  +  ^j‘'“y 

•^8“x“z 

(4.  46a) 


where 


Ag 


’ll  ^12 


gi  9  -  2g 


44 


(4.  46b) 


6  8  69 

Notice  that  when  Ag  =  0  (Cauchy  relation),  '  the  medium 
becomes  "isotropic,  "  i.  e. ,  x,  y,  z  axes  are  determined  by  the  dc  magnet¬ 
ization  (2:  oc  )  and  small  signal  spin  wave  amplitudes  a  and  a  only, 
z  X  y 

Substituting  Eqs.  (4.  41)  and  (4.46a)  into  Eq.  (4.  31b),  as  in  Example  1, 
we  obtain; 
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v± 


H± 


=  c 


m 

2M 
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I<girgi2-"2g4.j)a^  ±  f 


j(g],-gl2+2g44)az  +f 


-^g  O,  i 


'H 


(4.47) 
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Except  for  the  different  g's.  Eg,  (4.47)  is  similar  to  (4,44),  and 
therefore  most  of  the  characteristics  discussed  there  are  applicable  here, 
When  Ag  =  0,  Eq.  (4,  47)  is  identical  to  Eq.  (4,  41), 
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CHAPTER  V 

EXPERIMENTAL  RESULTS  OF  LIGHT  DIFFRACTION  FROM 
COHERENT  ELASTIC  AND/OR  SPIN  WAVES  IN  YIG 

In  this  chapter,  experimental  results  of  light  diffraction  from  coherent 

shear  elastic/spin  waves,  and  longitudinal  elastic/spin  waves  in  YIG  are 

presented.  Of  particular  interest  is  the  strong  infrared  Bragg  scattering 

7  8 

from  coherent  spin  waves  reported  in  section  5.  3.  The  experimental 
results  in  this  case  appear  to  be  in  good  agreement  with  the  theory  on 
light  scattering  from  coherent  spin  waves  derived  in  the  last  chapter. 

The  chapter  begins  with  a  review  of  previous  experimental  work  on 
light  diffraction  in  YIG,  followed  by  a  description  of  the  optical  experi¬ 
mental  setup.  In  section  5.  2,  results  of  Bragg  diffraction  from  shear 

73  74  75  76 

elastic /spin  waves  similar  to  the  work  of  Dixon  and  Smith  ’  ’  are 

described.  In  the  last  section  of  the  chapter,  5.  4,  we  report  on  measure¬ 
ments  of  temporal  longitudinal  elastic/spin  wave  conversion  efficiencies, 
which  are  difficult  with  conventional  microwave  techniques  whenever 
intrinsic  non-uniform  fields  due  to  demagnetizing  effects  occur  inside  a 
ferrimagnetic  cample. 

71 

Dixon  and  Matthews  first  employed  coherent  light  scattering  from 

shear  elastic  waves  to  measure  the  photoelastic  constants  of  YIG.  In  a 

72 

later  paper,  Dixon,  using  an  optical  heterodyne  system,  reported  large 

Faraday  rotations  of  shear  elastic  waves  and  a  substantial  frequency 

dispersion  of  shear  elastic/spin  waves  in  YIG.  Smith  reported  light 

scattering  from  shear  elastic/spin  waves  that  had  been  excited  both 
7475  76 

acoustically  '  and  magnetically.  Of  necessity,  he  employed  rf 
magnetic  and  electric  field  excitation  at  high  microwave  power  inputs  of 
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5  to  300  watts;  therefore,  nonlinear  interactions  probably  occurred.  He 

also  observed  the  dependence  of  the  scattered  light  intensity  on  the 

incident  light  polarization  and  diffracted  light  frequency  shifts.  Based 

56 

on  Auld  and  Wilson's  (A-W)  theory,  he  interpreted  this  dependence 

qualitatively  as  due  to  the  interference  of  the  scattered  light  from  elastic 

77 

and  spin  waves.  Subsequently,  A-W  observed  the  field  dependence  of 
scattered  light  intensity  for  longitudinal  elastic/spin  waves  in  an 
obliquely  magnetized  YIG  bar,  but  their  experimental  data  were  not 

C  C 

accurate  enough  to  either  prove  or  disprove  their  calculation.  Collins 
57 

and  Wilson  studied  the  light  scattered  from  magnetostatic  waves;  they 

reported  an  "anomalous  behavior"  unexplained  by  the  A-W  theory  which, 

however,  can  be  interpreted  quite  satisfactorily  by  our  new  theory 

79 

described  in  section  4.  4.  Most  recently,  Desormiere  et  al.  reported 
on  related  experiments  on  light  diffraction  from  magnetostatic  waves  and 
observed  "anomalous  behavior"  similar  to  what  Collins  and  Wilson  had 
reported. 

5.  1  Setup  for  Light  Diffraction  Experiments 

The  experimental  setup  shown  in  Fig. 5.  la  comprises  a  microwave 
system  which  is  used  for  injecting  and  monitoring  coherent  elastic  waves 
via  piezoelectric  transducers  in  YIG,  and  an  optical  system  for  sending 
a  coherent  laser  beam  and  detecting  the  diffracted  light  from  the  elastic 
and/or  spin  waves  in  the  sample. 

The  microwave  system  is  similar  to  that  in  Fig.  3. 1,  except  that 
a  TWT  (HP  489A,  maximum  power  output  =  2.0  watts  at  L  band),  capable 
of  amplification  of  fisec  pulses,  is  used  following  the  PIN  modulator  to 
boost  the  power  output  of  the  signal  generator.  When  using  bulk 
tr  ansducers,  which  arc  less  efficient,  have  narrower  bandwidth,  but  can 
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Fig.  5.  1(a)  Block  diagram  of  the  optical  experimental  setup. 
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handle  more  power  than  the  thin  film  transducer,  a  power  oscillator 
(a  modified  AIL  125,  externally  pulse  modulated), capable  of  20  to  50  watts 
power  output,  was  sometimes  used  to  facilitate  the  location  of  the 
diffracted  light  signal.  A  peak  power  meter  (Narda  6SA3A)  is  used  to 
measure  the  input  pulse  power  to  the  transducers.  The  microwave 
receiver  system  is  essentially  the  same  heterodyne  system  described  in 
s(;ction  3.  1,  but  with  a  limiter  and  an  attenuator  preceding  the  mixer/ 
preamplifier  to  prevent  an  excessive  power  loading  of  the  mixer.  The 
receiver  system  is  used  to  measure  the  transducer  losses,  optimize  the 
stub  tuner  (which  was  as  close  to  the  load  as  feasible),  monitor  the 
acoustic  echo  from  the  sample  and  provide  supplementary  information 
regarding  the  interpretation  of  the  optical  signal. 

The  samples  used  in  the  optical  experiments  are  rectangular  YIG 

bars  with  two  end  surfac  s  as  well  as  two  opposite  long  surfaces  optically 

pjlished.  The  orientation  of  the  crystals  is  described  separately  in  each 

section.  No  optical  antireflection  coating  is  used  on  the  long  polished 

surfaces.  The  piezoelectric  transducers  used  in  the  experiments  for 

injecting  elastic  waves  into  the  sample  are  either  bulk  transducers  or  a 

sputtered  ZnO  thin  film  transducer.  The  bulk  transducers  are  X-cut 

LiNbO^  plates  for  generation  of  linearly  polarized  shear  elastic  waves, 

or  Z-cut  Li- doped  ZnO  plates  for  generation  of  longitudinal  elastic 

waves.  The  dimensions  of  the  bulk  transducers  are  0.  005"  in  thickness 

and  about  0.  120"  in  diameter.  For  the  bulk  transducer  bonded  with 
80 

phenyl  benzoate,  typical  insertion  loss  (i.  e.  ,  the  power  ratios  of  the 
input  pulse  to  the  first  echo)  is  50  db  or  more  at  1  GHz.  The  sputtered 
thin  film  ZnO  transducer  (courtesy  of  Dr.  N.  F.  Foster  of  Bell 
Telephone  Labs.  ,  Inc.  ),  on  the  other  hand,  has  a  very  large  bandwidth 
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and  initially  exhibited  only  a  20  db  insertion  loss  but,  unfortunately, 
deteriorated  later  to  about  36  db. 

The  optical  system  consists  of  the  following  components  (Fig.  5.1a): 

(a)  A  laser:  Spectra -Physics  Model  125A  modified  to  give  maximum 
power  output  at  1150  nm;  power  output,  50  mW;  transverse 
mode,  TEMqqJ  dc  excitation.  ^ 

(b)  A  set  of  three  mirrors  mounted  on  micro-positioners  for 
coarse  and  vernier  movements  in  the  horizontal  and  vertical 
directions.  The  mirrors  are  dielectric-coated  and  give  99% 
reflectance  at  45°  incident  angle  for  "p-  and  s-polarized" 
light.  ^  One  of  the  three  mirrors  is  not  used  at  45°  angle  of 
incidence  and  probably  should  be  replaced  by  a  more  suitable 
mirror. 

(c)  Two  Glan-Thompson  polarizers  mounted  on  rotating  stages 
to  allow  selection  of  the  polarization  of  the  transmitted  light. 

The  polarizers,  with  an  extinction  ratio  of  10  :1,  prove  to  be 
ex*,remely  useful  in  improving  signal-to-noise  ratio  when  the 
diffracted  light  is  polarized  at  90°  with  respect  to  the  incident 
light,  as  in  the  cases  of  diffraction  from  shear  elastic  waves 
and  spin  waves. 


A  1  mW  1150  nm  laser  was  used  in  the  early  stages  of  the  work  and 
we  did  observe  diffracted  light  signals  with  10  watts  microwave 
power  into  a  LiNbOs  transducer  (~50  db  insertion  loss)  after  three 
hours  of  integration  time  with  a  PAR  160  boxcar  integrator. 

^  An  rf  excitation  gives  more  pov  er  output  and  stability  but  may  cause 
electromagnetic  interference  with  other  instruments  in  use.  There¬ 
fore  a  dc  excitation  sometimes  is  more  desirable. 

+  The  phrase  "p-polarized"  refers  to  light  polarized  in  the  plane  of 
incidence,  while  "s-polarized"  refers  to  light  polarized  perpendicular 
to  the  plane  of  incidence.  One  should  beware  that  some  manu¬ 
facturers  may,  without  warning,  supply  dielectric  mirrors  which 
give  99%  reflectance  at  45°  angle  of  incidence  for  only  s-polarized 
light. 
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(d)  Three  lenses  are  employed,  as  shown  in  Fig.  5.  la.  One  is  used 
for  focusing  tl  e  5  mm  laser  beam  onto  the  surfaces  of  the 
sample,  the  other  two  for  collecting  diffracted  light  and 
focusing  it  into  a  slit  in  front  of  the  photomultiplier. 

(e)  An  iris  in  front  of  the  first  collecting  lens  and  a  slit  in  front  of 
the  photomultiplier,  as  shown  in  Fig. 5.  la,  are  found  to  be  very 
nelpful  in  reducing  the  scattered  background  light. 

(f)  An  Amperex  150  CVP  photomultiplier  (S-1  photocathode, 
operating  at  1,  500  V)  is  cooled  to  the  dry  ice  temperature  in  a 
housing  (Products  for  Research  Model  TE-200).  Cooling  of 
the  photomultiplier  to  dry  ice  temperatures  is  necessary  here 
because  of  the  high  dark  current,  low  quantum  efficiency  of  the 
S-1  cathode,  and  the  low  level  diffracted  light  signal  obtained 
in  the  experiments.  A  heated  ring  is  installed  cn  the  outer 
surface  of  the  viewing  window  of  the  photomultiplier  housing  to 
prevent  the  formation  of  dew  on  the  cold  surface.  A  low  pass 
filter  (transmittance  80%  at  1150  nm  and  above,  from 

IR  Industries)  is  also  placed  in  front  of  the  viewing  window  to 
prevent  any  stray  room  light  from  getting  into  the  photo¬ 
multiplier. 

(g)  A  mica  half-wave  plate  for  1150  nm  is  also  used  between  the 
laser  and  the  polarizer  to  rotate  the  polarization  of  the  incident 
light.  A  Kodak  IR  Phosphor  Card,  which  converts  the  near 
infrared  to  visible  red  light,  is  useful  for  optical  alignment 

but  is  not  sensitive  enough  to  observe  the  diffracted  light. 

(h)  A  sample  stage  placed  in  the  airgap  of  the  pole  pieces  of  a 
Varian  6"  magnet  can  be  used  to  position  the  samples  properly 
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with  respect  to  the  incident  light  beam.  When  experimental 
data  are  taken  the  sample  is  not  moved  with  respect  to  the 
applied  magnetic  fields.  Other  details  of  the  optical  experi¬ 
mental  setup  are  shown  in  Fig.  5.  lb. 

The  electronic  system  following  the  output  of  the  photomultiplier 
contains  (^ig.  5.  la): 

(a)  Load  Resistor  R^.  Either  a  1  K  ohm  resistor  for  better 
signal-to-noise  ratio  or  a  50  ohm  resistor  for  faster  risetime 
of  the  diffracted  light  pulse  is  used. 

(b)  For  R^  =  1  K,  an  amplifier  with  10  MHz  bandwidth  and  an 
amplification  of  100  is  adequate  for  the  input  of  the  boxcar 
integrator.  However,  for  R^  =  50  ohm,  an  amplification  of 
1,  000  and  a  larger  bandwidth  (70  MHz)  are  needed  to  observe 
the  fast  diffracted  pulse. 

(c)  A  PAR  Model  160  boxcar  integrator  is  used  to  measure  the 
diffracted  light  signal  and  take  the  experimental  data.  The 
output  of  the  boxcar  integrator  is  fed  into  a  digital  voltmeter 
for  display  and  an  x-y  recorder  for  permanent  records.  The 
display  of  the  digital  voltmeter  is  found  to  be  very  helpful  for 
optimizing  the  signal  in  a  dark  room. 

The  experimental  difficulties  in  performing  optical  probing  of 
elastic  and/or  spin  waves  are  dominated  by  the  poor  signal-to-noise 
ratio,  worsened  in  particular  by  the  poor  quantum  efficiency  of  the 
S-1  photocathode  and  comparatively  low  power  laser  available  at 
1150  nm.  If  one  is  working  in  the  visible  spectrum,  a  factor  of  1, 000 
improvement  in  photocathode  quantum  efficiency  can  be  obtained  and,  in 
addition,  perhaps  a  factor  of  10  increase  in  laser  power  output  before 


*  I 
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(a)  Laser  beam  incident  system,  the 
laser  (not  shown)  is  at  the  left 
cf  the  picture. 


(b)  The  magnet  and  sample  holder 


(d)  Microwave  apparatus 


(c)  Optical  detection  eysttm 


Fig.  5.1(b)  Optical  experimental  setup. 
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problems  of  heating  the  sample  under  investigation  occur.  Wilson, 

82 

following  Anderson,  gave  a  simple  calculation  of  the  minimum 
detectable  intensity  of  an  optical  signal  diffracted  from  either  a  square 
wave  or  a  pulse  modulated  clastic  or  spin  wave.  There  is  nothing 
further  we  want  to  add,  except  to  point  out  that  under  our  experimental 
situation  (with  laser  power  50  mW  and  acoustic  power  inside  the  medium 
less  than  1  mW)  the  noise  is  always  dominated  by  the  background  noise, 
which  comes  either  from  scattering  from  imperfections  on  the  surface 
or  from  inside  the  crystal.  While  the  signal  intensity  is  proportional  to 
the  incident  laser  power,  the  background  noise  due  to  scattering  from 
the  crystal  imperfections  is  proportional  only  to  the  square  root  of  the 
laser  power.  Consequently  a  high  power  laser  is  still  helpful  even  in 
this  case.  Of  course  when  the  noise  is  dominated  by  the  dark  current, 
the  high  power  laser  is  even  more  helpful. 
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5.  2  Infrared  Bragg  Scattering  from  Shear 

Elastic/Spin  W  aveti 

Bragg  scattering  from  shear  elastic /spin  waves  in  YIG  has  been 
studied  extensively,  first  by  Dixon,  then  by  Smith.  Dixon 

has  exhausted  almost  all  the  experimental  cases  of  interest  with  this 
kind  of  light  scattering,  with  the  exception  of  the  dependence  of  the 
intensity  of  the  diffracted  light  on  the  incident  light  polarization  and  the 
the  sign  of  the  frequency  shift.  Smith  ’  °  first  reported  this  depend¬ 
ence  and,  based  on  the  A-W  theory,  interpreted  it  qualitatively  as  the 
interference  between  the  light  diffracted  from  elastic  waves  and  from 
spin  waves.  However,  he  found  an  "anomaly"  which  his  interpretation 
fails  to  explain.  The  "anomaly"  occurs  when  the  light  is  scattered  by 
the  shear  elastic/spin  waves  with  a  large  mixture  of  spin  v/aves.  He 
thought  this  "anomaly"  might  be  due  to  the  large  amplitude  of  the  shear 
elastic/spin  wave  signal  used  in  his  experiments.  In  fact,  he  did 

report  some  nonlinear  behavior  of  light  scattering  from  shear  elastic/ 

,  75 

spin  waves. 

In  this  section  we  report  some  results  on  light  scattering  from 
shear  elastic/spin  waves  under  essentially  the  same  experimental 
conditions  as  Smith’s,  but  employing  a  much  reduced  microwave  power 
input.  Within  the  range  of  microwave  power  which  we  operated,  no 
nonlinear  behavior  in  the  light  scattering  was  observed. 

Our  experiments  were  performed  using  a  rectangular  YIG  bar 
which  measures  4  mm  x  4.  6  mm  x  12.  7  mm  and  is  oriented  at  a  <100> 
axis  (Fig.  4.4  or  4.6).  A  I.lNbOj  shear  wave  piezoelectric  transducer 
with  a  50  db  insertion  loss  was  bonded  to  the  YIG  bar  and  used  to  excite 
shear  elastic  waves.  The  peak  power  of  a  1.0  ^ec,  1.3  GHz  micro- 
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wave  pulse  used  was  between  100  mW  and  500  m\\  (versus  Smith's 
5-300  W). 

Figure  5.2  shows  typical  diffracted  light  pulses  observed  with  the 
light  beam  near  the  free  end  of  the  bar.  A  diffracted  light  pulse  results 
each  time  the  shear  elastic  wave  packet,  which  is  reflected  back  and 
forth  in  the  sample,  crosses  the  path  of  the  lichf.  beam.  In  the  figure 
the  insert  is  a  trace  from  a  Tektronix  585  scope,  while  the  graph  was 
plotted  from  the  signal  conditioned  with  the  PAR  boxcar  integrator. 

The  quantities  and  ,  defined  in  Eq.  (4.  30),  denote  the  polariz- 
ations  of  the  incident  and  diffracted  light;  the  minus  sign  following  the 
subscript  v  indicates  the  downshifted  diffracted  light.  The  quantity 
is  the  external  applied  magnetic  field.  Figure  5.  3  is  a  plot  of  the 
diffracted  light  intensity  as  a  function  of  the  external  applied  magnetic 
field  Hq,  when  the  light  beam  is  incident  near  the  mid-point  of  the  long 
axis  of  the  sample  at  the  time  when  the  shear  elastic/spin  v/ave  packet 
crosses  the  path  of  the  incident  light  beam  for  the  first  time.  The  solid 
curve  was  obtained  when  the  incident  light  was  polarized  in  the  horizontal 
plane  of  incidence  (hence  E^)  and  the  downshifted  diffracted  light, 
polarized  perpendicular  to  the  plane  of  incidence  (hence  )•  On  the 
other  hand,  the  dashed  curve  was  obtained  when  the  incident  light  was 
polarized  vertically  (E^)  and  the  downshifted  diffracted  light,  polarized 
horizontally  (E^_).  The  Bragg  angle  for  the  shear  elastic  waves 
outside  the  sample  was  13.3°  (Eq.  (4.20)). 

The  sample  was  then  rotated  180°  around  an  axis  perpendicular  to 
the  horizontal  plane  of  incidence  so  that  the  upshifted  diffracted  light 
’'ould  be  obtained.  The  position  of  the  sample  was  readjusted  so  that 
the  light  beam  was  incident  near  the  mid-point  of  the  long  axis  again. 
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wave  pulse  propagating  along  a  <100>  direction  and 
polarized  along  a  <110>  direction  in  YIG.  The  graph 
was  recorded  with  the  signal  from  the  PAR  boxcar 
integrator.  The  insert  is  the  scope  trace. 
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A  plot  of  the  diffracted  light  intensity  as  a  function  of  the  external  applied 

fields  with  either  horizontally  or  vertically  polarized  incident  light  gave 

exactly  identical  curves  in  reverse  order.  In  contrast  to  the  previous 

case,  the  solid  curve  was  obtained  when  the  incident  light  was  vertically 

polarized  (E^)  and  the  upshifted  diffracted  light,  horizontally  polarized 

(Ejj^),  while  the  dashed  curve  was  obtained  when  the  incident  light  was 

horizontally  polarized  (Ejj)  and  the  upshifted  diffracted  light,  vertically 

polarized  (E^^).  These  results  are,  in  fact,  identical  to  those  Smith^^ 

has  obtained  with  high  microwave  power  input.  Therefore  we  conclude 

that  the  "anomaly"  he  reported  is  not  due  to  the  nonlinear  effect,  but  is 

intrinsic  with  the  light  diffraction  from  spin  waves. 

S  S 

Notice  in  Fig.  5.  3  that  when  H  >  H  (H  is  the  applied  field  value 
®  o  c  c 

for  which  the  shear  crossover  point^^'  is  at  the  mid-point  of  the  long 

axis),  the  diffracted  light  intensity  becomes  independent  of  the  incident 

light  polarization,  which  is  expected,  since  only  the  negative  rotating 

72  S 

shear  elastic  wave  is  present  here.  In  the  case  of  (and  the 

downshifted  diffracted  light),  E^  gives  a  systematically  larger  diffracted 

signal  E^_  than  the  E^.  (The  intensity  oscillations  in  Fig.  5.  3  result 

72 

from  the  acoustic  Faraday  rotation.  )  This  is  consistent  with  Smith's 

theory,  in  which  he  pointed  out  that  this  effect  is  due  to  the  interference 

of  the  light  diffracted  from  elastic  waves  and  spin  waves.  This  also 

check.r  with  our  new  theory,  derived  in  section  4.  4  (Example  1),  when 

our  calculations  (rather  than  A-\V's)  of  the  contribution  of  the  diffracted 

light  intensity  from  spin  waves  arc  used  in  Smith's  theory. 

When  11  i  H  (and  the  downshifted  diffracted  light),  the  shear 
o  c 

elastic/spin  waves  contain  a  large  mixture  of  .‘-.pin  waves.  In  this  case, 
e|j  gives  1  larger  diffracted  signal  (a  peak  marked  "spin"  in  Fig.  5.  3) 
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EXTERNAL  MAGNETIC  FIELD  (Oe) 

Fig.  5.  3  Diffracted  light  intensity  (in  arbitrary  units)  of  shear 


elastic/spin  waves  as  a  function  of  external  magnetic 
field;  the  light  beam  is  incident  near  the  mid-point 
of  t*  ,  long  axis.  The  insert  is  the  scope  trace  swept 
f:  right  to  left  when  s?  h®,  indicating  small 

group  velocity  dispersion.  A  full  explanation  is 
given  in  the  text. 
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which  Smith  was  not  able  to  explain  (the  "anomaly").  According  to  our 
theory  for  spin  wave  scattering,  detailed  in  section  4.  4,  Example  1 ,  the 
incident  electric  field  eJj  should  indeed  give  a  larger  diffracted  signal 
for  the  downshifted  diffracted  light. 

All  the  phenomena  of  the  converse  case,  i.  e.  ,  the  upshifted  case, 
can  be  similarly  explained. 

In  order  to  ascertain  that  the  peak  (marked  "spin")  in  Fig.  5.  3  does 
result  from  light  diffracted  from  coherent  spin  waves,  another  plot  of 
the  diffracted  light  intensity  as  a  function  of  external  magnetic  fields 
was  made  (Fig.  5.4).  In  obtaining  Fig.  5.4,  the  incident  light  was 
polarized  at  an  angle  of  ^  from  the  x-axis  (Fig.  4.  4  or  4.  6)  and  the 
diffracted  light  was  always  downshifted  in  frequency.  The  solid  curve 
was  obtained  when  the  polarizer  in  front  of  the  photomultiplier  was  set 

3 

to  transmit  only  the  diffracted  light  polarized  at  an  angle  of  from  the 
x-axis,  while  the  dashed  curve,  only  the  light  polarized  at  an  angle  of^ 
from  the  x-axis.  According  to  Example  1  in  section  4.  3,  the  diffracted 
light  from  the  elastic  waves  must  be  polarized  in  the  same  direction  as 
the  incident  light;  consequently  the  diffracted  signal  of  the  solid  curve 
could  only  result  from  the  light  diffracted  from  spin  waves.  Further¬ 
more:  (1 )  the  applied  field  value  at  the  peak  of  the  solid  curve  is  equal 
to  (2)  there  is  a  small  observable  decrease  in  the  Bragg  angle  Qg 
due  to  the  decreasing  wavenumbers  of  the  spin  waves  in  the  non-uniform 
internal  field;  the  inserts  in  Fig.  5.  4  indicate  a  small  increase  of  the 
delay  time  of  the  diffracted  light  pulse  due  to  the  decreasing  group 
velocities  of  the  waves;  and  (3)  the  diffracted  signal  disappears  when 

o  c 

Figure  5,  5  and  Fig.  5.6  show  plots  of  diffracted  light  intensity  as  a 
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SHEAR  ELASTIC  WAVES  (l.3GHz) 
POSITION  NEAR  THE  CENTER 


600  620  640  660  680  700  720  740  760  780 

EXTERNAL  MAGNETIC  FIELD  (Oe) 

Fig.  5.  4  Diffracted  light  intensity  (in  arbitrary  units)  of  shear 


elastic/spin  waves  as  a  function  of  external  magnetic 
field;  the  light  beam  is  incident  near  the  mid-point 
of  the  long  axis.  The  insert  is  the  scope  trace  swept 
from  right  to  left  when  b:  H®,  indicating  diffracted 
light  intensity  (marked  "spin"  on  the  curve)  and  a 
slight  decrease  of  the  group  velocity. 
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functlon  of  the  external  applied  field.  The  solid  and  dashed  curves  in 
both  figures  correspond,  respectively,  to  the  waves  propagating  parallel 
and  antiparallel  to  the  applied  field.  The  curves  in  each  figure  were 
obtained  by  reversing  the  applied  field  only  and  leaving  all  other  experi¬ 
mental  conditions,  including  the  sample  location  in  the  air  gap  of  the 
magnet,  unchanged.  Notice  in  these  figures  that  the  shape  of  the 
curves  is  the  same  for  a  given  incident  light  polarization  and  diffracted 
light  polarization  and  frequency  shift.  Thus,  light  diffraction  from  spin 
waves  is  not  affected  by  the  direction  of  the  applied  field,  in  agreement 
with  our  theory.  However,  the  magnitude  of  the  diffracted  signal  or 
the  height  of  the  curves  in  the  same  figure  differs  slightly.  Since  the 
applied  field  was  the  only  condition  that  had  been  changed,  the  conclusion 
is  that  the  sample  in  the  air  gap  experienced  different  applied  field 
homogeneity  for  the  two  opposite  directions  of  the  applied  field. 

Finally,  from  all  these  observations  we  conclude  that  the  intensity 
of  light  diffracted  from  coherent  spin  waves  depends  on  the  incident 
light  polarization  and  frequency  shift  of  the  diffracted  light.  This 
characteristic  is  intrinsic  for  the  interaction  of  light  and  spin  waves 
and  is  not  due  to  the  interference  between  the  light  diffracted  from 
elastic  waves  and  spin  waves. 


KATVfilio  k«MT  mttmttti 


-IM- 


•OOttO««OMOMO  700rt0740rt0  7*0 

tmiMAL  IMilWTiC  nCLO  CO*} 


Fif .  5.  ft  Diffracted  light  inlenalty  of  ahear  cLaatic/apIn 
aavta  aa  a  ftinctloo  of  external  mafnettc  ftelda, 
iUuatratlng  the  effect  of  homogeneity  of  the 
applied  field. 


600  620  640  660  680  700  720  740  760  780 

EXTERNAL  MAGNETIC  FIELD  (Oe) 

Fi^.  5.6  Diffracted  light  intensity  of  shear  elastic /spin 
waves  as  a  function  of  external  magnetic  fields, 
illustrating  the  effect  of  homogeneity  of  the 
applied  field. 


-132- 


78 

•*>.  3  Strong  Infrared  Bragg  Scattering  from  Coherent  Spin  Waves 

In  this  section  we  report  the  first  direct  observation  of  strong  Bragg- 
scattered  infrared  light  from  coherent  spin  waves.  The  maximum 
scattered  light  intensity  is  found  to  be  at  least  five  times  stronger  than 
that  from  longitudinal  elastic  waves  of  comparable  power.  The 
coherent  spin  waves  were  generated  through  efficient  space -gradient 
conversion  from  longitudinal  elastic  waves;  the  latter  were  excited  by 
the  sputtered  ZnO  transducer.  The  scattered  light  intensity  is,  again, 
found  to  be  dependent  on  the  incident  light  polarization  and  the  sign  of 
frequency  shift  (Stokes  or  anti-Stokes  lines),  and  independent  of  the 
direction  of  the  applied  magnetic  field.  All  these  observations  are  in 
agreement  with  our  ne^r  theoretical  results  derived  in  section  4.4, 
Example  2, 

Before  outlining  our  experiments,  it  is  deemed  helpful  to  review  the 
underlying  physics  of  microwave  spin/elastic  wave  conversion. 

A  ferromagnet  subject  to  a  spatially  varying  internal  magnetic  field 

has  a  non-uniform  refractive  index  insofar  as  spin  wave  propagation  is 

11  21 

concerned.  As  discussed  in  detail  by  numerous  authors,  '  a  spin 
wave  traveling  in  such  a  medium  will  undergo  variation  of  its  wave- 
number.  In  crossover  regions  where  the  frequency  and  wavenumber  of 
the  spin  wave  match  those  of  an  appropriate  elastic  wave,  conversion  to 
the  elastic  wave  or  vice  versa  will  occur  with  an  efficiency  tj  given 
approximately  by:^^' 


ri 


cr 


(3.6) 


The  quantity  H'(z^j,)  is  the  value  of  the  magnetic  field  gradient  in  the 
crossover  region;  h'  ..  is  a  critical  field  gradient  that  is,  in  simple 
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cases,  proportional  tc  the  square  of  the  magnetoelastic  splitting  frequency. 

Longitudinal  elastic  v/aves  launched  at  an  end  surface  of  an  axially 
magnetized  bar,  experience  an  increasing  magnetic  field  as  they  approach 
the  center  of  the  bar.  25  jf  appropriate  external  magnetic  field 

is  applied  such  that  there  are  crossover  regions  inside  the  sample, 
conversion  to  spin  waves  will  occur  at  those  regions,  with  efficiencies 
being  determined  by  Eq.  (3.  6).  After  conversion,  the  wavenumbers 
start  to  decrease  as  the  spin  waves  encounter  the  increasing  magnetic 
fields.  Subsequent  conversion  from  spin  waves  into  shear  elastic  waves 
is  therefore  not  possible.  On  the  other  hand,  if  one  starts  from  shear 
elastic  waves,  subsequent  conversion  from  spin  waves  into  longitudinal 
elastic  waves  would  be  possible  (Chapter  HI). 

For  YIG  at  room  temperature  and  with  the  magnetization  lying  in  a 
{100}  plane  at  an  angle  of  (22.  5°)  from  a  (100)  axis,  the  values  of 
for  the  longitudinal  and  shear  crossover  regions  are,  respectively, 

150  Oe/cm/GHz  and  19,000  Oe/cm/GHz.^^  Owing  to  the  small  longi¬ 
tudinal  magnetoelastic  splitting,  the  transition  from  longitudinal  elastic 
to  magnetic  character  is  very  abrupt.  Therefore,  in  a  fairly  uniform 
magnetic  field  the  spin  waves  can  be  locally  excited  through  the  spatial 
conversion  and  essentially  frozen  with  a  well  defined  wavenumber, 
before  they  eventually  spread  into  spin  waves  with  dissimilar  wave- 
numbers. 

Our  experiments  with  longitudinal  elastic/spin  wave  scattering 
utilized  the  geometry  shown  in  Fig.  5.  la.  The  rectangular  single 
crystal  YIG  bar  measures  4x4x12  (mm)  and  was  cut  so  that  its  long 
axis  is  oriented  in  a  {100}  plane  at  an  angle  of  (22.  5^^)  from  a  (100)  axis 
(Fig.  4.5).  Two  end  surfaces,  as  well  as  two  opposite  long  {1 00} 


-134- 


.surfac’t?s,  of  the  bar  were  polished  optically  flat.  The  sputtered  ZnO 
transducer  w’*!.  one-way  conversion  loss  of  approximately  15  db  at 
1 .  5  GHz  was  used  to  excite  longitudinal  elastic  waves.  The  peak  power 
of  a  0.  5  ^sec,  1 .  5  GHz  microwave  pulse  used  was  between  10  mW  and 
200  mW.  The  half-wave  plate  and  Glan-Thompson  polarizers  were  used 
to  obtain  various  combinations  of  incident  and  diffracted  light  polariz¬ 
ation.  The  Bragg  angle  for  the  longitudinal  elastic  waves  outside  the 
sample  was  6.  9°. 

Initially,  experiments  (Fig.  5.7)  were  done  at  fairly  high  magnetic 
fields  (2700  Oe)  where  at  1. 5  GHz  the  magnetoelastic  coupling  is 
virtually  zero.  Notice  that  as  the  long  axis  of  the  sample  is  not  along 
any  principal  crystallographic  direction,  the  photoelastic  parameters 
involved  are  a  linear  combination  of  the  photoelastic  constants,  and  the 
polarization  of  the  light  diffracted  from  the  longitudinal  elastic  wave  will 
not  be  in  the  same  direction  as  that  of  the  incident  light.  The  electric 

H  f\ 

fields  of  the  diffracted  light  (E  ,  E„)  are  related  to  those  of  the  incident 

V  n 

light  (E^,  E^),  according  to  Eqs.  (4.  31b)  and  (4.  35),  by: 


e'^ 

-P, ,  +Pi  o'*'2P. . 

E^ 

V 

2 

11  12  44 

V 

-  -^c 

"  4  ^ 

^^11  ■’■^12'^^  ^44 

L 

•  « 

where  c  is  a  constant  that  can  be  calculated,  if  necessary,  from 

Eq.  (4. 12a).  The  absolute  values  of  the  photoelastic  constants  P..  for 

71 

YIG  are  reported  by  Dixon  and  Matthews  as: 


IP^jl  =  0.025 

IPjgl  =  0.073 

IP44I  =  0.041 
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elastic  wave  pulse  propagating  along  a  ”22.  5^”  YIG 
bar.  The  graph  was  recorded  with  the  signals  from 
the  PAR  boxcar  integrator.  The  insert  is  the  scope 
trace. 
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The  signs  of  the  photoelastic  constants  are  not  known.  From  the 
measurements  of  the  intensity  of  the  diffracted  light  with  an  incident 
light  polarized  either  horizontally  (in  the  plane  of  incidence)  or  virtically 
(perpendicular  to  the  plane  of  incidence),  we  obtain  (from  Eqs,  (4.  36) 
and  (5. 1 )); 

(P,l+3Pl2-2P^^)2:(3Pjj^.Pj2+2P^/:(-P,j+Pj2+2P^/=,  55:4.1 

If  one  assumes:  P,  =  25:75:(-41),  then 

11  12  44  ’ 

(Pi1+3Pi2-2P^4)2:(3Pii+Pi2+2P4^)2:(-Pji+Pj2+2P4/  =  91:4:1 
On  the  other  hand,  other  possibilities  exist  such  as; 

PirPl2-P44  =  25:(-73):41 
then 

(Pi1+3Pi2-2P44)2:(3Pii+Pi2+2P44)^(P,i+P,2+2P44)2  =  300:27.5:1 

orifPii:Pj2:P44  =  25:73:41 

then 

(Pi1+3Pi2-2P44)^(3Pjj+Pj2+2P44)-(-Pii+Pi2+2P44)^  =  1.56:3,0:1 

or  if  Pii:Pi2:P44  =  -25:73:41. 
then 

<Pll+3Pi2-2P44>'^<3Pll+Pl2*2P44'^(-Pn^Pl2*2P^/  =  2:1:5 
From  these  results  we  conclud  that  it  is  most  likely  that: 

^ir^l2'^44  "  (0.025):(0.073):(-0.041). 

We  do  not  expect  our  measurements  to  be  in  good  agreement  with 
Dixon's  since  the  total  error  in  the  measurements  tends  to  become 
larger  when  linear  combinations  of  the  measured  quantities  are  formed. 
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Nevertheless,  in  all  cases  using  either  horizontally  or  vertically 
polarized  incident  light,  the  component  of  the  scattered  light  (either 
downshifted  or  upshifted  in  frequency)  with  polarization  perpendicular  to 
that  of  the  incident  light  is  always  much  smaller  than  the  one  with  the 
same  polarization  as  that  of  the  incident  light. 

Next,  the  magnetic  field  was  slowly  decreased  until  the  longitudinal 
crossover  region  occurred  near  the  middle  of  the  bar.  Under  these 
conditions  H'(z^^)  is  small  and  efficient  longitudinal  elastic  wave/spin 
wave  conversion  occurred,  as  indicated  by  the  strong  attenuation  of  the 
elastic  wave  echoes  observed  with  the  microwave  circuitry  (Fig.  5.  8(b)). 
When  the  laser  was  focused  on  the  crossover  region,  high  intensity 
diffracted  light  was  observed  to  coincide  with  the  time  of  arrival  of  the 
longitudinal  elastic  waves  at  the  crossover  point,  as  is  shown  in  the 
lower  trace  of  Fig.  5.  8(c). 

As  observed  in  the  experiments,  the  polarization  of  the  diffracted 
light  was  perpendicular  to  that  of  the  incident  light  only  when  the  latter 
was  either  horizontally  or  vertically  polarized.  This  indicates  that 
Ag  ^  0  in  Eq.  (4.47).  There  were  no  observable  changes  of  the  Bragg 
angle  from  that  of  the  longitudinal  elastic  waves.  This  excludes  the 
possibility  that  the  scattering  was  due  to  shear  elastic  waves,  for 
otherwise  the  Bragg  angle  would  be  almost  twice  as  large.  Further¬ 
more,  the  downshifted  diffracted  light  intensity  was  maximum  when  the 
incident  light  was  horizontally  polarized  and  minimum  when  the  incident 
light  was  vertically  polarized,  the  ratio  of  maximum  to  minimum 
diffracted  light  intensities  being  approximately  8.  For  the  upshifted 
case,  the  reverse  was  true.  Again  this  kind  of  asymmetry  cannot  be 
expected  from  shear  elastic  waves,  but  can  be  easily  explained  in  terms 
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(c)  Ho  =  736  Oe 

HORIZONTAL  TRACE  I SEC/CM 

Fig.  5.  8  Scope  traces  indicated  with  signals  from  the  micro- 

wave  circuitry  and  the  photomultiplier  as  a  function 

of  external  magnetic  field.  The  upper  and  lower 

traces  are  synchronized;  the  upper  traces  are  from 

the  microwave  receiver  and  the  lower  traces  are 

from  the  photomultiplier.  The  diffracted  light 

signal  is  polarized  at  ir/2  with  respect  to  the  incident 
i  d 

light  (Ejj,  ^v-^‘  upper  trace  indicates  the 

first  leakage  pulse,  and  two  longitudinal  elastic 
echoes  where  amplitudes  saturate  the  receiver.  In 
(b),  the  longitudinal  echoes  are  strongly  attenuated 
due  to  the  magnetoelastic  coupling.  In  (c),  the 
longitudinal  echoes  are  still  attenuated;  the  lower 
trace  now  indicates  a  strong  diffracted  light  signal 
from  coherent  spin  waves  (compare  the  diffracted 
light  pulse  with  that  in  Fig.  5.  7). 
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of  Eq.  (4.47)  if  one  assumes  that  f  and  (gj  j  ■6i2^^®44^  same 

sign  at  1150  nm.  When  the  magnetic  field  was  reversed,  all  the  above 
situations  remained  unchanged  in  agreement  with  Eq.  (4.  47). 

Figure  5.  9  is  a  typical  plot  (in  arbitrary  units)  of  the  intensity  of 
the  light  diffracted  from  spin  waves  and  elastic  waves  as  a  function  of 
the  external  magnetic  field;  the  intensity  of  the  light  diffracted  from 
both  waves  was  proportional  to  the  microwave  input  power  within  the 
power  range  (10  to  200  mW)  used  in  the  experiment  (Fig.  5. 10).  The 
incident  light  here  was  vertically  polarized,  the  light  spot  position  was 
approximately  in  the  middle  of  the  bar,  and  the  diffracted  light  was 
upshifted  in  frequency.  The  polarization  of  the  light  diffracted  from  the 
spin  waves  was  horizontal;  from  the  elastic  waves,  vertical.  There 
was  negligible  contribution  from  longitudinal  elastic  waves  to  the 
horizontally  polarized  diffracted  light  intensity,  as  supported  both  by  the 
calculations  and  experiments  done  at  higher  magnetic  fields. 

Specifically  in  this  case,  the  polarization  of  the  incident  light  and  the 
frequency  shift  of  the  diffracted  light  were  arranged  so  as  to  give  the 
strongest  possible  diffracted  light  intensity,  not  only  from  spin  waves 
(Eq.  (4.47)),  but  also  from  longitudinal  elastic  waves  (Eq.  (5.1)).  Note 
that  the  scale  has  been  amplified  by  a  factor  of  5  in  the  case  of  much 
weaker  vertically  polarized  diffracted  light  from  the  elastic  waves. 

Also,  the  minimum  and  maximum  of  these  two  diffracted  light 
intensities  do  not  coincide  due  to  the  very  small  conversion  region  and 
the  comparatively  large  size  of  the  light  spot.  The  minimum  of  the 
vertically  polarized  diffracted  light  is  not  exactly  zero,  which  may  be 
due  to  the  non-zero  magnetic  field  gradient  at  the  crossover  point 
caused  either  by  a  large  chip  in  the  sample  at  the  opposite  end  of  the 


HO 


A1ISN31NI  IHOn  a3d311V0S 


O) 

m 


bo 


I 


141 


5. 10  OiffractiMi  li^iht  Intrmilty  mm  a  function  of  input 
m*crowav<»  power. 


li'JinHiliifiT  or  I»y  lh«*  iraoKvorno  flelJ  ^radirntH.  The  trailinf*  edge  of 

Ihr  curve  itibeih'd  "ioi)|<itudtoAl  eUKtic  huvoa"  of  Kl^.  •*».  0  Ir,  of  courRc, 

Uue  lo  rcHidufli  convepRion*®  to  Rpln  uiive  form  ar  the  eroRRover  polnl 

moveR  tounrd  the  end  of  the  RAmple.  In  thiR  chrc,  the  Rpin  wnvcR  had 

nirendy  Ix'en  cut  off  and  did  not  diffract  any  ii^ht.  ARRutning  fuii 

converRion  from  elARtic  wavcR  into  Rpin  wavor,  and  URing  the  relative 

diffrat  te(i  intenRity  of  einRtic  wavcr  and  Rpin  wavea,  with  either 

horizontally  polarized  or  verticaily  polarized  incident  ii/ht,  we  obtain, 

at  \  ■  1 1  f»0  nm,  either  f  •  1 . 08  x  lO”*,  J(f;i  |  ‘R] 2. 1 5  x  lO’**  or 

vice  veraa  (bt'caufle  the  Rign  of  jr—r^ — --"ti  ■  %could  not  be  determined 

•»!  -<8||-R|2 

front  the  intenRity  data).  ItcRoiution  of  thiR  ambiituity  and  the  values  of 
niuRt  await  precise  magneto -optic  iiteaRurenients  of  f  and 
^*^1 1  I'ftraday  rotation  and  Voij^  effect.  We  also 

observed  the  intensity  of  light  diffracted  from  coherent  spin  waves  as  a 
function  of  the  light  position  along  the  bar  and  external  applied  field 
strength.  For  spin  waves  of  a  fixed  frequency,  the  strength  of  the 
external  applied  fields  has  to  be  of  such  a  value  that  the  light  beam 
position  coincides  with  the  longitudinal  crossover  point.  Therefore  in 
principle  one  has  a  good  way  of  measuring  the  internal  field  profile  by 
using  light  scattering  from  coherent  spin  waves.  Of  course  if  the  field 
gradient  at  the  crossover  point  is  large,  weak  conversion  from  elastic 
waves  to  spin  waves  occurs  and  the  scattering  intensity  is  reduced,  as 
was  indeed  observed.  Since  there  is  a  large  chip  on  the  sample, 
theoretical  calculations  of  the  internal  field  profile  are  not  feasible. 
Therefore  no  results  on  internal  field  profile  probing  are  presented  here. 
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5,  4  Measurements  of  Temporal  Longitudinal  Elastic/Spin 

Wave  Conversion  Efficiencies 

Optical  probing  of  longitudinal  elastic/spin  waves  subjected  to  pulsed 
magnetic  fields  was  used  to  measure  the  temporal  longitudinal  elastic/ 
spin  wave  conversion  efficiencies  by  directly  observing  the  scattered 
light  intensity  from  the  converted  spin  waves.  It  is  difficult  to  make 
the  measurements  with  conventional  microwave  techniques  whenever 
intrinsic  non-uniform  fields  due  to  demagnetizing  effects  occur  inside  a 
ferrimagnetic  sample.  The  experiments  were  performed  under  identical 
conditions  and  using  the  identical  sample  as  in  section  5.  3.  The  sample 
used  here,  however,  has  a  38  Oe/Amp  coil  uniformly  wound  around  it, 
except  for  a  space  near  the  center  of  the  long  axis  where  the  light  beam 
was  focused.  The  lower  trace  of  Fig.  5.11(a)  shows  the  light  pulse 
diffracted  from  coherent  spin  waves  (H^  =  726  Oe).  The  upper  traces 
(synchronous  with  the  lower  traces)  in  Fig.  5. 10  show  a  positive  field 
pulse  that  is  in  the  direction  of  the  applied  field.  In  Fig.  5. 11(a)  the 
small  field  pulse  is  applied  after  the  spin  waves  cross  the  light  beam 
and  therefore  has  no  effect  on  the  diffracted  light.  In  Fig.  5. 11(b)  the 
dc  ’^“'as  field  is  decreased  by  15  Oe  and  the  field  pulse  is  increased  by 
15  Oe  and  is  applied  before  the  elastic  waves  are  launched  at  the  bar  end. 
The  result  is  that  the  elastic  waves  still  convert  to  spin  waves  at  the 
light  beam  position  through  a  spatial  gradient  and  therefore  a  diffracted 
light  pulse  from  coherent  spin  waves  can  be  seen,  as  in  5.  11(a).  Now 
the  field  pulse  is  increased  by  another  23  Oe  while  still  being  applied  at 
the  same  time,  as  in  Fig.  5. 11(b).  Therefore  no  spin  waves  could  be 
observed  at  the  original  light  beam  position,  as  shown  in  the  lower 
trace  of  Fig.  5. 11(c).  We  then  changed  the  timing  of  the  field  pulse 


Diffracted  light  intensity  (E  ,  E^. )  from  coherent 
spin  waves  as  a  function  of  the  timing  of  the 
pulsed  bias  field;  the  spin  waves  are  generated 
through  temporal  conversion  from  longitudinal 
elastic  waves  (time  scale  0.  5  ^ec/cm),  A  full 
explanation  is  given  in  the  text. 
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while  maintaining  its  amplitude.  In  Figs.  5. 11(d)  through  5. 11(g),  the 
diffracted  light  signal  was  observed  at  a  fixed  delay  following  the 
application  of  the  field  pulse.  Figure  5, 12  is  an  expanded  version  of 
Fig.  5.11. 

These  two  figures  not  only  give  further  endorsement  to  the  interpre¬ 
tation  of  the  light  scattering  from  coherent  spin  waves,  as  discussed  in 
the  last  section,  but  also  enable  us  to  measure  the  efficiencies  of 
temporal  spin/longitudinal  elastic  wave  conversion.  The  results  are 
indicated  in  Fig.  5. 13  where  the  field  gradients  are  determined  from  the 
risetime  and  amplitudes  of  the  field  pulse.  These  results  do  not  agree 
very  well  with  the  theory.  The  discrepancy  may  be  caused  by  the 
non-uniform  internal  fields  inside  the  sample. 
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(f) 

Fig.  5. 12  (a)  diffracted  light  pulses  from  longitudinal 

elastic  waves  (E^,  {time  scale  0.  5^sec/cm). 
(b)  to  (f)  -  an  expanded  version  of  Fig.  5. 11 
(E^,  )  (time  scale  0.  2  usec/cm). 
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Appendix  I 

Effective  Anisotropy  Fields  and  Small  Anisotropy 
Tensor  N  for  Cubic  Crystals 


We  here  discuss  the  effective  anisotropy  fields  and  small  signal 
anisotropy  tensor  N  for  cubic  crystals  in  a  coordinate  system  in  which 
the  static  magnetization  M  is  along  the  z-axis  (referred  hereon  to  the 
transformed  coordinate  system),  but  otherwise  arbitrary  with  respect  to 
the  crystallographic  directions. 

In  a  coordinate  system  (^,  rj,  {)  with  axes  along  the  cube  edges,  the 
anisotropy  energy  density  of  a  cubic  crystal  can  be  written  as; 


W 


an 


3 

L 

i,  m  =  l 
i  >  m 


S.  m 


mJ 

i  ri  C 


(Al.  1) 


As  we  are  only  concerned  with  the  effective  anisotropy  field  portion 
Sri 

H.  which  is  at  most  linear  in  the  transverse  small  signal  magnetization 
components  in  the  transformed  coordinate  system,  '  the  term  of  W 
containing  K2  in  Eq.  (Al.  1)  can  be  safely  neglected.  Equation  (Al.  1) 
can  then  be  written  as; 
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an 
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(A1.2) 


Suppose  the  transformed  coordinate  system  is  oriented  at  the  Euler 
Angles  ((p,  6,  0)  from  the  "cubic-edge"  coordinate  system.  Expressing 
in  terms  of  the  transformed  coordinate  system,  one  obtains; 

K.  K,  _  3  3 

^an  =  T - ^ 

2M  ‘-je  =  lf,  i,j,  k=l 


r-t  rp  rp 

"if 


T„,  M.M.M.M 
ik  f  1  j 


k. 


(Al.  3) 
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where  ^  is  given  in  Eq.  (2.  4).  The  effective  anisotropy  field  in  the 
transformed  coordinate  system  (where  the  static  magnetization  M  is 
along  the  z-axis)  can  then  be  written  down  as; 


aw 


an 


2K, 


3 

£ 


3 

£ 


i=l  f,  j.  k=l 


^^ir^  33  ^  2  2 

— 4[/^i  T, .  T,3  M3  +  3  £  ^  T,. 


(AI.4) 
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je=i 


+  o(m2.  ml  m^m^) 


j?=i 


Only  terms  which  are  at  most  linear  in  M  and  M  are  retained  in 

X  y 

2  2 

Eq.  (Al.  4)  and  higher  order  terms  like  0(Mj,  Mg,  M^Mg)  are  to  be 

3.0  30 

dropped.  Furthermore,  in  the  expressions  H.  and  H„  (i  =  1,2  in 

2^1  2  3  3 

Eq.  (A1.4)),  the  term - 2^  T-.  T-o  M,  (i  =  1,  2)  does  not  contribute 

-an 

to  h  in  the  small  signal  equations  of  motion  Eq.  (1.2).  On  the  other 
2^1  3  4  3  an 

hand,  the  term - 2  S  T-„M„  in  H„  is  the  only  one  that  enters  the 

u  M^  je=l 

equations  of  motion  and  is  in  fact  equal  to  in  Eq.  (1 . 3a).  We 
therefore  write  (using  M3  —  ivl,  M^  =  m^,  and  Mg  =  nig): 
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Recall  the  definition  of  the  small  signal  anisotropy  term  N(Eq.  (1. 3b) 


as; 


hf^  =  L  -Nf"  m.  (i  =  1,  2) 
^  j=l  '  y  J 


(A1.6) 


Using  (Al.  5),  (Al.  6)  and  the  T  from  Eq.  (2.  4)  (with  0  =  0),  one  obtains: 

2K. 
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an 
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=  + - ^  sin^  2  0  (1  -  ^sin^  2(p)j  (Al.7a) 
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(A1.7b) 
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^  ^  sin© sin  20 sin  4(p]  (Al.7c) 


2K 

Hg"  =  t  2e  +  sin'^0sin^2(p)J  (Al .  8) 


**  stn 

Since  in  the  small  signal  equations  of  motion,  N  is  always 

£tri 

multiplied  with  the  magnetization  M  and  Hg  /M,  with  m^  or  m^.  In 

view  of  this,  sometimes  it  is  found  convenient  to  define  an  effective  N 
siri  t*  sin 

such  that  Hg  /M  and  N  are  combined  as  follows; 
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(A1.9a) 


Ran 

^^22  M  *^^22 


!K 

— 1^1  ■  i(sin^20+  sin^0sin^2(p)J 


u  M 

“n 


3  •  2„  .  2„  \ 

-  ^sin  0  sin  2(pj 


(Al.9b) 


N 


an 

12 


=  N, 


an 

21 


=  N 


an 

12 


=  n: 


an 

21 


■1  51- 
2K 


■^2^- I  sin  6  sin  20  sin  (Al.9c) 


Using  this  effective  anisotropy  tensor,  N  ,  the  internal  magnetic 
field  H.  in  the  small  signal  equations  of  motion  is  just  the  sum  of  the 
applied  and  demagnetizing  fields  only. 
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Appendix  II 

The  First  Order  MsijnetoeUstic  ConKtnnts  for  an  Arbitrarily 
Oriented  Coordinate  Syatem  in  Cubic  Cryataln 


The  first  order  (in  strain)  maRnetoelaatic  constants, 
pertaining  to  the  small  signal  magnctoelastic  equations  of  motion  in  an 
arbitrarily  oriented  coordinate  system  in  cubic  crystals,  are  derived 
here.  The  connection  between  the  magnetoelastic  constants  b|j|^ 
defined  in  Eq.  (1. 3d)  aiid  ifl  nlso  discussed. 

In  general,  the  first  order  (in  strain)  magnetoelastic  energy 
density  of  a  cubic  ferromagnet  can  be  written  in  terms  of  the  first 
order  magnetoelnstic  constants  *>8; 


W 


me 


(f,  i,  J  and  k  ■  I.  2,  3) 


(A2.1) 


where  Sjj^  is  the  symmetrical  finite  strain  tensor  and  therefore  the 
magnetoelastic  constants  have  the  symmetry  properties; 


^’fijk  *  ^Ifjk  “  ‘’fikj  '  ‘’ifkj 
With  cubic  symmetry,  the  b^^j^'s  have  the  additional  symmetry 
properties; 


(A2.2a) 


Tijk 


’jkfi 


(A2.2b) 


Consequently  there  are  at  most  18  independent  b^jji^'s  for  an  arbitrarily 
oriented  coordinate  system  in  cubic  crystals. 

Since  we  are  mainly  concerned  here  with  the  small  signal 
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I'qtintions  of  motion,  the  finite  strain  tensor  S..  can  be  approximated  by 
the  small  signal  strain  tensor  and  can  be  written  as; 


’Jk 


,  5e!,\ 

2  JtXj/ 


(A2.3) 


where  p  is  the  clastic  displacement. 

Suppose  the  transformed  coordinate  system  is  oriented  at  the  Euler 
Angles  (g,  0,  0)  with  respect  to  the  "cubic-edge"  coordinate  system.  We 
now  want  to  |)rocecd  and  find  the  terms  of  the  ordinary  first 

order  magnctoelastic  constants  b^  and  b2  And  the  Euler  Angles. 

If  W^^^p  is  expressed  in  the  "cubic-edge"  coordinate  system  there 
arc  only  two  different  i.  c. ,  bj  and  bj  since 


^1  *  ^iii 


(i  '  1.  2,  3) 


(A2.  4r) 


2^2  *  ^jij  '  ^jiji  *  %ji  *  ^jiij  j.  2,  3)  (A2.4b) 


All  other  b 


fijk 


(A2. 4e) 


We  can  then  rewrite  (A2. 1)  as; 

3  „  3 


W 


me 


„2  1  I  «  ,  2  1  m  ImJ 


<r 


jf  m 

^2/^2  ^  N 

if  m 

(b,  -  b,)  3  2 

r  Mjs 

1=1  *  ** 


(A2.  5n) 


(A2.  5b) 


Since  the  first  term  of  Eq.  (A2.  5b)  is  invariant  under  rotations,  only 
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the  second  term  has  to  be  transformed  into  the  new  coordinate  system. 
Let  uo  define  as  the  invariant  part  and  as  the  transformed  part 
of  the  in  the  new  coordinate  system. 


'nil 


=  b 


2222 


=  b' 
3333 


b, 

_ 4 

M 


(A2.  6a) 


and 


*^1212  ^^1221'  ^^^12  *^2121^ 
^1313^*^1331'  ^3113*  ^31 3p 

”  *^^23  ^*^2332'  *^^23'  *^^32^ 

=  A_ 

2M^ 


all  other  (b')'s 


(b.  -  b„)  3 

1.11  ^  i ^  ^  'T  'T  'T 

°fijk  “  o  E  T,,  T„,  T, 


M 


2  "j  "jff  "jfk 


(A2.6b) 


(A2.6C) 

(A2.7) 


where  T  is  given  in  Eq.  (2.4).  Note  that  from  Eq.  (2.  7),  the  b'^^jj^'s 
are  invariant  with  the  interchange  of  all  four  indices  and  thus  at  most 
take  15  different  values.  The  first  order  magnetoelastic  constants 
bfijk  can  then  be  easily  obtained  from  Eqs.  (A2.  6)  and  (A2.  7).  The 
results  are  indicated  in  Table  A2. 1. 

For  small  signal  magnetoelastic  problems,  b^jjj^'s  are  not  very 
convenient  to  use.  Now  let  us  choose  a  particular  sot  of  the  Euler 
Angles  itp,  6,0)  such  that  the  z-axis  of  the  transformed  coordinate  system 


coincides  with  the  static  magnetization  M,  i.e.,  M=  i^  M^.  Then 
working  in  this  particular  coordinate  system,  we  note  from  Eqs.  (1.  3d), 
(1. 5)  and  (A2. 1),  that  the  only  energy  terms  that  contribute  to  the 


The  First  Order  Magnetoelastic  Constants  i"  a  Coordinate  System 

at  the  Euler  Angles  (<p,  $,  0)  from  the  "Cubic-Edge"  Coordinate  System 


Table  A2. 1  (continued) 


2221  2212  2122  1222 


Table  A2.  ]  (continued) 
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effective  magnetoelastic  torque  are  those  with  one  and  only  one  Mg. 
Consequently  we  may  define  a  small  signal  magnetoelastic  energy  density 
w  with  the  corresponding  small  signal  magnetoelastic  constants  b...  as 

1X10 

sp. 

w__  =  ^  E  b,„  M,  ^  (A2.8) 


me 


1  3 

^  E  b...  M.  ^ — 


Rewrite  (A2. 1 )  as; 


W. 


me 


M,  3  p 

E  [boo...  M.,  +  2 


M 


2.  “  jL“33jk  ^’^3 


^  l  1  /^Pi  +  ®Pk\ 

E  b„..,  M.  1  •  T 
j  j  3i]k  ij  2  Sx^/ 


2  3 

L  E  b...,  M.  M.  S., 

i.f=lj.k-l  ^  ^ 

M  3  2  ap 


(A2.9) 


1 


2  3 


+  E  E  bfiikM^M  S 

M^  i,  f=lj,k=l  ^  ^ 


By  comparing  (A2.  8)  and  (A2.9),  and  letting  Mg  —  M,  one  obtains 


'’Sjk  °  '’33jk 

V  =  ”  (i  =  l.  2;K  =  >.2.3) 

and 


(A2.  10) 


b...  =  b.,  , 

ijk  ikj 


(A2.  11) 


According  to  Eq.  (A2.11),  there  are  in  general  1  8  different  b.  .,'s 

IJK 

with  cubic  symmetry,  at  most  12  independent  b... 's.  These  are  the 

ijk 

ones  corresponding  to  the  first  12  in  Table  A2. 1.  Using  Eq.  (A2.  10) 
and  Table  A2. 1,  we  obtain  Table  A2.  2. 


Table  A2.2  (continued) 
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Finally,  using  Eqs.  (A2.7),  (A2. 10),  and  (2.  5),  we  want  to  express 
specifically  and  in  a  coordinate  system  which  lies  at  the  Euler 
Angles  (tp,  0,  Ip)  with  respect  to  the  "cubic-edge"  coordinate  system; 

^^l”^2^  f  r  2  2  "1 

b233  =  — sin  0l- I  sin 26' (cos  20  +  ^  sin  ©sin  2<p)J 


+  cos 


sin^0  sin  4<pj  j- 


(A2.12) 


(bi-b)f  r  3  n 

=  - — -isin^j-^sin  0  sin  4<p  1 

-  cos  0  ^  -  I  sin  20  (cos  20-»-  |  sin^0sin^  2<p)J  | 


(A2.  13) 


It  may  be  interesting  to  note  from  Eqs.  (A2. 12),  (A2. 13),  and 

2  2 

Table  A2.  2  that  the  term  (bj^33  +  b233)  is  independent  of  ip  and  can  be 


written  as 


.2  2  _ 

^133  ^233  "  |_4Sin 


r  2  2  2 

l^sin  40+ sin  2<p(7  sin  0sin20sin  40 


(A2. 14) 


2  4  2  6  2  *1 

+  I  sin  2(p  sin  0  sin  2  0)  +  ^  sin  0  sin  4<pJ 


It  is  also  useful  to  express  some  b.jj^'s  in  the  present  coordinate 


system; 


‘’232  '  ‘’223  ° 


+  sin  2 sin  0  sin  20  sin  4(p] 


(A2.15) 


,  2  Ti  .2^  .2 

+  cos  ^  sin  0  sin 


2„]} 


162- 


^31  ~  ^113 


bg  b,  -b. 


^2  f  2  r  2  2  ” 

+  — — -l^sin  1  sin  0  sin  2(p 


r  1 

+  sin  2il)  I  -  ^  sin  0  sin  20  sin  4^J 
+  cos^^p  sin^  20  (1  -  sin^  2(p)  J  I 


'^123'  ^^132  ^  ^^213  ^^231 


(b^-b2) 

M 


sin0  sin  20  sin 

2  r  2  2  2  I 

+  sin  20  sin  0  \  sin  2(p(l  +  cos  0)  -  cos  0  I 

+  cos^^j^- IjSin©  sin  20sin  4^j  I 


(A2. 16) 


(A2.  17) 
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Appcndlx  111 


Dielectric  Tensor  of  a  '*Cublc"  Magnetic  Medium 


59 


As  pointed  out  by  many  authors"'^' and  discussed  briefly  in 
Chapter  IV,  magneto-optics  can  always  be  described,  macroscopically, 
in  t^rms  of  an  effective  dielectric  tensor  f  by  setting  the  permeability 
p  =  1.  Once  this  f  is  specified,  one  can  substitute  it  into  Maxwell's 
equations  and  proceed  to  solve  the  equations  with  suitable  boundary 
conditions. 

The  dielectric  tensor  7  is  defined,  for  the  electric  field  of  a  single 

frequency  a).  ,  as 

^  3 


(I  *1,2,  3)  (A3.1) 


In  a  non-polarized  medium,  is  symmetrical;  however,  in  a 

magnetically  saturated  medium  with  magnetization  M  the  generalized 

59 

principle  of  symmetry  of  kinetic  cojjfficient  requires: 


(A3.  2) 


Furthermore,  since  we  are  only  concerned  with  sinusoidally  varying 
fields,  we  can  use  complex  notation,  i.  e. ,  l5,  E  e  ^  .  Therefore  c 
is  now  complex.  Assuming  media  is  optically  lossless,  f  must  then  be 
Hermitian. 


where 


Cjj  =  cJ'j  (A3.  3a) 

€jj  =  cjj  +  jc’jj  (rt3.3b) 


Using  (A3. 1),  (A3.  2^  and  (A3. 3),  we  then  obtain; 


-164- 


€[.(M)  =  cJi(M)  =  fJjC-M) 


(A3.  4) 


e;'(M)  =  -c];(M)  =  -c;|j(-M) 


(A3.  5) 


Following  (A3.  4)  and  (A3.  5),  one  concludes  that  c!.  is  symmetrical 

..  ^ 

and  an  even  function  of  M,  while  c-  is  antisymmetrical  and  an  odd 
function  of  M. 

Using  antisymmetrical  unit  tensor  e... ,  c"  can  then  be  written  in 


terms  of  a  vector  3. 


=  e...  G, 
ij  ijk  k 


(A3.  6) 


1.  e. , 


•f  =  f 

zy  'yz 


*xz  '  ^zx 


=  G. 


X 


=  G. 


V 


(A3.  7) 


'*zy  ■  Vz 


=  G. 


D.  =  (cl  .  +  icV.)  E. 
1  ij  iJ  J 


=  c'. .E.  +  i  e. G,  E. 
ij  J  k  j 


(A3.  8) 


=  ,;.E:.c(ex3>j 

A  medium  in  which  the  relation  between  and  6  is  of  this  form  is 
said  to  be  gyrotropic. 

As  magneto-optical  effects  are  fairly  small,  one  can  expand  c  into 

power  series  of  M  or  direction  cosines  (a  »0  ,0_)  of  M 

X  V  2 


(A3.  9) 
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(A3. 10a) 


®ijkl 


g 


jikjf 


g 


ijik 


=  g 


jiik 


(A3. 10b) 


whore  and  gyj^jj  contain  all  the  symmetry  properties  of  a  crystal. 

For  example,  when  M  is  along  a  [00  *]  direction  in  a  cubic  crystal,  one 

can  simplify  (A3.  9)  by  keeping  only  the  linear  term  of  at  and  using  the 

z 

symmetry  properties  of  the  crystal.  One  obtains,  in  the  "cubic-edge" 
coord  unto  system  (or  principal  axes): 


-jfa. 


0  Cj  +  g^a^ 


(A3.  11) 


This  is  the  well  known  dielectric  tensor  whose  off-diagonal  terms  give 
rise  to  Faraday  rotation  and  diagonal  terms  magnetic  birefringent  effect 
(Voigt  effect). 

In  the  case  of  spin  wave  scattering  with  the  static  magnetization  M 

along  z-axis  direction,  one  should  keep  all  linear  terms  of  ot  and  a 

X  y 

which  may  also  come  from  second  order  or  higher  order  terms  such  as 

2  3 

a  a  anda  a  . 

X  Z  X  Z  X  z 

Notice  from  Eqs.  (A3.  9)  and  (A3. 10)  that  f.  is  a  2nd  rank  polar 

Kft 

59 

tensor  and  g^jj^jj  ^  4th  rank  polar  tensor.  For  a  cubic  crystal 
(class  m3m),  like  YIG,  using  the  "cubic  edge"  coordinate  system  one 
obtains 


(A3. 12) 
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where  6^^^  is  the  kronecker  delta  function,  and 


[g]  = 


11 

^12 

^12 

0 

0 

0 

12 

^11 

^12 

0 

0 

0 

12 

^12 

^11 

0 

0 

0 

0 

0 

0 

^44 

0 

0 

0 

0 

0 

0 

g44 

0 

0 

0 

0 

0 

0 

^4' 

(A3. 13) 


where  is  written  in  the  matrix  form  as  in  the  case  of  the  photo¬ 

elastic  tensor  in  Eq.  (4.25a).  The  higher  order  terms  which  contain 

2  3 

a  first  order  term  of  a  or  a  such  as  a  a  and  a  a  do  not  exist  in  a 

X  y  z  X  z  X 

cubic  crystal  and  therefore  need  not  be  considered.  In  a  crystal  with 

2  3 

less  symmetry,  where  terms  like  a„0(^  andCK  do  exist,  we  can 

Z  X  Z  X 

always  lump  them  into  effective  or  f. 

Combining  Eqs.  (A3.  6)  (A3.  9),  and  (A3. 10a),  the  small  change  of 
the  dielectric  tensor,  6c -  =  6c|j  +  3  'jj,  can  be  expressed  as: 


11 

^12 

Si  2 

0 

0 

0 

12 

Sll 

Si  2 

0 

0 

0 

12 

Si  2 

^11 

0 

0 

0 

0 

0 

0 

®44 

0 

0 

0 

0 

0 

0 

^44 

0 

0 

0 

0 

0 

0 

s^, 

0  a 

j 

0  a" 
3 

0  a! 


(A3. 14) 


where  6c'..  =  6c'.. 
1]  31 
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where  6c'.'.  =  -6c''-.  The  quantity  f  is  related  to  the  Faraday  effect, 

J  ^ 

while  the  matrix  is  related  to  the  birefringence  effect.  Besides, 

f  and  g  ,  in  general,  are  a  function  of  optical  angular  frequency  "^r 
m.n  Lj 

wavelength 

59 

The  derivation  here,  following  Landau,  is  similar  to  that  of  Auld 

and  Wilson,  °  but  the  result  is  quite  different.  They  failed  to  take 

account  of  the  terms  g^^^  in  Eq.  (A3. 14);  therefore  the  effect  of  on 

6Cxz  is  nil.  We  will  see  that  these  terms  are  essential  in  interpreting 

57 

our  experimental  results  as  well  as  Collins  and  Wilson's  "anomalous" 
experimental  observation.  Besides,  the  origin  of  g  comes  from 
birefringence.  In  an  experiment  like  coherent  spin  wave  light 
scattering,  by  looking  at  the  geometry  of  the  interaction  one  would 
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